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PREFACE 


This  report  describes  the  development  of  the  digital  com- 
puter real-time  program  for  the  mathematical  model  of  a JEFFB 
Amphibious  Assault  Landing  Craft  (ACV-LC) . The  mathematJcal 
molel  is  presented  in  tlio  report  NAVTIL^EQUIPCEN  73-C-0138-1. 

The  real-time  program  is  designed  to  be  implemented  in  the 
Sigma  7 computer  facility  of  the  Naval  Training  Equipment  Center. 
The  system  integrating  the  program  and  other  components  is  being 
used  as  an  ACV-LC  experimental  training  device. 

In  the  areas  of  Vehicle  Generated  Wave,  Offshore  Wave, 
and  the  compartment  pressure  models,  special  program  techniques, 
simplification  and  remodeling  of  the  mathematical  model  were 
made  in  order  to  make  the  simulation  run  in  real-time. 
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SECTION  1 
INTRODUCTION 

This  report  details  an  all  digital  real-time  simulation  developed  for  the  Naval 
Training  Equipment  Center (nav/TRAEOUIPCEN)  of  the  Bell  Aerospace  JEFF-B  design,  an  air 
cushioned  vehicle  landing  craft  (ACV-LC).  This  simulation  is  .intended  for  use 
as  a pilot  trainer;  therefore,  only  gross  vehicle  characteristics  and  performc  nee 
can  be  expected.  This  report  contains  a nresentation  of  the  model  and  a description 
of  tlie  implementation.  The  mathematical  model  used  has  been  detailed  in  the 
report;  Mathematical  Model  of  an  Air  Cushion  Vehicle. 

A complete  description  of  specific  model  areas  as  they  differ  is  included  for 
correlation  purposes  where  possible  and  necessary.  Some  of  the  areas  of 
modelling  completely  differ  from  the  given;  other  areas  differ  only  in  format 
or  form.  In  all  cases,  however,  the  symbology  has  been  changed  but  is  adequately 
defined. 

This  report  discusses  techniques  used  in  the  simulation  development.  In 
several  cases,  special  techniques  were  developed  to  make  the  simulation  run 
in  real-time. 

Complete  operational  procedures  are  given  in  detail  in  order  to  provide  the 
operator  with  all  the  tools  to  run  the  simulation  and  to  make  on-line  changes 
as  necessary  (i.  e.  initial  conditions). 

This  report  is  not  intended,  necessarily,  for  use  in  duplicating  the  simulation. 

Although,  a complete  description  of  all  components  is  provided,  no  programming 
techniques  and  methods,  as  such  have  been  included, 

•f; 

This  all  digital  simulation  was  developed  to  use  the  equipment  and  be  supported 
by  the  personnel  at  NAVTRAEQUIPCEN’s  general  purpose  computer  facility.  This  facility 
includes  a Xerox  Sigma  7 Computer;  complete  hybrid  interface  of  analog  to  digital 
converters,  digital  to  analog  converters, and  discrete  inputs;  aind  trvmking  station 
with  removable  patchboard.  Also  connected  to  the  Sigma  7 is  a line  drawing  CRT 
with  appropriate  interface.  Necessary  for  simulation  of  an  ACV-LC  is  the  mockup 
of  the  pilot  cabin  with  a set  of  controls  duplicating  placement  and  operating 
characteristics  as  are  located  on  the  vehicle. 

The  CRT  display  is  used  to  provide  the  visual  effects  necessary  to  dock  the 
ACV-LC  inside  the  stern  well  of  a mother  ship.  This  display  is  external  to  the 
cockpit  and  provides  a three  dimensional  picture  for  approach  and  docking 
maneuvers. 

SECTION  II 
ACV-LC  JEFF-B 

The  Bell  Aerospace  Air  Cushioned  Vehicle  Landing  Craft,  JEFF-B,  has  been 
modelled  and  simulated  real-time  in  a pilot  trainer.  There  are  two  very  important 
areas  in  understanding  the  operation  of  the  simulation  and  the  content  of  this  report. 

First,  the  size  and  shape  of  the  vehicle  and  type  of  effectors  determine  the  Aip's 
handling  characteristics.  Second,  the  type  of  controls  and  displays  available  to 
the  pilot  define  the  method  of  driving  the  ship. 
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VEHICLE  CONFIGURATION 

The  ACV-LC  is  virtually  rectangular  in  shape  with  a flat  bottom  and  a semi- 
rigid inflatable  skirt  surrounding  the  periphery  of  the  hull.  This  skirt  hangs 
below  the  hull  bottom  approximately  4.  5 ft.  allowing  adequate  space  for  an  air 
volume  or  mass  for  flotation  purposes.  Air  is  pumped  into  each  quarter  of 
this  cushion  compartment  by  large  fans,  providing  the  pressure  necess^  *y  for 
the  lift  to  keep  the  hull  above  tvater.  Leakage  of  air  occurs  at  the  boundary 
between  skirt  and  water  surface  and  is  vented  through  moveable  nozzles  at  the 
top  of  the  hull.  (See  Figure  1) . 

These  2 thrust  nozzles  (1  each,  port,  starboard)  are  located  in  the  forward 
part  of  the  ship  at  the  top  of  the  super  structure  and  provide  large  amounts 
of  thrust.  These  nozzles  are  capable  of  directing  this  thrust  through  360° 
and  move  in  unison. 

Located  at  the  stern  of  the  ACV  are  2 huge  ducted,  variable  pitch  propellers 
(1  each,  port,  starboard)  that  provide  thrust  in  the  longitudinal  direction. 

Each  propeller  is  driven  directly  from  its  output  power  shaft  which  is  governor 
controlled  and  coupled  to  3 gas  turbines.  The  air  inlet  fans  are  geared  to  the 
power  shaft. 

Located  behind  each  propeller  are  large  air  vanes  (rudders)  which  move 
togetlier  and  provide  steering  control.  When  used  in  conjunction  with  the  bow^ 
thrust  nozzles,  the  ship  can  make  coordinated  turns,  sharp  turns,  and  sloppy 
turns;  also  the  ship  can  be  moved  laterally. 

The  pilot  cabin  is  located  at  the  bow  on  the  starboard  side  at  the  top  of  the 
super  structure  allowing  good  peripheral  vision  from  the  ACV-LC.  The 
mother  ship  and  its  stern  well  can  be  seen  clearly  for  docking  purposes. 

The  air  pressure  under  the  hull  tends  to  flatten  the  water  surface  and 
to  create  a depression  when  hovering.  This  effect  causes  water  to  be  pushed 
outward  from  beneath  the  ship  causing  vehicle  generated  waves  (VGW). 

These  waves  are  of  little  consequence  to  the  vehicle  performance  except 
for  the  bow  wave,  if  moving  forward.  However,  the  pressure  distribution  is 
greatly  effected.  This  bow’wave  is  pushed  by  the  flat  noa  V hull  and  causes 
the  humping  speed  phenomenon.  This  bow  wave  causes  sufficient  drag  that  once 
created  at  constant  low  speeds  it  is  not  possible  to  accelerate  the  vehicle  through 
it.  This  pressure  wave  has  a i laximum  height  in  the  velocity  range  16-20  knots. 
At  slowly  increasing  forward  speed,  there  is  a region  of  decreased  acceleration; 
then,  as  the  vehicle  climbs  over  the  wave,  a sudden  increase  in  acceleration 
occurs  as  the  vehicle  rides  down  the  wave. 

However,  if  the  vehicle's  trajectory  is  such  that  it  crosses  its  own  path,  then 
interference  of  the  past  VGW  \vith  the  present  VGW  does  occur,  impacting  the 
vehicle  maneuverability. 
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PROPELLER  DUCT 


FIGURE  1.  ACV-LC  JEFF-B  Eff 
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PILOT  CONTROLS  AND  DlSPLi\YS 

The  pilot  needs  all  controls  to  maneuver  the  vehicle.  Since  the  vehicle 
has  no  automated  control  system,  the  individual  control  devices  for  all 
effectors  must  be  available  to  the  pilot.  (See  Table  1).  All  C gas 
turbine  throttles  (3  PORT,  3 STBD)  command  turbine  speed.  The  STBD 
and  PORT  governor  controls  set  the  output  power  shaft  speed  according 
to  the  command  and  the  related  gas  turbine  speeds.  The  power  shafts 
govern  the  propeller  speeds  (propulsion),  fan  speeds  (compartment  air 
flows),  and,  consequently,  the  tlirust  nozzle  air  flows.  There  are  3 
propeller  pitch  control  devices,  one  each  PORT  and  STBD  to  individually  set 
each  propeller  pitch,  and  a vernier  to  change  both  propeller  pitch  commands 
simultaneously  bj'  mo'/ing  the  steering  wheel  fwd/aft.  The  thrust  nozzles  are 
directed  by  the  steering  wheel  and  by  a svdtch  to  direct  the  thrust  fwd/aft. 

The  rudders  are  directed  by  foot  pedal  commands. 

For  forward  motion,  the  pilct  commands  equally  the  6 gas  turbines  and  the 
2 pov/er  shafts.  To  command  the  correct  forward  thrust  for  a given  speed 
the  pilot  directs  the  nozzle  fwd  or  aft  and  sets  the  two  propeller  pitch  commands 
equally  to  the  proper  angle.  The  pilot  can  then  vary  the  propeller  pitch  command 
by  tilting  the  steering  wheel  to  adjust  the  vehicle's  desired  forward  speed. 

T o execute  a turning  maneuver,  the  pilot  can  command  several  of  the  effectors 
in  a variety  of  ways.  Differential  commands  may  be  set  into  the  power  shaft 
governors  and/or  propeller  pitch.  The  steering  wheel  may  be  turned  to  deflect 
the  nozzle  thrust  and/or  the  rudder  pedals  may  be  used  to  deflect  the  propeller 
air  flow.  Normally,  a turn  would  be  executed  by  coordinated  commands  to  the 
nozzle  thrust  and  rudder  deflections.  A tighter  turn  could  be  csade  by  also 
commanding  the  propellers  pitch  differentially. 

Visual  feedback  is  provided  to  the  pilot  so  that  he  will  know  the  state  of  the 
effectors  and  the  vehicle.  Gas  turbine  and  power  shaft  speeds  are  examples 
of  effector  state.  Forward  speed  and  heading  are  examples  of  vehicle  state. 

See  T^le  2 Tor  the  complete  list  of  ACV-LC  displays. 

A pilot  station  configuration  of  controls  and  displays  is  shown  in  Figure  2. 
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TABLE  1.  ACV-LC  PILOT  CONTROLS 


STBD  gas  turbine  1,  2,  3 throttles 
PORT  gas  turbine  1, 2,  3 throttles 
STBD  power  shaft  governor 
PORT  power  shaft  governor 
STBD  propeller  pitch 
PORT  propeller  pitch 
Vernier  propeller  pitch 
Thrust  nozzle  angle 
Thrust  nozzle  switch 
Rudder  pedsds 
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TABLE  2.  ACV-LC  PILOT  DISPLAYS 


Forward  velocity 
Lateral  velocity 
Heading 
Pitch  angle 
Roll  angle 

Apparent  wind  speed 
Apparent  wind  angle 
Hull  height  at  CG 
STBD  propeller  pitch 
PORT  propeller  pitch 
STBD  rudder  angle 
PORT  rudder  angle 
Thrust  nozzle  angle 
STBD  power  shaft  speed 
PORT  power  shaft  speed 
STBD  gas  turbine  speeds  (3) 

PORT  gas  turbine  speeds  (3) 

STBD  gas  turbine  fuel  flow  rates  (3) 
PORT  gas  turbine  fuel  flow  rates  (3) 
STBD  gas  turbine  torques  (3) 

PORT  gas  turbine  torques  (3) 
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Figure  2.  Pilot  station  controls  and  displays 
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SECTION  III 
MATHEMATICAL  MODEL 

The  complete  mathematical  model,  as  used  in  the  simulation,  is  presented 
here  for  completeness.  Most  equations  differ  little  from  the  given  model. 
Symbology  is  different  and  the  programming  does  follow  the  symbology  as 
outlined  in  this  section.  Solution  techniques  are  described  where  special 
programming  techniques  are  needed  to  solve  equations.  No  other  programming 
techniques  are  described.  In  general,  however,  one  programming  technique, 
modularity,  is  employed  throughout  the  simui.ation  program.  A separate  module 
or  subroutine  is  developed  for  each  separate  topic. 

EQUATIONS  OF  MOTION 

This  section  describes  the  equations  of  motion  and  the  particular  force  and 
moment  components  due  to  each  effector  or  environmental  effect  as  acting 
upon  a point  mass.  These  individual  force  or  moment  components  are  then 
summed  to  yield  the  total  acting  about  each  a>-is  centered  at  the  pilot. 


FORMAT.  The  standard  form  of  the  equations  of  motion  for  the  ACV-LC  are  shown. 
Manipulations  of  these  equations  were  performed  in  order  to  eliminate 
acceleration  cross -coupling  terms.  This  produced  a set  of  nonsimultaneous 
equations  in  6 variables  wMch  could  be  solved  routinely. 

The  equations  of  motion  show  the  summation  of  forces  and  moments  about 
each  axis.  This  summation  has  been  broken  down  into  the  various  components 
and  performed  in  a separate  subroutine. 
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ur  Body  Coordinates  Centered  at 

\ Pilot  Cabin 

(X  ^ ),  = Center  of  Gravity 

CG  CG  CG  ^ 

2 

Acceleration  (Linear)  in  direction  (Ft/Sec  ) 

2 

Acceleration  (Linear)  in  direction  (Ft/Sec  ) 

2 

Acceleration  (Linear)  in  direction  (Ft/Sec  ) 

Acceleration  (Angular)  about  cixis  (Rad/Sec^  ) 

2 

Acceleration  (Angular)  about  axis  (Rad/Sec  ) 

Acceleration  (Angular)  about  axis  Z (Rad/Sec  ^) 

b 

Velocity  (Linear)  in  direction  X (Ft/Sec) 

b 

Velocity  (Linear)  in  direction  Yj^  (Ft/Sec) 

Velocity  (Linear)  in  direction  (Ft/Sec) 

Velocity  (Angular)  about  axis  X^^  (Rad/Sec) 

Velocity  (Angular)  about  axis  Yj^  (Rad/Sec) 

Velocity  (Angular)  about  axis  (Rad/Sec) 

Position  of  CG  measured  along  = -30.  0 Ft. 

Position  of  CG  measured -along  Y = -18.  0 Ft. 

b 

Position  of  CG  measured  along  Z^  = + 8.  0 Ft. 

ACV  MASS  = 10879.  5 Slugs 

6 2 

Moment  of  Inertia  about  Axis  X,  =5.  672  x 10  SLUG-FT 

^ 7 2 

Moment  of  Inertia  about  Axis  Y =1.  629  x 10  SLUG-FT 

^ 7 2 

Moment  of  Inertia  about  Axis  Z =2.  057  x 10  SLUG-FT 

b 

Summation  of  Forces  Acting  on  Vehicle  in  Direction  X (Lbs) 

b 

Summation  of  Forces  Acting  on  Vehicle  in  Direction  Y (Lbs) 

b 

Summation  of  Forces  Acting  on  Vehicle  in  Direction  Z (Lbs) 

b 

Summation  of  Moments  Acting  on  Vehicle  about  Axis  X,  (Lb-Ft) 

b 

Summation  of  Moments  Acting  on  Vehicle  about  Axis  Y (Lb-Ft) 

b 

Summation  of  Moments  Acting  on  Vehicle  about  Axis  Z (Lb-Ft) 

b 
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If  ,1F  .IF  . S M . I M . 2 M 
- X y z X ^ z 


Subscripts 

Forces  & Moments  due  to: 

PROP 

Both  Propellers 

NOZ 

Both  Thrust  Nozzles 

RUDDER 

Both  Rudders 

AERO 

Aero  Dynamics  - Windage 

AIR 

Air  Momentum 

GRAY 

Gravity 

SKIRT 

Skirt  and  Spray  Drag 

DUCT 

Both  Propeller  Ducts  - Windage 

CUSH 

Cushion  Pressure 

DAMP 

Yaw  Aerodynamic  Damping 

SEA 

Seaway  and  Vehicle  Generated  Waves 
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RUDDERS.  First,  tue  axial  velocity  of  air  through  each  propeller  duct  is 
calculated.  Second,  the  axial  air  pressure  exerted  at  the  rudders  is  calcu- 
lated. Third,  the  total  lift  and  drag  coefficients  are  calculated  based  on 
rudder  angle.  Then,  the  total  forces  and  moments  due  to  the  rudders  are 
calculated. 
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Subroutine:  RUDDERS 
RUDDERS  (Forces  & Moments) 


AXIAL  j 


^^DUCT 


AWIND 


cos  13  J ^ -i-  '^PROPi 


2 /?A 


DUCT 


^PROPj 


^PROP  j ^ ° 


AXIAL  j 


LIFT 


DRAG 


= 1/2  PV 


AXIAL  j 


• 053  tp  ij,  < 20 

K ^R  ~ 


1.  06  15_ 
I^r| 


I > 20  ° 


.02  + (.422  X 10"^)^  ^ 


^RUDDER  ■ "^DRAG  ^ RUDDER  ^ ^AXIALS  ^^AXIALP^ 


\uDDER  “ ^LIFT  ^ RUDDER  ^ ^AXIALS  ^ AXIALP^ 

K = Z Y =0 

RUDDER  R RUDDER 

M = X " = 0 

RUDDER  R RUDDER 

N =C  A (YP  +YP  \ 

RUDDER  DRAG  RUDDER  ' RS  AXIALS  RP  AXIAL  P' 

+X  Y 

R RUDDER 


V 

AWIND 

Apparent  Wind  Velocity  (Ft /Sec) 

Apparent  Wind  Angle  (DEG) 

^PROPj 

Thrust  on  Propeller  j (Lbs) 

j 

1 for  STBD,  2 for  PORT 

P 

Air  Density  = .00237  Slugs /Ft  ^ 

A 

DUCT 

. 2 

Duct  Area  = 123  Ft 

! 

\ 


II 
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^AXIAL  j 
^AXIAL  j 
R 

^LIFT 

‘"drag 

^RUDDER 


Axial  Velocity  Inside  Duct  at  Rudder  j (Ft/Sec) 

2 

Dynamic  Pressure  Inside  Duct  at  Rudder  j (Lbs/ln  ) 
Rudder  Angle  (RAD) 

Lift  Coefficient 
Drag  Coefficient 

2 

Rudder  Area  = 47.  2 Ft 

Position  of  Rudder  along  = 0 Ft. 

Position  of  Rudder  along  Y STBD  = -4  Ft 

b 

Position  of  Rudder  along  Y PORT  = -32  Ft 
• b 

Position  of  Rudder  along  = -67. 1 Ft. 
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Subroutine:  NOZZLES 


Thrust  Nozzles  (Forces  & Moments) 

■^NOZj'P- 

^NOZ  ^ ^^NOZS^^NOZ 
^NOZ  ~ ^'^NOZS^'^NOZ  ^N 

^NOZ  ' "^NOZ^N 
^NOZ  ^ ^NOZ^N 


NOZ  ' ^noz^n"'^nozs^ns^’^noz  P^NP^  ^'l 


j 1 for  STBD,  2 for  PORT 

Qnoz  j Flow  Through  Nozzle  j (Ft  /Sec) 

^NOZ  j Nozzle  j (Lbs) 


Position  of  Nozzle  (DEG) 


Position  of  Nozzle  along  X =-E3,  16  Ft. 

b 

Position  of  Nozzle  along  Z =-3.  0 Ft. 

b 

Position  of  Nozzle  along  Y =-.  25  Ft  STBD 

b 

Position  of  Nozzle  along  Y =-35.  75  Ft  PORT 

b 
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PROPELLER  DUCTS.  First,  the  wind  angle  of  attack  of  each  duct  is  calculated. 
Second,  the  side  force  drag  coefficient  is  calculated  for  each  duct.  Third, 
the  wind  velocity  acting  on  each  duct  is  calculated.  Fourth,  the  wind  velocity 
causes  pressure  on  each  duct.  Fifth,  the  total  forces  and  moments  on  both 
propeller  ducts  are  calculated. 

A Computer  plot  by  the  real-time  subroutine  of  the  side  force  coefficient  is 
shown  in  Figure  4 for  verification  purposes. 


DUCT  j 
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Subroutine;  DUCTS 

PROPELLER  DUCTS-  WINDAGE  (Forces  & Moments) 
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^'AWIND 

r 

'^PROPj 

3 

X 

‘^DUCT  3 
^'dUCT  3 

^DUCT  3 
^DUCT  3 
%UCT  3 
•^DUCT 


Apparent  Wind  Velocity  (Ft/Sec) 

Apparent  Wind  Angle  (DEG) 

Angular  Velocity  about  Axis  Zj^(RAD/SEC) 
Thrust  of  Propeller  jCLbs) 

1 for  STBD,  2 for  PORT 

Angle  of  Attack  - Relative  Wind  to  Duct  (DEG) 
Air  Velocity  Around  Duct  (Ft/Sec) 

Air  Density  = .00237  Slugs/Ft^ 

2 

Dynamic  Pressure  Acting  on  Duct  (Lbs/Ft  ) 
Duct  Side  Force  Coefficient 
Duct  Diameter  = 11.25  Ft. 

Duct  Chord  = 4.  67  Ft. 

Position  of  Duct  along  = -67, 1 Ft. 

Position  of  Duct  along  Zj^  = 0 Ft. 
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Figure  4.  Propeller  Duct  Side  Force  Coefficient 


NAVTRAEQUIPCEN  75-C-0057-1 


AERODYNAMICS  DUE  TO  i INDAGE.  From  the  apparent  wind  angle  (wind  angle  of 
attack  acting  on  the  vehicle)  the  drag  coefficients  are  calculated.  The 
forces  and  moments  are  then  calculated  from  the  apparent  wind  velocity. 

The  aerodynamic  drag  coefficients  are  plotted  by  the  real-time  subroutine  for 
verification  purposes.  The  sway  drag  coefficient  is  shown  in  Figure  5 . 

The  surge  and  yaw  drag  coefficients  are  shown  in  Figures  6 and  7 . 
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'DRAG  Y 


DRAG  X 


DRAG  N 


FRONT 


A WIND 


CG 


CG 


CG 


PLAN 


CUSH 


Apparent  V/ind  Angle  (DEG) 

Side  Force  Drag  Coefficient 

Frontal  Drag  Coefficient 

Yaw  Moment  Coefficient 

2 

Frontal  Area  = 836  Ft 

Air  Density  = ,00237  Slugs/Ft  ^ 

Apparent  Wind  Velocity  (Ft/Sec) 

Position  of  CG  measured  along  = -30.  0 Ft. 

Position  of  CG  measured  along 

Position  of  CG  measured  along  = + 8.  0 Ft- 

2 

Planform  Area  = 3200  Ft 
Length  of  Cushion  = 77  ft. 
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Figure  5.  Aerodynamic  Drag-Sway  Coefficient 
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Figure  7 . Aerodynamic  Drag-Yaw  Coefficient 
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GRAVpY.  From  the  roll  and  pitch  angles  and  weight  of  the  vehicle  the 
total  forces  and  monents  due  to  gravity  are  calculated. 


NAVTRAEQUIPCEN  75-C-0057-1 


Subroutine:  GRA\r[Ty 

Gravity  (Forces  & Moments)  (When  Integrating,  include  with  Equations  of  Motion) 

Y - -Z  0 

GRAY  " GRAY 

^GRAY  " ^'GRAy'^ 

^GRAY  = 

^gray  ""  ^gray\g‘^gray^cg 
^gray  ^ ^gray^cg*^gray^cg 

^GRAY  "^GRAY^CG'^GRAY^CG 


r 

m 

ACY  Mass  = 10879.5  Slugs 

1 

; ^ 

Acceleration  of  Gravity  = 32.  17  Ft/Sec  ^ 

^CG 

Position  of  CG  measured  along  = -30.  0 Ft 

1 

^CG 

Position  of  CG  measured  along  = -18.0  Ft 

1 

1 ^CG 

Position  of  CG  measured  along  Z,  = +8.  0 Ft 

D 

j 

I 

Euler  Angle  Roll  of  ACY  (Rad) 

j 

1 ^ 

Euler  Angle  Pitch  of  ACY  (Rad) 

1 
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Subroutine:  AIR 

Air  Momentum  (Forces  and  Moments) 

^AIR 

^^AWIND  ^FANS  '*'^FAN 

^AIR  = 

F cos  6 

AIR 

^AIR 

sin  p 

AIR 

^AIR  = 

-7  Y 

CG  AIR 

^AIR 

^CG  ’^AIR 

^AIR 

<='cg-='f>  ""air  -^cg  *air 

/> 

2 

Air  Density  = .00237  Slugs/Ft 

^AWIND 

Apparent  Wind  Velocity  (Ft/Sec) 

Pa 

Apparent  Wind  Angle  (DEG) 

^FANS 

3 

Total  Flow  Through  STBD  Fan  (Ft  /Sec) 

^FAN  P 

3 

Total  Flow  Through  PORT  Fan  (Ft  /Sec) 

^AIR 

Force  Acting  on  ACV  due  to  Air  Drag  (Lbs) 

^CG 

Position  of  CG  measured  along  = 

-30.0  Ft 

^CG 

Position  of  CG  measured  along  = 

-18.0  Ft 

^CG  • 

Position  of  CG  measured  along  = 

+8.0  Ft 

X 

F 

Position  of  Fans  measured  ulong  X 

b 

= -12.0  Ft 
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CUSIHON  PRESSURE.  Air  pressure  under  tlje  hull  causes  the  primarv  iift 
the  vehicle.  Unbalances  in  pressures  cause  moments  to  be  applied  to  ti.i 
vehicle. 
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CUSH 
i=l 
i=2 
i=3 
i=4 
X . 

Cl 

Y . 

Cl 

^CUSH  i 


Cushion  Compartment  Area  = 800  Ft 

STBD  Fwd  Cushion  Compartment 

STBD  Aft  Cushion  Compartment 

PORT  Aft  Cushion  Compartment 

PORT  Fwd  Cushion  Compartment 

Position  of  Center  of  Cushion  Compartment  along 

Position  of  Center  of  Cushion  Compartment  along 

2 

Pressure  within  Cushion  Compartment  (Lbs/Ft  ) 


cl 


X 


cZ 


X 


c3 


X 


c4 


c 1 


c2 


’c3 


‘c4 


-10  Ft. 

-50  Ft. 
-50  Ft. 
-10  Ft. 
-8  Ft. 
-8  Ft. 
-28  Ft. 
-28  Ft. 
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VEHICLE  GENERATED  WAVE  DRAG  APPROXIMATION.  A constant  drag  coefficient  of 
.5  in  both  longitudinal  and  lateral  directions  on  the  vehicle  is  used  to 
calculate  v*-  diag  forces  and  moments.  This  approximation  is  only  used  when 
the  vehicle  generated  wave  package  is  not  running.  Primarily,  this  approxi- 
mation was  made  to  facilitate  development  of  the  simulation. 


f 

m 
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Subroutine:  SKIRT 


WAVE  - CUSmON  PRESSURE  DRAG  APPROXIMATION  (Forces  & Moments) 
(For  use  without  Vehicle  Generated  Waves) 


^ SKIRT  ^ 

5 u 

I“1 

SKIRT  ^ 

5 V 

lv| 

^ SKIRT  " 

~^CG^  SKIRT 

^^SKIRT^ 

.^CG^  SKIRT 

N = 

SKIRT 

X Y -Y 

CG  SKIRT 

Velocity  (Linear)  in  direction  (Ft/Sec) 

Velocity  (Linear)  in  direction  (Ft/Sec) 

Position  of  CG  measured  along  X^  = -30.  0 Ft 
Position  of  CG  measured  along  = -18.0  Ft 

Position  of  CG  measured  along  = +8.  0 Ft 
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SKTRT  AND  SPRAY  DRAG.  The  composite  ei fects  of  the  skirt  taking  wave 
slap  and  spray  hitting  the  bow  and  sides  of  the  vehicle  are  taken  into 
consideration  by  a type  of  drag  with  constant  coefficient  of  .25. 
Over  land,  there  is  no  wave  slap  and  spray,  so  the  drag  is  zero.  Tlie 
forces  and  moments  are  then  calculated. 


i j 

[ I 
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i 


Subroutine:  SKIRT  VG 

SKIRT  AND  SPRAY  DRAG  (Forces  & Moments) 
(For  use  \vith  Vehicle  Generated  Waves) 


^P  ” *^w4  ^HULL 

^SKIRT  '^D  “ 

u\ 

^SKIRT  " "^D 

|v| 

^SKIRT  ~ ‘^SKIRT  ^P 

^skirt"  ^skirt 

^p 

N = Y 

skirt  skirt 

^CG 

C 


D 


u 

V 


X 


CG 


Y 


CG 


^HULD 


Drag  Coefficient  = . 25  , =0  over  land 

Vehicle  Surge  Velocity  along  X (Ft/Sec) 

• b 

Vehicle  Sway  Velocity  along  Y (Ft/Sec) 

b 

Height  of  Pilot  Cabin  above  water  (Ft) 

Position  of  CG  measured  along  X^^  = -30.  0 Ft. 

Position  of  CG  measured  along  Y,  = -18.  0 Ft, 

b 

Height  over  water  of  point  No.  4 (Ft) 

Height  of  Pilot  Cabin  above  Hull  Bottom  = 12  Ft. 
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StmNY  AND  VEHICLE  generated  WAVES.  Tl.e  sum  of  the  wave  heights  at  sneci- 
ik.  points  under  the  hull  of  the  waves  before  the  vehicle  arri»ea  ( ' ^ 

thd  caused  by  the  vehicle  (VOW  is  used  hfa  cSd 

o arri.e  at  the  total  torces  and  Boments  pushing  on  the  skirt  or  the  iLhlcle. 
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i 


^CUSH  1 
^CUSH  2 
^CUSH  3 
^CUSH  4 


NAVTRAEQUIPCEN  75-C-0057-1 


1-21,  C ushion  P e riphex-y  Point  s 
Total  height  of  waves  due  to  seaway  and  VGW  (Ft) 

Pressure 

Pressure 

Pressure 

Pressure 

i 


2 

within  compartment  1 (Lbs /Ft  ) 

2 

within  compartmen:  2 (Lbs /Ft  ) 

O 

within  compartment  3 (Lbs /Ft") 

2 

within  compartment  4 (Lbs /Ft  ) 


I 
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riinr  CONTROLS 

The  pilot  controls  are  the  devices  that  the  pilot  manipulates  in  order  to 
properly  maneuver  the  vehicle. 

THRUST  NOZZLES.  The  pilot  has  a steering  wheel  that  directs  the  thrust  nuzzle 
positions  190°.  He,  also,  has  a switch  to  direct  the  thrust  forward  or  aft. 
These  two  pilot  commands  are  used  to  generate  an  actual  nozzle  angle  comiiiand 
that  sweeps  360°.  Using  a constant  rate  servo  system,  the  nozzles  are  rotated 
to  obtain  the  actual  nozzle  angle.  (See  Figure  8 ) • 


NOZZLE  CONTROL 
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Subroutine:  EFPOS 
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Figure  8 . Thrust  Nozzle  Control 
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RUDDERS.  I'oot  pedals  activated  by  Che  pilot  command  a rudder  angle.  Througli 
a constant  rate  servo  system  the  rudder  angle  is  obtained.  (See  Kigure  9 .) 
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PROPELLER  PITCH  CONTROL 


P 

PCL  3 
6 

“^PCW 


J 

1 for  STBD,  2 for  PORT 

p 

PCW 

Vernier  Wheel  Command  (+  10°  to  -10°) 

^PCLj 

Propeller  Pitch  Angle  Command  from 
pilot  (+  40°) 

^PCj 

Propeller  Pitch  Angle  Command  (+35°  to 

Propeller  Pitch  Angle  (+35°  to  -40°) 

I 

■ Figure  10  . Propeller  Pitch  Control 
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TURBINES.  The  pilot  has  the  controls  to  vary  all  6 gas  turbine  speeds  and 
to  independently  control  the  power  shaft  speeds.  The  speeds  ot  the  gas 
turbines  are  determined  by  a first  order  lag  with  a time  constant  of  3 
seconds.  However,  the  speed  of  the  power  shaft  is  governor  controlled  and 
if  set  too  low  will  tend  to  decrease  the  operating  speeds  of  the  gas  tur- 
bines through  a first  order  lag  of  2 seconds. 

For  the  particular  operating  speeds  of  the  gas  turbines  and  power  shafts,  the 
optimum  speed  of  the  power  shaft  and  the  gas  turbine  optimum  shaft  horse- 
power are  calculated.  From  these  the  actual  shaft  horsepower  of  each  gas 
turbine  is  computed  and  it  is  a simple  calculation  to  arrive  at  the  shaft 
horsepower  for  each  power  shaft.  The  power  absorbed  by  each  propeller 
and  fan  is  determined  from,  its  characteristics.  The  angular  acceleration 
of  each  power  shaft  is  then  calculated  and  integrated  to  obtair  .ne  power 
shaft  speed,  port  and  starboard.  The  fan  and  propeller  speeds  are  geared 
directly  to  the  port  and  starboard  power  shafts. 
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NAVTRAEQUIPCEN  75-C-0057-1 


Time  Constant  of  Turbine  Throttle  = . 3 Sec 


Time  Constant  of  Power  Shaft  Control  = . 5 Sec 


Sampling  Period  (Sec) 

1,  n = 1,  2,  3 STARBOARD 

2,  n = 4,  5,  8 PORT 
Turbine  Speed  (RPM) 


Turbine  Speed  Command  (RPM)  (15500  <N^^<  18700) 
Power  Shaft  Speed  (RPM) 


AT 

SCj 


'SOPTn 


Power  Shaft  Speed  Command  (RPM)  (9000<Ng^j  < 16000) 
Optimum  Power  Shaft  Speed  for  each  Turbine  (RPM) 


OPTn 


Optimum  Shaft  Horsepowe-^  xor  each  Turbine  (HP) 


Ac^ual  Horsepower  Output  for  each  Turbine  (HP) 


Torque  produced  for  each  Turbine  (Ft-Dbs) 


Power  Sh-ft  Torque  STBD  (Ft-Lbs) 


Power  Shaft  Torque  PORT  (Ft-Lbs) 


FAN  j 


MAN  j 


FAN  J 


Power  Absoibed  by  Cushion  Fans  (HP) 
Manifold  Pressures  (Ltos/^Ft^) 


Fan  Air  Flows  (Ft  /Sec) 
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Fan  Efficiency  = . 85 


WIND 


Pitch  Angle  of  Propeller  (DEG) 


Apparent  Head  Wind  Velocity  along  X (Ft/Sec) 

b 


Calculated  Constant 


HP„„  . Power  Absorbed  by  Propellers  (HP) 
PROP  j 


TQ  . Torque  used  by  Fans  (Ft-Lbs) 
F AN  j 


^ Torque  used  by  Propellers  (Ft -Lbs) 


Power  Shaft  Accelerations  (RPM/Sec) 


^MACH  Machinery  Moment  of  Inertia  = 3.  533  Slug -Ft 


^FAN  j Speeds  (RPM) 

N Propeller  Speeds  (RPM) 
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COCrCPlT  CONTROLS  AND  DISPLAYS  - TEIIPORARY.  For  the  tcirporarv  cocr«it  no->i-r> 
at  NAV'TRAEOUIPCCN  tor  which  this  sirmiacion  was  developed,  -^11  controls  were 
implemented.  However,  the  pas  turbine  controls  were  lumped  together  ind  di- 
v.ded  ,-ORT,  stBD  on  the  vehicle  there  ere  3 fee  terblnL  on  LcS  PoJt 
Starboard  Side.  Instead  of  6 individual  throttles  onlv  2 were  provided  so 
tnat  t.ie  3 turbines  on  each  side  receive  the  same  command.  Table  3 shows  the 
implementecJ  controls  and  how  they  were  scaled. 


«ot  all  displays  were  implemented  as  they  would  be  on  the  vehicle.  Most 
maters  were  synchro  driven  instead  of  straight  analog  voltage  tvpe  naters. 
Only  1 ot  tne  6 gas  turbine  speeds  was  displayed.  Table  4 shows  the  imnle- 
^nted  displays  and  how  they  were  scaled.  Figure  11  shows  the  lavout  con- 
tiguration  of  the  controls  and  displays  which  were  implemented. 
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Subroutine:  DTA 
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COCKPIT  CONTROLS  AND  DISPLAYS.  Since  the  real-time  ACV-LC  simulation  was 
developed,  NAVTR<\EQUIPCEN  has  put  into  operation  a cockpit  which  more  closely 
• resembles  the  pilot  cabin  aboard  the  ACV-LC.  Table  5 shows  the  implemented 

controls  and  how  they  have  been  scaled.  Table  6 shows  the  implemented  dis- 
plays and  how  they  v;ere  scaled.  Figure  12  shows  the  current  configuration  of 
the  controls  and  displays  as  implemented. 
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Subroutine:  ATD 
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AIR  FLOK  AKD  PRF.SSURi-S 

Because  this  vehicle  is  air  cushioned,  air  is  the  dependent  n:ediu!n.  The 
flow  of  air  must  be  directed  into  the  various  compartments  with  sufficient 
volume  and  pressure  to  lift  the  vehicle  above  the  water  surface.  The 
intent  of  this  section  is  to  describe  the  solution  of  air  pressures  in  the 
compartments.  This  solution  is  not  an  easy  ta~.k  because  air  flov.-  causes 
air  pressure  and  air  pressure  causes  air  flow.  Although  the  intent  1*=  to 
describe  the  pressures,  it  is  nec.'ssarv  to  describe  in  several  sections 
component  parts  of  the  air  flow-pressure  relationship. 

FAN  AND  Ti.RUST  NOZZLE  AIR  FLOW.  Air  flow  out  through  each  thrust  nozzle 
is  calculated  directly  from  the  manifold  pressures.  Air  flow  through  the 
fans  is  calculated  from  the  manifold  pressure,  but  is  also  dependent  upon 
the  fan  speed. 
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Subroutine:  AIRFLOW 
AIR  FLOW  RATES 
^MANS 


^NOZS 


NOZ  P 


^PmANsI  pMANsj 

^]V1AN  P 
\1  |^!v:aN  pj  P MAN  p| 


QpANS  = 

2000 


= ‘^FAN  P 
^ 2000 


-640  'll^MANS"^°°l 

^MANS-^°®  - 15.  8{P,  . , ...-300) 

- 

P -3001 

MANS  1 -• 

- 640  J P 


I I 

5 -300i  P -300 

MANP  I MANP  -15.8  (P .,-300) 

IP 13001 

MANP  J 


P 

MANS 

Pressure  in  STBD  Manifold  (Lbs/Ft^) 

^MANP 

Pressure  in  PORT  Manifold  (Lbs/Ft^) 

N 

FANS 

Angular  speed  of  STBD  Fan  (RPM) 

N 

FANP 

Angular  Speed  of  PORT  Fan  (RPM) 

^NOZS 

STBD  Nozzle  Flow  Rate  (Ft^/Sec) 

°NOZP 

3 

PORT  Nozzle  Flow  Rate  (Ft  /Sec) 

^FANS 

STBD  Fan  Flow  Rate  (Ft^/Sec) 

^FANP 

PORT  Fan  Flow  Rate  (Ft^/Sec) 
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HEIGHT  OF  HULL  BOTTOM  ELATING.  The  height  of  each  hull  bottom  plating  point 
above  water  Is  calculated  from  the  height  above  mean  water  of  the  pilot, 
the  total  wave  height  (seaway  + VGiv) , the  roll  and  pitch  angles,  and  the 
X,Y,Z  positions  of  each  hull  point.  Figure  13 shows  the  definition  of  the 
hull  bottom  platirife  points  as  they  are  located  on  the  planform.  Table  7 
defines  the  coordinates  of  the  hull  bottom  plating  points  relative  to  the 
pilot  cabin. 
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Subroutine:  HULL  HT 

HEIGHT  OF  HULL  BOT"'OM  PLATING  ABOVE  WATER 

H = -Z  +0X  -«^Y  -Z 

HULL  k i HULL  k HULL  k HULL  k 

n = n +T1 

'k  SEAk  VGW  k 

h = H -n 

wk  HULL  k • k 

Number  of  Hull  Bottom  Plating  Points 
Inertial  Position  (Vertical)  of  Pilot  Cabin  (Ft) 

Euler  Angle  Roll  (RAD)  1 

Euler  Angle  Pitch  (RAD)  ^ 

% 

Position  of  Hull  Bottom  Plating  Pts  along  Xj^  (Ft)  j 

Position  of  Hull  Bottom  Plating  Pts  along  Yj^  (Ft)  | 


k 

Z. 

1 

<#> 

e 

^HULL  K 
^HULL  k 


^HULL  k Position  of  Hull  Bottom  Plating  Pts  along  (Ft) 

^SE^  k Height  of  Waves  due  to  Seaway  (Ft) 

^ VGW  k Height  of  Waves  due  to  Vehicle  Generated  Waves  (Ft) 

Total  Height  of  Waves  (Ft) 

Hhull  k Height  of  Hull  Points  above  Mean  Sea  Level  (Ft) 

h Height  of  Hull  Points  above  Water  (Ft) 
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1/  -LE  7 . "LANFOR^;  POINT  I.OCATIONS  (25) 


k 

^iULL 

^HULL 

1 

10 

-18 

12 

2 

10 

-8 

12 

o 

o 

10 

2 

12 

4 

-10 

2 

12 

5 

-30 

2 

12 

6 

-50 

2 

12 

7 

-70 

2 

12 

8 

-70 

-8 

12 

9 

-70 

-18 

12 

10 

-70 

-28 

12 

11 

-70 

-38 

12 

12 

-50 

-38 

12 

13 

-30 

-38 

12 

14 

1 

o 

-38 

12 

15 

10 

-38 

12 

16 

10 

-28 

12 

17 

-10 

-18 

12 

18 

-30 

-18 

12 

13 

-50 

-18 

12 

20 

-30 

-8 

12 

21 

-30 

-28 

12 

22 

-10 

-8 

12 

23 

-50 

-8 

12 

24 

-50 

-28 

12 

25 

-10 

-28 

12 

(Ft) 
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CUSHION  COMPARTMENT  AIR  VOLUMES.  Tho  average  height  of  each  custiion 
compartment  hell  bottom  section  is  computed  by  a parabolic  fit  in  both  X 
and  Y directions  to  ^he  individual  heights  above  water  for  that  coanart- 
ment.  The  air  volume  of  each  compartment  is  simply  the  product  of  this 
average  height  and  the  constant  planfonn  area. 





"CUSH  1 


"CUSH  2 


^CUSH  3 


"CL^H  4 


'CUSH  n 
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Subroutine:  VOLCUSH 
Average  Cushion  Compartment  Volumes 

= . 0277778  ^^"w2^\4’**w20'^"wl7^ 

+ . 4444444  h^„22 

. .0277778  (h^,5  +\-7’^*".v9"*wl8^ 

+ . 4444444  h^^.23 

= .0277773  (h,.i8-^\^9-"NvirNvl3^  llll“l<Nv2l'*wl9%'10‘'*'wl2> 

-- 4444444  h^^,24 

= .0277778  (\.l+^vl8■'^vl3"\•l5^  ^ - llllli’-<^vl6"^vl7"^v2^\'14^ 
+ .4444444  h^25 

""  "^^CUSH^CUSHn 


n = 1 4 

k = 1 25 


"CUSH 


"CUSH  n 


CUSH  n 


Number  of  cushion  compartments 

Number  of  huU  bottom  points 

Height  of  huU  bottom  points  over  water  (Ft) 

2 

Cushion  compartment  area  = 800  Ft 

Average  height  above  water  of  cushion  compartment  (Ft) 

3 

Air  volume  of  cushion  compartment  (Ft  ) 
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UAV.  1 WIPING.  The  wave  pumping  air  flows  are  simply  the  rate  of  air  volume 
chamgc  of  each  cushion  compartment.  In  other  words,  wave  pumping  is  the 

time  derivative  of  each  cushion  volume  and  is  calculated  bv  a traVe^oidai 
technique.  ' 
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Q 


PUMP  n = 


n = 1---4 

^CUSH  n 

‘^PUMP  n 
t 

A t 


Subroutine;  WV  PUMP 


^ n 

^CUSH  n^'^CUSH  n^_ 
a1 


- 1- * ''cUSH  n/^-°^CUSHn,.  ^CVSH  n,. 


At 


Number  of  cushion  compartments 

O 

Air  Volume  of  Cushion  Compartment  (Pt"") 

Rate  of  Compartment  Volume  Change  (Ft^/Sec) 
Time  (Sec) 

Sampling  period  (SEC) 
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CUSHION  COMl^ARTMENT  ESCAPE  AREAS.  The  clearances  between  the  bottur..  oI  the 
skirt  and  tlie  water  are  calculated  fron;  the  height  over  water  of  each  huJi 
point  around  the  skirt  periphery  using  a constant  skirt  height  oi  'i.b  it. 
These  cle.arances  are  parabol ically  fit  along  the  X,  Y directions  to  obtain 
an  average  water  clearance.  This  clearance  is  multiplied  by  the  average 
cushion  length  (60  ft.)  to  get  the  escape  area  for  each  cushion  compartment. 
However,  the  areas  solved  in  the  program  are  multiplied  by  a constant  for 
later  use. 
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Subroutine:  SKA  REA 
Cushion  Compartment  Escape  Area 


CLR,  = h , - h >0 
K \vk  s 


AREAj 


AREA„ 

O fi  \ 


L 

2C 

^'-'D 

6 \ 

P 

L 

rici 

6 

iCLRj  +4  CLR^  + 3 CLR3  + 8 CLR^  + 2 CLR, 


8 


2C 


D 


AREA.  -L 

4 


_£ 

2C 


CLRg-!-  4 CLRjq+  3 CLRjj  + 8 CLR^g’^  2 CLR^^ 


D 


|cLRj 


+ 4 CLR^g+  3 CLR,  8 CLR, ,+  2 CLR 


15 


‘'14 


43 


k = 1 16 


wk 


CLRj^ 

AREA^ 

AREA, 

I 

AREA, 

AREA^ 

L 

P 


Number  of  cushion  periphery  points 
Height  of  hull  bottom  points  over  water  (Ft) 

Height  of  skirt  =■4,5  Ft 

Air  gap  between  skirt  and  water  at  hull  bottom  points  (Ft) 
Cushion  compartment  1 escape  area  (Ft^) . 

' ” n 

2 

Cushion  compartment  2 escape  area  (Ft  ). 

9 

Cushion  compartment  3 escape  area  (Ft*').  ^ 

Cushion  compartment  4 escape  area  (Ft“). 

Width  of  Cushion  Compartment  = 20  Ft. 

Density  of  Air  = .00237  slugs/Ft^ 

Discharge  coefficient  = . 42 


2C 


D 


2C 


D 


P 

2C 


D 
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CL'SHFON  AND  MANIFOLD  AIR  PUElFSURES.  ..ix  simultaneous  e'luatiuns  are  written 
for  air  flow  between  the  A cushion  conipartments  and  the  2 manifolds  and 
include  the  4 wave  pumping  air  flows.  These  air  flow  equations  are  rewritten 
by  equating  Che  sum  of  air  flows  with  zero  and  by  substituting  the  6 pressure 
expressions  in  place  of  the  air  flows. 

Because  direct  solution  is  not  possible,  a Nev/ton-Raphson  technique  is  de- 
veloped. Setting  the  zero  side  of  the  equation  to  E,  a sot  of  error  etjua- 
tions  is  defined.  Taking  the  partial  derivatives  of  each  error  v^ith  resoect 
to  each  of  the  6 pressures  yields  a 6 x 6 matrix.  The  change  in  rre^.sure  is 
then  computed  by  inverting  the  6x6  derivative  matrix  and  multiplying  by 
the  error  m.atrix.  This  change  in  pressure  is  simply  added  to  the  existing 
pressure.  Iterations  on  this  solution  are  necessary  until  the  changes  in 
pressure  become  very  .'>mall. 
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Subroutine:  CUSH  P 
Compartment  Pressures 


= 

^CUSH  1 

^2 

= 

^CUSH  2 

^3 

= 

^CUSH  3 

^4 

= 

^CUSH  4 

^5 

^MANS 

P6_ 

= 

^MANP 

[dedp]  = 

fDEDli  = , 

dE. 

1 

L AJ 

AF’ 

[dedp]  W 

F 

F'+ap’ 

CRIT 

AP  ^P 

P • P 

If  CRIT  f. 

0.0004  Solution  Obtained,  Otherwise  Repeat  Solution 
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QPUMP^  + 589  ASQRT  (Pg-P^)  + 675  ASQRT  (P4-P1) 

-338  ASQRT  (P1-P2)  " AREA^  ASQRT (P^) 

+1128  Cgj^jj^^SQRT  (Pj)  (109. 0-P^) 

QPUMP2  + 589  ASQRT  (Pg-P2>  + 338  ASQRT  (Pj^  -P2) 

-675  ASQRT  (P2-P3)  ' AREA2  ASQRT{P2) 

+1128  Cgj^jj^^jASQRT  (Pg)  (109.0  - Pg) 

QPUMPg  + 589  ASQRT  (Pg-Pg)  + 675  ASQRT  (P2-P3) 

-338  ASQRT  (P3-P4)  -AREAg  ASQRT(P3) 

+1128  CgjQ^^  ASQRT  (P3)  (109.  0 - P3) 

QPUMP^  + 589  ASQRT  (P6-P4)  + 338  ASQRT  (P3-P4) 

-675  ASQRT  (P4-P1)  -AREA^  ASQRT  (P^) 

+1128  Cgj^iRT  (P4)  (109.0  - P4) 

-589  ASQRT  (P5-P4)  "589  ASQRT  (Pg-Pg)  "346  ASQRT  (P^) 

- NFANS  (640  ASQRT  (P,-300)  + 15.  8 (P,:-300)) 

200U  ■ ® 

-589  ASQRT  (Pg-Pg)  "589  ASQRT  (Pg-P4)  "346  ASQRT  (Pg) 

- NFAN  P (640  ASQRT  (P„-300)  + 15.  8 (P„-300)) 
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dEJdP,  = -589  DASQRT  (P--PJ  -675  DASQRT  (P.-P,) 

I i O 1 4 1 

-338  DASQRT  (Pj-P^)  -AREA^  DASQRT  (P,) 

-1128  ASQRT  (P^) 

+1128  DASQRT  (P^)  (109.0  - P^) 

dE^IdP^  = 338  DASQRT  (P^-Pg) 

aEj/ePg  = 0 

aE,/aP,  = 675  DASQRT  (P  , -PJ 
14  4 1 

aEj/aPj.  = 589  DASQRT  (P^  -P^) 

oE^/dPg  = 0 

d'E^ldV^  - 338  DASQRT  (P^-P2)  = ^E^/dPg 

DASQRT  (Pg-Pg)  -338  DASQRT  (Pj-P^)-675  DASQRT(P2-p3) 
-AREAg  DASQRT  (Pg)  -1128  ASQRT  (P2) 

+1128  DASQRT  (P2)  (109.0  - P2) 

8E2/aP3  = 675  DASQRT  (P2-P3) 

dE^ldF^  = 0 

dE^IdP^  = 589  DASQRT  (Pg -P2) 
aE„/3P.  - 


0 
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9E3/  aPj 

aE3/aP2 


= 675  DASQRT  (P2-P3)  = 


aEg/aPs 


-589  DASQRT  (Pg-^s)  "675  DASQRT  (Pg-Pg) 

-338  DASQRT  (P3-P4)  -AREA3DASQRT  (P3) 

-1128  A§QRT  (P3)  +1128  Cgj^j^^DASQRTCPgXlOQ.  O-P3) 


aEg/uP^ 

9^3/ 3 P5 


= 338  DASQRT  (P3-P4) 


aE3/aPg 


589  DASQRT  (Pg-P3) 


aE^/aPj 


= 675  DASQRT  (P4-P1)  =aEj/aP4 


oE^/aPg 


338  DASQRT  (P3-P4)  = 3 Eg/ 8?^ 


= -589  DASQRT  (Pg-P4)  -338  DASQRT  (P3-P4) 

-675  DASQRT  (P4-P1)  -AREA^DASQRT(P^) 

-1128  ASQRT  (P^)  +1128  Cgj^j^^DASQRT  (P4)(109. 0-P^) 


?E4/3P6 


589  DASQRT  (P6-P4) 
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aE./dP.  ^ 589  DASQRT  (P.-PJ  -dEjdP. 

’3  1 D 1 1 3 

dE^IdP^  = 589  DASQRT  (P5-P2)  =dE^ldP^ 

dE^IdP^  --  0 

dEg/aP^  = 0 

dEg/aPg  = -589  DASQRT  (P5-P1)  -S89  DASQRT  (Pg-P2)-346  DASQRT(Pg) 

- NFANS  .(640  DASQRT  (P.-300)  + 15.8) 

2000  ^ 

aEg/aPg  = 0 

aEg/aPj  = 0 

aEg/aP2  = 0 

aEg/aPg  = 589  DASQRT  (Pg-Pg)  ^dE^/dP^ 

aEg/aP^  = 589  DASQRT  (Pg-P4)  -dE^/dP^ 

aEg/aPg  = 0 

= -589  DASQRT  (Pg-Pg)  -589  DASQRT  (Pg-P4) 

-346  DASQRT  (P„)  - NFAN  P (640  DASQRT  (P.,-300)+l5.  8} 


aEg/aPg 
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^CUSH  1 

^CUSH  2 

^CUSH  3 

^CUSH  4 

^MANS 

^MAN  P 
E 

Ap 

^SKIRT 

‘^PUMP  1 

^PUiVIP  2 

^PUMP  3 

^PUMP  4 

%ANS 

^FAN  P 
AREAj 

AREAg 

AREAg 

AREA^ 


2 

Pressure  in  cushion  compartment  1 (Lbs/Ft  ) 

9 

Pressure  in  cushion  compartment  2 (Lbs/  Ft“) 

2 

Pressure  in  cushion  compartment  3 (Lbs/  Ft  ) 

9 

Pressure  in  cushion  compartment  4 (Lbs/  Ft"") 

Pressure  in  ifTBD  manifold  (Lbs/  Ft^) 

Pressure  in  PORT  manifold  (Lbs/Ft^) 

Compartment  air  flow  errors;  ideally  = 0 

2 

Change  in  6 compartment  pressures  per  iteration  (Lbs/Ft  ) 
Skirt  rigidity  coefficient  =.01 


Rate  of  cushion  compartment  1 volume  change  (Ft^/  Sec) 

O 

Rate  of  cushion  compartment  2 volume  change  (Ft  / Sec) 

O 

Rate  of  cushion  compartment  3 volume  change  (Ft  /Sec) 

3 

Rate  of  cushion  compartment  4 volume  change  (Ft  / Sec) 
Angular  spee'’  of  STBD  fan  (RPM) 

Angular  speed  of  PORT  fan  (RPM) 


2 

1 escape  area  (Ft  ) 

• Constant 

2 

2 escape  area  (Ft  ) 

• Constant 

2 

3 escape  area  (Ft  ) 

• Constant 

2 

4 escape  area  (Ft  ) 

• Constant 
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Subroutine:  WIND  PAR 


WIND 


"wind 

= " + ^WIND  ^^'^WIND^ 

^WIND 

" ■ ®WIND  ^^'^IND^ 

''  awind 

2 2 
""  " WIND  WIND 

^AWIND 

= Vv^ 

• AWIND 

Pa 

= tan'^  '"wind 

"wind 

®wind 

Wind  Velocity  (Ft /Sec) 

^WIND 

Set  of  Wind  (DEG) 

Heading  of  Vehicle  (DEG) 

u 

Velocity  (Linear)  in  direction  (Ft /Sec) 

V 

Velocity  (Linear)  in  direction  (Ft /Sec) 

"wind 

Apparent  Head  Wind  Velocity  along  Xj^  (Ft /Sec) 

^WIND 

Apparent  Side  Wind  Velocity  along  (Ft /Sec) 

^AWIND 

Apparent  Wind  Velocity  (Ft/Sec) 

Pa 

Apparent  Wind  Angle  ( + 180*^) 

1 

I 
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SEAWAY.  These  waves  are  generaced  as  a single  component  sinucoid  moving  nortiu;ard . 
These  waves  are  meant  to  simulate  swell  consisting  of  relatively  long  period 
and  wave  length.  Computation  limits  in  real-time  prohibit  the  addition  of 
other  frequency  components. 

Because  this  vehicle  is  meant  to  be  a landing  craft,  considerable  effort  in 
this  dynamic  seaway  simulation  is  devoted  to  the  creation  of  shoaling  and 
breaking  waves.  Shoaling  waves  naturally  blend  into  a single  frequency 
making  a realistic  offshore  many  component  wave  generation  superfluous. 

As  the  wave  approaches  the  beach,  which  runs  east-west,  the  wave  length 
shortens  and  the  wave  amplitude  increases  causing  the  wave  slope  to  steepen. 

Once  the  wave  slope  becomes  sufficiently  steep,  the  wave  breaks  causing 
the  amplitude  to  decrease  sharply.  Then,  as  the  wave  continues  to  encroach 
upon  the  beach,  the  amplitude  continues  to  decrease  until  it  becomes  zero  at 
the  boundary  of  the  beach. 

The  waves  are  created  completely  in  the  earth -fixed  coordinate  system.  It 
is  necessary  to  obtain  the  wave  heights  at  all  the  hull  bottom  plating  points 
for  use  in  calculating  forces  and  moments.  Consequently,  it  is  first  necessary 
to  calculate  the  inertial  components  of  all  hull  points  and  the  earth-fixed  position 
of  each  hull  point. 

Over  land  there  is  no  seaway,  so  for  the  particular  hull  point  that  is  over  land, 
its  associated  wave  height  is  set  to  zero. 

Over  the  deep  ocean  region  (offshore)  a simple  cosine -function  is  evaluated 
using  a phase  angle  based  on  position  and  the  wave  frequency. 

Over  the  shoaling  region,  the  computational  problem  begins  to  get  more  difficult . 
From  the  depth  of  water  directly  beneath  each  hull  point,  the  wave  number  and 
phase  angle  at  that  point  must  be  determined.  From  the  water  depth,  wave 
number,  and  offshore  wave  number  and  amplitude,  the  actual  wave  amplitude 
under  each  hull  point  is  calculated  and  limited.  The  wave  height  of  the  first 
order  is  a cosine  function  of  the  local  amplitude  and  phase  and  offshore  wave 
frequency.  The  second  order  wave  height  is  a cosine  function  of  twice  the 
local  phase  and  offshore  wave  frequency  multiplied  by  a polynomial  in  offshore  and 
local  wave  numbers  and  multiplied  by  the  square  of  the  local  wave  amplitude.  The 
final  wave  height  under  each  hull  bottom  plating  point  is  the  sum  of  the  first  and 
second  order  wave  heights. 

Figure  14  defines  the  ocean  boundary  regions  in  the  horizontal  plane  relative 
to  the  earth-fixed  coordinate  system.  Figure  15  defines  the  offshore  para- 
meters and  relates  the  shoaling  region  parameters  geometrically  in  the  ver- 
tical plane. 
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The  wave  number  and  phase  angle  functions  are  piecewise  fit  (f{x)=A./x^B.x+C.) 

of  9 equal  segments  and  2 half  segments  using  a least  squared  fit  and  adjusted 
for  function  continuity  at  the  segment  boundaries.  See  Figures  16  and  17 
for  computer  generated  plots  of  these  functions  in  the  shoaling  region. 

Figure  18  shows  a wave  profile  (snapshot  in  time  versus  position)  as  it  travels 
from  the  deep  ocean  onto  the  beach.  Figure  19  shows  an  expanded  section  at 
the  beach  of  Figure  18. 
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Subroutine:  CVIHQLL 
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^ULL  k ^ ■ ^HULL  k ^ I 
^ULLk^^"  -^^HULLk  I 

^“^ULL  k ^ ^HULL  k ^HULL  k 


k 

< 

^ULL  k 
^HULL  k 
^HULL  k 

4> 

0 

t' 

^iHk 

YiHk 

^iHk 


Number  of  Hull  Bottom  Plating  Points 
Vehicle  to  Inertial  Transformation 

Position  of  HuU  Bottom  Plating  Pts.  along  (Ft) 

Position  of  Hull  Bottom  Plating  Pts.  along  Yj^  (Ft) 

Position  of  HuU  Bottom  Plating  Pts.  along  (Ft) 

Euler  Angie  RoU  (RAD) 

Euler  Angle  Pitch  (RAD) 

Heading  of  Vehicle  (DEG) 

Components  of  Hull  Bottom  Points  along  (Ft) 
Components  of  Hull  Bottom  Points  along  (Ft) 
Components  of  Hull  Bottom  Points  along  Z.  (Ft) 
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Dynamic  Seaway  Generation 


Subroutine:  SEAWAY 
Initial  Calculations 
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Number  of  Hull  Bottom  Plating  Points 
OFFSHORE  Wave  Length  (Ft) 

OFFSHORE  Wave  Period  (SEC) 

1^  OFFSHORE  Wave  Frequency  (RAD/SEC) 

2 

Gravitational  Acceleration  = 32. 17  ft /sec 

) Depth  at  Which  Shoaling  Begins  (Ft) 

s 

OFFSHORE  Wave  Number  (Rad/J't) 

[ Peak  to  Peak  Wave  Height  (Ft) 


ig  Slope  of  the  Ocean  Bottom 

[g  Position  at  which  Shoaling  begins  along  X.  (Ft) 

Inertial  Position  (North)  of  Pilot  Cabin  (Ft) 

\ Inertial  Position  (East)  of  Pilot  Cabin  (Ft) 

..  Inertial  Position  (Vertical)  of  Pilot  Cabin  (Ft) 


iHk 

W 

W 

IHULL  k 

4hULL  k 

'iHULL  k 
) 

X 

X 


Components  of  hull  bottom  points  along  (Ft) 
Components  of  hull  bottom  points  along  Y.  (Ft) 
Components  of  hull  bottom  points  along  (Ft) 

Position  of  hull  bottom  plating  pts.  along  X.  (Ft) 
Position  of  Hull  bottom  plating  pts.  along  Y.  (Ft) 
Position  of  hull  bottom  plating  pts.  along  (Ft) 
Ocean  depth  ^mde^  hull  point  along  Z.  (Ft) 

Wave  number  at  hull  point  (Rad /Ft) 
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Subroutine:  SYSPIIINT 
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2 TT 


X 


X.  - C 
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w_ 


K 


K X.  - u 
XI  o 


X 


C 

A, 


X. 

1 


CJ, 


Offshore  wave  frequency  (Rad/Sec) 

Wave  number  at  pilot  station  (Rad/Ft) 

Wave  velocity  along  X^  (Ft /Sec) 

Wave  length  at  pilot  station  (Ft) 

Vehicle  velocity  along  Xj  (Ft /Sec) 

Wave  encounter  frequency  (Rad/Sec) 
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Ocean-Land  Definition 


North 


Figure  15.  Ocean  Bottom  Profile 


Figure  17.  Wave  Phase 
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VEHICLE  GENERATED  WAVES 

This  section  of  the  mathematical  model  is  the  most  complex  part  of  the  whole 
simulation.  The  implementation  and  the  concept  are  reasonably  straight 
forward,  but  the  mathematics  of  the  kernel  function  and  the  development  of 
the  model  are  sophisticated  and  complex.  Although  the  implementation  is 
straightforward,  the  number  of  calculations  is  voluminous.  This  fact  presents 
great  difficulty  in  doing  these  ca’-  ilations  in  real-time. 

The  vehicle  rides  on  a cushion  -ir  which  is  obviously  pressurized  to 
counteract  the  force  of  gravit;) . The  higher-than-atmospheric  pressure  causes 
a depression  in  the  water  surface.  If  the  craft  is  not  moving,  i.  e. , just 
hovering,  this  surface  depression  is  static  in  nature.  This  situation  can 
be  explained  as  a sequence  of  pressure  impulses  applied  on  the  free  surfaces 
for  a prolonged  period  of  time.  If  the  duration  over  which  these  impulses 
are  applied  is  short,  a small  number  of  radiating  waves  would  be  generated. 

In  the  limit,  however,  the  shape  under  the  cushion  would  be  static  in  shape. 

The  present  implementation  only  allows  for  a finite  length  of  time  for 
application  of  impulses  (To  = 8 sec).  Hence  the  "static"  froe-surface 
elevation  will  still  be  wavy  outside  of  the  air  cushion.  It  is  this  wave 
pattern  which  constitutes  the  vehicle  generated  waves.  No  concern  is 
given  to  the  influence  of  these  waves  beyond  the  confines  of  the  hull. 

After  a wave  pattern  is  created  and  the  vehicle  moves,  its  forward  motion 
will  create  another  wave  pattern.  As  the  vehicle  moves  over  the 
established  wave  pattern,  distortion  of  the  existing  wave  pattern  occurs 
under  the  vehicle  because  of  the  newly  generated  pattern.  This  changing 
of  wave  shape  from  one  pattern  to  the  next  mathematically  can  be 
expressed  as  a convolution.  Therefore,  to  create  a complete  profile  of  the 
vehicle  generated  wave  pattern,  it  is  necessary  to  convolve  all  wave  shapes 
generated  in  the  past  with  the  present. 

The  pressure  distribution  under  the  craft  used  for  rigid  body  dynamics 
is  dependent  upon  the  air  volume  and  leakage  gaps  across  the  complete 
planform.  In  order  to  know  the  complete  distribution  of  the  air  volume, 
it  is  necessary  to  determine  the  height  between  planform  and  water  at  all 
points  under  the  hull-  This  is  synonymous  with  knowing  the  height 
of  the  wave  at  each  point  with  respect  to  mean  water  level. 

The  model  of  the  vehicle  generated  waves  is  only  concerned  with  determining 
these  wave  heights  by  assuming  a time-independent  pressure  distribution. 

The  change  in  pressure  and  air  volume  distributions  are  considered  elsewhere 
using  the  vehicle  generated  wave  heights. 


WAVE  HEIGHTS.  The  vehicle  generated  wave  heights  are  determined  under  a selective 
number  of  points  of  the  hull  planform.  This  determination  is  made  by 
performing  a convolution  integral  over  the  kernel  function  at  specific 
time  intervals  in  the  past. 
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The  kernel  function  is  a three  diirjcnsional  wave  height  function  of  position 
(X,  Y)  and  time.  This  function  is  pretabulated  versus  X and  V for  specific 
sampling  periods  past  (t-nAt).  In  other  words,  a wave  pattern  is  established 
at  each  time  period  in  the  past  by  table  lookup  as  a function  of  position 
(how  far  the  vehicle  has  travelled  since  being  at  the  previous  point).  These 
kernels  are  then  convolved  over  all  time  yielding  the  wave  height. 

Over  land  there  is  no  water  to  be  displaced,  the  associated  VGW  height 
is  set  to  zero.  If  the  complete  vehicle  is  over  land,  the  skirt  and  spray 
drag  coefficient  is  set  to  zero. 

The  history  of  the  trajectory  is  stored  and  updated  at  each  time  step.  This 
trajectory  history  consists  of  the  change  in  inertial  position  of  the  vehicle 
from  the  present  and  the  heading  of  the  vehicle  at  each  time  step  in  the  past. 
The  change  in  inertial  position  of  each  hull  point  is  calculated  from  the 
vehicle  change  in  position  and  the  inertial  hull  point  component.  This 
change  in  position  is  transformed  from  earth-fixed  coordinates  to  body-fixed 
coordinates  for  each  time  period.  This  change  in  position  is  then  referenced 
to  the  vehicle  center  of  gravity.  Because  the  /GW  is  symmetrical  about  the 
body  center,  only  the  magnitude  of  this  change  in  X,  Y position  is  needed  to 
arrive  at  the  kernel  value.  All  kernel  values  are  integrated  by  Simpson's 
Rule  and  dimensionalized  in  time  and  position  to  yield  the  final  VGW  height 
under  each  planform  point. 

Complete  plots  were  generated  to  show  the  physical  shape  of  the  Vehicle 
Generated  Wave  (VGW).  These  plots  show  the  wave  for  vehicle  straight 
line  motion  and  constant  forward  speed  for  all  time.  Three  sets  of  plots 
are  included  for  vehicle  forward  speed  of  0,  30,  and  84  ft/ sec  corresponding 
to  0,  18,  and  50  knots.  Each  set  consists  of  three  plots.  The  first  plot 
shows  the  VGW  height  (ri)  along  the  longitudinal  direction  of  craft  motion. 

The  second  plot  shows  the  VGW  height  (h)'  along  the  lateral  direction  of  the 
vehicle.  The  third  plot  shows  the  VGW  height  .(ri)  in  time  at  a fixed  point 
on  the  ocean  surface  as  the  vehicle  passes  over  it. 

The  first  i"  ->  wave  profiles  (see  Figures  20,  21,  23,  24,  26,  and 
27)  are  > dimensional  cuts  (slices)  along  the  vehicle  body  axes  (x,  y) 
centered  . at  the  CG  cb  snapshotted  at  a particular  instant  in  time; 
although,  u wave  profile  as  seen  from  the  vehicle  does  not  change  when 
its  speed  is  constant.  The  first  wave  profile  (see  Figures  20,  23, 

26)  versus  x position  shows  the  CG  at  x=0  with  the  bow  at  x=40  and  the 
stern  at  x=-40.  The  second  wave  profile  (see  Figures  21,  24,  27  ) 

versus  y position  shows  the  CG  at  Y =0  with  the  starboard  side  at  y=20  and 
the  port  side  at  y=20. 

The  third  wave  profile  (see  Figures  22,  25,  28  ) shows  the  two 
dimensional  wave  height  effect  of  a fixed  point  on  the  ocean  surface  as 
the  centerline  of  the  ACV  approaches  and  passes  over  the  point.  At 
t=0  the  CG  is  directly  over  the  point. 
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Subroutine:  VG  WAVE 

Over  Land 
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k Number  of  Hull  Bottom  Plating  Points 

n =1, . . . , N Number  of  History  states,  N is  an  even  positive  integer 


VGW  k 


'D  SKIRT 


^iHk 

^iHk 


Ay 
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VGW  Wave  Heights  (Pt) 

Skirt  Drag  Coefficient 

Current  inertial  component  of  Hull  Point  along  (Ft) 

Current  inertial  component  of  Hull  Point  along  Y-  (Ft) 

Current  inertial  component  of  Hull  Point  along  Z.  (Ft) 

Vehicle  Change  in  inertial  position  of  each  history  state 

along  X.  (Ft) 

Vehicle  Change  in  inertial  position  of  each  history  state 

along  Yj  (Ft) 

Vehicle  Change  in  inertial  position  of  each  history  state 

along  Zj  (Pt) 

Change  in  inertial  position  of  hull  pt.  along  X-  (Ft) 

Change  in  inertial  position  of  huE  pt.  along  Y^  (Ft) 

Change  in  inertial  position  of  hull  pt.  along  Z.  (Ft) 

Inertial  to  vehicle  transformation  for  each  history  state 

Change  in  position  of  hull  pt.  along  (Ft) 

Change  in  position  of  hull  point  along  Y^  (Ft) 

Change  in  position  of  huU  point  along  (Ft) 

Magnitude  of  change  of  position  of  huH  point  referenced  to 
CG  along  (Ft) 

Magnitude  of  change  of  position  of  hull  point  referenced 
to  CG  along  Y^,  (Ft) 

Magnitude  of  change  of  position  of  hull  point  referenced  to 
CG  along  (Ft) 
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Position  of  CG  measured  along  Xj^  = 30.  0 Ft. 

Position  of  CG  measured  along  Yj^  =-18. 0 Ft. 

Position  of  CG  measured  along  = +8.0  Ft. 

Kernel  •''alue  for  history  states 
State  or  kernel  at  t-NAt 

e 

Nominal  Pressure  in  Cushion  Compartment  =109.375  Lbs/Ft 
Density  of  Water  = 1.98  slugs/Ft^ 

2 

Gravitational  acceleration  = 32.17  Ft /Sec 
Vehicle  half  length  = 40  Ft. 
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Figure  23,  VGW  Profile  - u = 30  ft/ sec,  Longitudinal 


Profile  - u = 30  ft/ sec.  Time 


Profile  - u = 84  ft /sec.  Longitudinal 


igure  27.  VGW  Profile  - u = 84  ft/sec,  L/ateral 


Figure  28.  VGW  Profile  - u = 84  ft/sec.  Time 
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Kr-RNEL  VALVE.  The  vehicle  generated  wave  is  created  In  the  present  and  has  existed 
for  all  time.  The  kernel  function  represents  the  nondimensionallxed  wave  height  at 
a particular  instant  in  time  due  to  an  impulse  in  pressure.  Therefore, 
the  table  consists  of  kernel  values  t-nAt) 

versus  X,  Y position  for  a discrete  time  (t-nAt).  With  n known,  it  is 
held  constant  for  the  kernel  value  K ^vn  From  the 

position  of  where  the  wave  height  is  to  be  oon^med,  the  kernel  value  can  be 
obtained. 

In  the  table,  a particular  kernel  exists  for  discrete  X,  Y positions.  The  kernel 

function  is  assumed  to  be  linear  between  tabulated  points  of  known  X,  Y position. 

The  particular  R , R position  may  lie  between  other  X,  Y positions  in  the 

table.  Linear  interpolation  is  performed  on  the  4 kernel  values  at  the  4 X,  Y 

positions  immediately  surrounding  the  desired  R ,,  R . 

X y 

Computer  plots  for  the  kernel  function  were  generated  for  each  unit  of 
time  from  the  VGW  table.  A family  of  6 curves  was  plotted  for  each  sub- 
table  of  kernel  (Kxy)  versus  longitudinal  position  (Rx)  and  lateral  position 
(Ry).  These  plots  were  created  by  choosing  alternate  Ry's  (Rx's)  and  hold- 
ing the  value  of  Ry  (Rx)  constant  and  then  sweeping  through  each  Bx  (Ry) 
in  the  subtable.  Nbndimensionalized  kernels  were  plotted  versus  position 
for  times  2,  4,  6,  8 seconds  corresponding  to  subtables  numbered  22,  42, 
62,  82  (see  Figures  29,  30,  31,  32). 
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Subroutine:  IKINTRP 


Retrieve  K . for  the  function  K (R  ^ ) 

xy^  xy  xn,  yn  n 

For  each  r values  of  the  kernel  function  are  pretabulated  for  each 


X,Y  position. 


Constant 

n 


K (R  .R  ) 
xy  xn'  yn 


Linear  interpolation  on  Kj  , Kg , Kg,  K ^to  obtain  K^ 
Weighting  factors: 

iTrfn  _ X^”R..  irrnr*  ^ f \TTrrt 


K-K  1 

p p-i 
n _ Y -R 

q q-1 


WT^g  = l-WT^,^ 


WT  „ = 1-WT  , 
y2  yi 


^xy  = 'VT^2  ^Tyl  ^2  ^ ^^'^x2^^^y2^4 
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I’OoITlOS  HISTORY  OF  ACV-LC.  As  has  already  been  discussed,  it  is  necessary  to  kr.'-.f 
all  past  positions  of  the  vehicle  in  order  to  determine  the  vehicle  generated  wave 
( VGW)  heights.  However,  it  is  not  actually  necessary  to  know  the  exact 
position,  but  the  change  in  position  for  each  time  step  up  to  the  present. 

This  position  has  2 horizontal  components,  north  and  east.  The  variables 
are  stored  in  a table: 

6X.  change  in  north  position  for*  each  time  step 

^ Yj  change  in  east  position  for  each  time  step 

The  VGW  model  needs  the  boefy  a>ds  change  in  position,  necessitating  the 
saving  of  angles.  With  roll  and  pitch  angles  small,  assume  a change  only 
in  the  horizontal  plane  and  only  save  vehicle  heading.  Transformations 
need  the  direction  cosines  of  the  angles,  however.  Therefore,  the  TOriables 
are: 


Cb  Cos  ^ at  each  time  step 
Si  Sin  i at  each  time  step 

The  table  of  past  positions  becomes  successive  entries  of: 


1 

Ci 

si 

These  positions  are  saved  for  t-nA  t,  n=  0, 1, 2, 3, , Nowhere  t is  the  presort 
time  step. 

The  position  history  over  all  time  would  consist  of  a voluminous  amount  of 
data  that  would  exceed  the  core  storage  of  the  compirter;  as  well  as,  the 
convolution  over  all  time  would  be  an  impossible  number  of  computrons. 

The  number  of  past  position  states  has  been  limited  to  =160. 

The  position  entries  to  tiie  table  are  performed  top  down  witii  the  present 
position  at  tiie  top.  The  true  heading  is  inserted  at  each  time  step,  but 
the  linear  change  in  {K>sition  at  the  present  time  step  has  to  be  added  to 
each  previous  change  in  position  to  yield  the  total  change  in  position  from 
the  present  to  that  particular  time.  For  example,  the  table  looks  like  the 
following  for  only  a change  in  north  position: 


pii  iw.'  F > - 


- . _ 
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Ax 

3At 

« 

AXg 

• 

• 

-X,At 

• 

V-1 

AX. 

N:-Nt  At 

T 

At  the  next  time  step,  t +At,  AX.  =A  X 

VAt 

is  summed  with  all  entries  and  the  previous  table  becomes: 


AXj+AX^ 

AX.  +AX^ 

AX^AX” 
o n 


+AX„ 


Where  AX^  =AX.  has  become  AX^  - +Ax  . The  previous  AX^^ 

H-NrAt  7- 

is  dropped  from  the  bottom  of  the  table  and  lost.  This  table  has  N^+1  entries, 

N.J.  past  positions  plus  the  present  position. 

However,  in  the  interest  of  limiting  computations,  the  convolution  integral 
need  not  be  performed  over  all  position  history,  but  for  selected  time  intervals. 
These  positions  are  stored  in  an  array  (STATE)  for  t-nA  t:  n=0, 1,  2, 3,  4,  6, 8. . . , 
The  array  length  consists  of  /2  +2  entries.  The  remaining  position 
histories  are  stored  in  a separate  entry  (ODDSTATE)  for  t-nA  t:n=5, 7,  9. . . , -1. 

The  array  length  consists  of  N-j-  /2 -2  entries.  This  table  storage  and  update 
is  performed  by  Subroutine:  HISTORY. 

Initialization 

The  position  table  histories  must  be  imtialized  in  order  to  begin  execution. 

All  entries  are  set  equal  to  the  initialized  values  of  change  in  linear 
position  and  heading.  These  4 quantities  are  stored  successively  throughout 
both  tables; 
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TABLE  GENERiXTION.  The  Vehicle  Generated  Wave  (VGW)  table  consists  of  82  subtabies  o! 
distinct  kernel  values.  Each  subtable  corresponds  to  a pressure  wave  versut»  x,  y 
position  for  a discrete  time  in  the  past.  The  last  subtable  (number  82)  corresponds 
to  8 seconds  in  the  past.  The  x,  y positions  are  referenced  to  the  CG  in  body- 
coordinates  and  are  magnitude  only. 

The  kernel  table  is  generated  by  a FORTRAN  program  TABKS,  The  spline 
curve  data  points  are  read  and  a set  of  spline  curve  coefficients  are  computed. 

At  which  point,  the  program  reads  the  data  to  begin  table  generation.  This 
data  consists  of  nondimensionalized  time  and  a list  of  nondimensionalized 
X and  y positions.  For  each  x position,  the  y position  list  is  swept  to  compute 
a kernel  value  for  each  combination  of  x,  y position.  The  number  of  kernel 
values  computed  equals  the  number  of  x values  times  the  number  of  y values. 
Consequently,  the  computational  ability  has  been  limited  to  30  x's  by  30  y's. 

After  the  subtable  is  computed,  it  is  formatted  and  written  on  magnetic  tape. 

Data  for  time  and  position  is  then  read  for  the  next  subtable  computation 
and  the  calculations  are  repeated  for  a new  time. 

These  subtables  were  checked  for  curve  smoothness  by  an  eyeball  method. 

Two  computer  plots  were  generated  per  subtable,  one  versus  R , and  the  other 
versus  Ry,  The  fifth  R^.  (Ry)  in  the  subtable  was  randomly  chosen  as 
represent  .tive  and  all  Ry's  (Rj^'s)  were  swept.  These  plots  were  visually 
selected  on  piecewise  linear  curve  fit.  Discontinuities  were  reduced  by 
selectively  adding  an  R^  (Ry)  to  the  input  position  data.  The  subtable  was 
recreated  by  only  calculating  the  kernel  values  which  did  not  already  appear 
in  the  subtable.  The  final  check  of  curve  smoothness  was  accomplished  by 
plotting  a family  of  6 curves  per  subtable  (see  Section  3.5.2)  and  visually 
determining  the  necessity  of  adding  other  points. 
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TABKS:  Compute  kernel  values  from  the  equation 


{R.  R.  t)  = / “ 

XV  Jo 


k dk 


Y sin  Y t 


2 i+j 

3:  (-1)  Q(k,R  4^  ) 

j = 1 


t is  non-dimensionalized  time 

where  R and  R are  the  non-dimension  coordinates  of  the  point  vhere  the 
X y 

kernel  is  being  computed  relative  to  the  vehicle  eg.  in  the  vehicle  coordinate 
system  at  non-dimensional  time  t 

k is  the  radial  component  of  the  conversion  of  the  two  dimensional  inverse 
fourier  transform  to  polar  coordinates 

a is  a diffusion  parameter  (set  to  0.5) 

Y = y g k tanli  (kh)  ~ 

g = acceleration  of  gravity 

h = non-dimensional  depth  of  water  (assumed  1.0) 

and  v are  polar  coordinates  of  (R^,  R^)  referred  to  the  four  corners 
of  the  craft 


2 i+3 

S (-1)  Q (k,  R..,  q'.y  integral  of  the  angular 

j = 1 13  i] 

component  of  the  inverse  Fourier 
transform 

FINDKO 

Determine  a limit  of  integration  (k^)  to  replace  the  infinite  limit,  such  that 

' ^ 2 ( I /k^  a -ri-  ‘'o 

where  t is  the  maximum  allowable  error  (set  to  .01) 


2 

y 

i = 1 
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INTDK2 

Perform  the  integration  over  k using  a Gauss -quadrature  method. 


WNING 


Compute  the  integrand 
k Y sin  ( Y t) 


2 

T 


i = 1 .1  = 1 


(-1)  Q (k,  R..,  ^ ..) 


Q24 

Perform  the  double  sum  and  compute  the  function 


Q (k,  V ij) 


= - v 


n = 1,  3,  5, 


Cn  (k’) 


where  k'  = 


and  «■'  = 


J2jj  = 2n  " order  Bessel  function 

(k')  are  coefficients  which  are  fitted  with  a spline  curve 


PREP:  Read  in  coefficients  of  spline  curves. 

GCNK4:  Compute  C^  (k’)  from  spline  coefficients 

BESSEL:  Compute  the  Bessel  function  (kR)  in  double  precision 

SESSEL:  Compute  the  Bessel  function  (kR)  in  single  precision 
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VGV:  CALCULATIONS.  It  has  already  been  stated  that  there  are  a volumir.ous 
number  of  calculations  to  determine  all  VGw'  heiehts  under  the  hull  plunforir.. 
So  many  tnat  the  calculations  cannot  be  performed  in  real-tip:e.  The  Vf-l' 
height  calculations  are  performed  in  the  background  and  are  processed  throuph 
a first  order  lag  in  real-time  to  prevent  a step  in  VGU  height. 


NAVTRAEQUIPCEN  75-C-0057-1 


Subroutine:  FILTER 


n 


VGW^  VCW  ® 1 ™WnEW 


At 


n 


VGW  "VGW  VGW^g^ 

k 


k 

1- 


-41_ 

^ VGW 


VGW 


k 

At 

T 


n 


VGW 


VGW. 


VGW 


Number  of  hull  bottom  plating  points 

Sampling  period  = . 05  Sec 

First  order  lag  time  constant  = . 05  Sec 

Vehicle  generated  wave  height  (Ft) 
Newly  calculated  VGW  height  (Ft) 


NEW, 
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STATE  INTEGRATIONS 

An  orthogonal  Euler  transformation  matrix  is  defined  (.-;ee  table  8 Yco  relate  body 
coordinates  of  the  vehicle  to  the  earth -fixed  coordinate  system.  With 
the  transformations  defined,  the  Basic  Motion  Equations  (see  table  9 )can  also  be 
defined  to  show  how  all  positions,  rates,  and  accelerations  are  related. 

Two  integrations  are  necessary  to  take  the  output  of  the  Equations  of  Motion 
(acceleration)  to  arrive  at  position.  Tlit,?.  are  accomplished  by  a second  order 
predictor -corrector  numerical  integratior  This  scheme  executes  very 
quickly  and  only  needs  a single  execution  or  the  Equations  of  Motion.  The 
stability  and  accuracy  of  the  predictor -corrector  are  very  good  for  such  short 
computation  times. 

The  second  order  predictor -corrector  method  used  is  shown  in  one  variable 
to  facilitate  ease  in  reading.  The  scheme  in  all  variables  is  also  shown 
for  completeness. 

The  integration  scheme  simply  predicts  position  and  velocity,  corrects 
position  and  velocity,  corrects  position,  solves  the  Equations  of  Motion,  and 
corrects  velocity.  The  velocity  integrations  are  performed  in  the  earth- 
fixed  coordinate  system.  Summing  initial  conditions  with  the  integrated 
values  is  also  necessary. 

Because  the  scheme  is  of  second  order,  two  past  values  are  needed  for 
each  integration.  A simple  Euler  method  is  used  to  get  the  scheme  started, 
along  with  an  initialization  pass. 
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TABLE  8 . EULER  ROTATION  DEFINITION 
Subroutine:  C THREE 


Coordinate  Transformation-  Inertial  to  Vehicl*^ 


r 1 

“ " 

— 

0 

4= 

X y z 


^11 

'^12 

A 

"13 

A^ 

A^ 

A 

21 

22 

23 

A,  , 

A,, 

A 

31 

32 

33 

~cec4 

; cos^> 

« -S  0 -1 

-cqisjp  +s0sez'4> 

; C9ic4  +s<pses4) 

> S^C  0 

S^Sd.  +C<^SQC4 

;-S0c#  +c<i>ses4 

j C5SC  0 

i.-_e 

-S4>  \ cii>  r <4 

^ j — 


0s^+ec^:  i 


For  Small 
6,  0 {<10°) 


0 Euler  Angle  Roll  of  ACV  (RAD) 

0 Euler  Angle  Pitch  of  ACV  (RAD) 

^ Euler  Angle  Heading  of  ACV  (DEG) 
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TABLE  9 . BASIC  MOTION  EQUATIONS 

Subroutine 


i 

f 

• 

U 

dt 

P' 

o J 

• 

f 

V 

dfc 

INTGRT 

= 

oj 

L 

ft 

• 

dt 

oJ 

w 

r* 
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SECOND  ORDER  PREDICTOR  - CORRECTOR  INTEGRATIONS 


Scheme  in  One  Variable 


For  t=0.  Where  X = U = 0 
t t 

X = X + X 
o t 


U = U + U 
o t 


X = C U 
t i 


X = X 
t -2At  t 


X = X 
t-  At  t 


APC  = 0 


Gravity  Forces  and  Moments 


Sum  Forces  and  Moments 


Equations  of  Motion 


U . = U 

t - 2 At  t 


U = 0 

t-  At 


APC  = 0 
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For 

X 

t 

U 


t = At 

A t 
2 

A t 


t - A t 


U 


t-  A t 


X = X + X 
o t 


U 


U + X^ 
o t 


X = C U 
t i 


X =4 
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[ ^ (X  + X ) + X ] 
5 2 ^ t- t-ZAt 


X = X +X 
o t 


Gravity  Forces  and  Moments 
Sum  Forces  and  Moments 
Equations  of  Motion 


= i t ( ij  + U ) + U 1 

^ 5 2 t t -At  . t - 2At 


u = u + u 

o t 

• V 

X = c u 

t i 


u = u 

t -At  t 


• • 

U ■ = u 

t - At  t 


X = X 

t-At  t 

• • 

X =x 

t-At  t 
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I I I 


X = X + X 
o t 

u = u + u 

o t 

• V 

X = C u 
t i 

c =x  +_M(X+X  ) 

t t-At  2 t t-At 

APC=  4/5  (C  - P ) 
t t 

X = P +aPC 
t t 

X = X + X 
o t 

X , = X 

t-2^t  t-At 

X = X 

t-At  t 

Gravity  Forces  and  Moments 
Sum  Forces  and  Moments 
Equations  of  Motion 
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• « 


APC  = 4/5  (C^  - P^) 


U = P +APC 
t t 


U = U 4 U 
o t 


• V 

X = C U 
t i 


^t-ZAt  ~ ^t-At 


U = U 
t-At  t 


U = U 
t-At  t 


X = X 
t-At  t 
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I 


X 

0 

U 

0 

Ufc 

X 

u 


X 


t-2At 


X 

t-At 

APC 


u, 


t-2At 


APC 
At 
U 


t 


u 


t-At 


t-  At 


* 

U 


t-  At 


Initial  Inertial  Position 
Initial  Body  Axis  Velocity 
Integrated  Inertial  Velocity 
Integrated  Body  Axis  Acceleration 
Inertial  Position 
Body  Axis  Velocity 

Vehicle  to  Inertial  Transformation  Matrix 
Inertial  Velocity 

Inertial  Position  at  2 Time  Periods  Past 

Inertial  Velocity  at  1 Time  Period  Past 

Position  Corrector  Value 

Body  Axis  Velocity  at  2 Time  Periods  Past 

Velocity  Corrector  Value 

Sampling  Period 

Body  Axis  Acceleration 

Body  Axis  Velocity  at  1 Time  Period  Past 

Inertial  Position  at  1 Time  Period  Past 

Predicted  Inertial  Position 
Predicted  Body  Axis  Velocity 
Corrected  Inertial  Position 
Corrected  Body  Axis  Velocity 

Body  axis  acceleration  at  1 time  period  past 


) 
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Chant^e  in  Inertial  Position  for  Vehicle  Generated  Wave 


For  t = 0 


AX.  = 0 

AY.  = 0 
1 

For  t =At 

A X.  -AX.  = X. 
^ H-At  ^ 


AY  =AY  =Y 
^ h-At 

For  t = nA  n > 1 

AX,=  i [2At(2Xt.^t  +X,) 


AX.  =AX. 
H-At  ^ 

AY,  =AY 
H-4t  ^ 


l> 

Change  in  position  along  X,  (Ft) 

Change  in  position  - ong  Y^  (Ft) 

AX, 

H-At 

A X^  at  one  sampling  period  past  (Ft) 

^Yi 

H-At 

A Y,  at  one  sampling  period  past  (Ft) 
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Scheme  in  Total 

For  t = 0.  Where 


Subroutine:  INTGRT 


Subroutine 

SYSACVIC 


SUMICS 


SUIVIICS 


SUMICS 


+ 


SUMICS 
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Gravity  Forces  and  Moments 
Sum  Forces  and  Moments 
Equations  of  Motion 


gravity 

forces 


motion 
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153 


mil  rr 


For  t=At 
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Subroutine 


SUMICS 


SUMICS 


SUlvUCS 


154 
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Subroutine 


SUMICS 


SUMICS 


CTHREE 


CTHREE 
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Gravity  Forcet)  and  Moments 
Sum  Forces  and  Moments 
Equations  of  Motion 

u'')  r Tu  + u.  .1  fu  ^ .SI 


Subroutine 

SUMICS 


SUMICS 


SUMICS 

GRAVITY 

FORCES 

MOTION 
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Subroutinp 


SUMICS 


SUMICS 


CTHREE 


CTHREE 
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A PC 


u 


A PC 

U3 

[a  PC 
A PC 

Uo' 


Subroutine 


SUMICS 


-180°<</<i80° 


SUMICS 
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Subroutine 

SUMICS 


SUMICS 


SUMICS 


CTHREE 


CTPIREE 
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Subroutine 


180°<  <180‘^ 


SUMICS 


SUMICS 
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Gravity  Forces  and  Moments 
Sum  Forces  and  Moments 
Equations  of  Motion 


Subroutine 

GRAVITY 

FORCES 

motion 
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Subroufcinci 


SITMICS 


SUMICS 


CTHREE 


CTHREE 
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TRANSCENDENTAL  FUNCTIONS 


The-  transcendental  functions  are  the  standard  arithmetic  functions  that 
are  implemented  as  library  subroutines  for  general  purpose  usage. 

ARCTANGENT.  This  function  was  taken  directly  from  the  Sigma  7 FORTRAN 
library  (6ATAN) . It  is  a two  argument  function.  The  one  argument  entry 
was  removed.  See  Figure  33  for  a coimjutor  generated  plot. 

COSINE.  This  function  was  taken  directl:^  from  the  Sigma  7 FORTRAN  library 
{6C0S).  The  sine  entry  was  removed.  See  Figure  34  for  a computer  generated 
plot. 

S!NE  AND  COSINE.  These  functions  were  specially  developed  as  curve  fits  for 
fast  execution.  The  cosine  is  simply  90®  phase  shifted  from  the  sine  with 
angle  inputs  limited  to  +360°.  This  function  is  a piecewise  linear 
(:(x)  - A;X  + Rj)  fit  of  equal  10°  segm.ents  usine  a nin-max  (minimir.e  the 
maximum  error)  fit  and  adjusted  for  function  continuity  at  the  segment 
boundaries.  The  maxim.um  absolute  error  equals  .0019  at  +80°  and  +110°. 

See  Figures  35  through  40  for  the  functions,  absolute  error,  and  per- 
centage error. 

SQUARE  ROOT.  This  function  was  taken  directly  from  the  Sigma  7 FORTRAJJ 
library  (6SQRT) . 

Other  function  approximations  were  developed  specially  for  fast  execution. 
This  function  is  a piecewise  linear  (f(x)  = AjX  + Bj)  fit  with  segment  length 
chosen  allowing  a maximum  error  of  .5%  using  a min-nax  fit.  The  maximum 
absolute  error  lies  at  the  lower  boundary.  This  function  (7SGNSQR)  is 
a signed  square  root: 


* V ifr 


It  is  implemented  as: 


f(x)  = 


A^X  + 


AiX  - Bi 


See  Figure  41  for  computer  plot  of  function. 
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Another  function  (7ASQRT)  is  the  same  except  that: 

f(x)  = X IX!  < 1 

See  Figure  42  for  computer  plot  of  function. 

DERIVATIVE  OF  SQUARE  ROOT. 

This  function  is  a piecewise  hyperbolic  (f(x)  = Aj/X+Bj)  fit  with  segment 
length  chosen  allowing  a maximum  error  of  2%  using  a min-max  fit  The 
maximum  absolute  error  lies  at  the  lower  boundary.  This  function 
(7  DASQRT)  is  a signed  square  root  derivative: 

{1 IXI  > 1 

2 yixT 

1 IXI  S 1 

See  Figure  43  for  computer  plot  of  function. 
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Figure  38.  Cosine  Absolute  Error 





Figure  AO.  Cosine  Percentage  Error 





Figui’e  41.  Signed  Square  Root  Approximation 


Figure  42  , Signed  Square  Root  Approximation  with  Linear  Region 


Figure  A3  . Signed  Square  Root  Derivative  Approximation 
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CRT  DOCKING  AID-VISUAL 

In  order  to  give  the  pilot  visual  affects  that  would  simulate  those  on  the  craft, 
a line  drawing  was  created  on  a cathode  ray  tube  (CRT).  This  display 
simulates  a three  dimensional  picture  of  the  mother  ship's  stern  well  to 
facilitate  docking.  Range  and  bearing  and  a target  showing  relative  size 
according  to  distance  are  shown  to  the  pilot. 

No  attempt  is  made  here  to  describe  the  method  or  development  of  tiiis 
visual  simulation.  The  display  was  developed  and  programmed  at  NTEC 
and  interfaced  by  CSDL  to  the  real-time  simulation. 

Although  an  attempt  to  display  a mother  diip  is  made  and  docking  maneuvers 
can  be  performed,  no  modelling  of  this  complex  two  body  problem  has  been 
considered. 


Another  Visual  Simulation 

As  a part  of  the  ACV-LC  trainer,  another  visual  simulation  has  also  been 
installed.  This  visual  simulation  needs  signals  from  the  simulation.  This 
part  of  the  trainer  was  developed  and  interfaced  to  the  real-time  simulation 
by  NAVTRAEOUIPCEN.  No  attempt  at  description  is  made  here. 
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INITIALIZATION 

In  order  to  begin  execution  of  the  mathematical  model  equations,  certain 
variables  must  be  initialized.  Since  this  model  was  developed  around  a quasi- 
operating  point  of  the  vehicle,  no  startup  procedures  have  been  defined.  Initial- 
ization of  the  vehicle  to  this  rtate  must  be  performed.  The  state  chosen  for 
initialization  is  a hovering  condition  at  zero  speed  with  the  vehicle  supported 
by  its  air  cushion. 

The  individual  parameters  to  be  initialized  are  all  encompassing,  but  their 
values  were  chosen  by  trial  and  error.  Other  parameters  that  need  to  be 
initialized  are  calculated  by  their  respective  equations  and  by  properly 
arranging  the  flow  of  execution. 

The  parameters  which  need  to  be  initialized  are  simply  the  effectors  and  the 
pressures.  This  establishes  the  proper  qua  si -operating  condition  of  the 
vehicle.  Correspondingly,  however,  the  pilot  commands  are  initialized  on 
a one  to  one  basis  with  the  actual  control  device. 

Subroutine:  SYSACVIC 


^Tn  ^ ^Tn<#) 


^ANj  ^ ^Sj  ^^FAN 
^ ^FANj^^p 
^CUSH  i " ^CUSH  i4> 


Pj 

*MANj 


= P 


= P 


MANjip 
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i=l 

i=2 

i=3 

i=4 

P 


CUSH  i 


CUSH  i4> 


N, 


Tn 


N, 


Tn<^ 


1,  n=l,  2,  3 Starboard 

2, n=4,  5,  6 Port 

STBD  Fwd  cushion  compartment 
STBD  Aft  cushion  compartment 
PORT  AFT  cushion  compartment 
PORT  FWD  cushion  compartment 

Pressure  within  cushion  compartment  (Lbs/Ft^) 

Initial  value  of  Pqusjj  j = 109.  37297  Ibs/ft^ 
Turbine  speed  (RPM) 

Initial  turbine  speed  = 16,  800  RPM 


N, 


Sj 


Power  shaft  speed  (RPM) 


No- . 
^FANj 

Np. 

^MANj 

^MAN3<^ 

% 

^%AN 

GRp 


Initial  power  shaft  speed  = 13, 200  RPM 
Fan  speed  (RPM) 

PropeUer  speed  (RPM) 

Pitch  angle  of  propeller  (Deg) 

Initial  pitch  angle  of  propeller  = 12.  9° 
Manifold  pressure  (Lbs /Ft  ) 

Initial  manifold  pressiire  = 130. 9 Ibs/ft^ 
Thrust  nozzle  angle  (Deg) 

Initial  thrust  nozzle  angle  = 180° 

Power  shaft  to  fan  gear  ratio=  . 1297 
Fan  to  propeller  gear  ratio  = . 6427 


i 
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Subroutine:  CMDINIT 


%cw 

= 0 

^NC 

= 0 

^RC 

= 0 

^NS 

= 0 

^PjCL 

^TCn 

^SC3 

^PCW 

^NC 

^RC 


NS 


^PjCL 

^3<P 


N, 


TCn 


N, 


N, 


Tn«#- 

SCj 


1, n=l,  2, 3 Starboard 

2, n=4, 5, 6 Port 

Propeller  pitch  vernier  wheel  command  (Deg) 
Nozzle  an^e  command  (Deg) 

Rudder  angle  command  (Deg) 

Nozzle  an^e  switch  (Fwd/Aft) 

Propeller  pitch  angle  command  (Deg) 

Initial  propeller  pitch  angle  = 12.  9° 

Turbine  speed  command  (RPM) 

Initial  turbine  speed  = 16, 800  RPM 

Power  shaft  speed  command  (RPM) 

Initial  power  shaft  speed  = 13, 200  RPM 
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SECTION  IV 
SOFTWARE  STRUCTURE 

The  software  structure  is  a modular  system  and  is  basically  comprised  of  5 
major  divisions: 


Program  load  and  initialization 

Backgroimd 

Foregroimd 

Trap  handling 

Common  storage 

The  first  three  tasks  are  linked  directly  by  program  flow  and  shown  in 
Figure  44  . Traps  are  handled  for  normal  and  abnormal  occurrences  in 
these  three  tasks.  The  separate  common  storage  (module  DATAPOOL) 
contains  data  (constants  and  variables)  which  is  used  or  manipulated  by 
the  4 above  tasks. 


LOAD  AND  INITIALIZATION 

The  ACV-LC  real-time  simulation  program  is  loaded  from  9 track  magnetic 
tape.  Once  loaded  program  execution  begins  with  all  interrupts  inhibited. 

The  traps  handled  by  the  program  are  initialized  and  are; 

Control  panel 

Counter  3 zero 

Counter  3 

Power  up 

Power  down 

CAL4 

Memory  parity. 

Non-allowed  operation 
Unimplemented  instmction 
Push  down  stack 
Fixed  point  aritiunetic  overflow 
Floating  point  fault 

Decimal  arithmetic  fatilt  ^ 

1 s 
' f 

Watch  Dog  timer  runout 
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PROGRAM  ENTRY 


READVGW 


Cue  Typewriter  for  9T  Unit  # 


Inhibit  Interrupts 

I 

Initialize  T raps 


READVGWI 


i 

Read  VGW  Table 

Initialize  I/O 

Initialize  Timing 


EXECACVl 


ACVSIM 


1 


BACKGROUND 


FOREGROUND 

1 1 


SIMULATION 


RDVGWBPM 

INITSTRT 


Figure  44  . Basic  Software  Flow  and  Structure 
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Restart 

CALI 

The  vehicle  generated  wave  table  of  Kernel  values  is  read  if  the  option  was 
selected  at  the  time  of  load  tape  generation  and  if  the  proper  response  is 
made  to  type%vriter  cue.  The  table  is  read  from  the  tape  and  stored  in 
sequence.  The  beginning  addresses  of  each  subtable  are  stored  in  a 
separate  table. 

Initialization  of  aU  Input/Output  (I/O)  is  performed.  All  I/O  software 
request  flags  are  reset.  The  time  of  day  and  the  date  are  reset  and  the 
type-in  flag  is  set  to  request  the  proper  time  and  date.  The  I/O  hardware 
and  the  audible  alarm  are  reset  to  their  usable  states. 

The  program  timing  initialization  is  performed  by  resetting  and  disabling  all 
interrupts  and  by  setting  the  Counter  3 interrupt.  The  interrupts  used  are 
enabled: 


Counter  3 
Counter  3 Zero 
Memory  Parity 
Control  Panel 


Program  Entry 

At  the  time  of  self  loading  tape  creation,  a transfer  address  is  specified 
which  ^vill  be  the  point  of  initial  execution  after  program  load.  There  are 
3 beginning  points  from  which  the  simulation  program  can  start  (see 
Figure  44): 


READVGW 

READVGWI 

INITSTRT 


READVGVV 

If  READVGW  is  selected  as  the  program  entry,  this  message  is  typed  on 
the  'console  typewriter: 

TYPE  9T  UNIT  NO.  FOR  VGW  TABLE 

OR  NL  TO  NOT  READ  TABLE: 

The  typewriter  input  light  is  turned  on  and  the  program  %vaits  for  an  input 
by  the  operator.  This  input  must  be  a decimal  digit  and/or  a nevr  line(NL) 
character. 
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If  a decimal  digit  I and  NL  is  typed,  the  vehicle  generated  wave  kernel 
table  will  be  read  from  9 track  (9T)  magnetic  tape  <MT)  xmit  ntimber  I. 
Transfer  will  be  made  to  the  program  entry  REIADVGWI. 

If  an  NL  character,  only,  is  typed,  the  vehicle  generated  wave  table 
will  not  be  read  and  the  program  will  transfer  to  the  program  entry 
miTSTRT. 


READVGWI 

if  READVGWI  is  selected  as  the  program  entry,  the  9T-MT  unit  number 
equals  0 and  the  vehicle  generated  wave  kernel  table  will  be  ret.d  from 
the  end  of  the  self  loading  program  tape. 

When  the  kernel  table  is  read  successfully,  these  messages  are  typed  to 
the  operator: 


VGW  TABLE  SUCCESSFULLY  READ 
NO.  OF  SUBTABLES  READ  = 82 

The  flag,  VGWFLAG,  will  be  set  equal  to  1 to  run  the  VGW  package  and 
the  program  will  transfer  to  the  program  entry  INITSTRT. 


INTTSTRT 

If  the  vehicle  generated  wave  kernel  table  has  not  been  read  by  the  time 
this  program  entry  has  been  reached,  the  VGW  package  will  not  be  run 
in  simulation.  The  program  entry  INITSTRT  initializes  the  traps,  I/O, 
and  the  timing  and  tranters  to  the  Backgroxmd  Task  (BACKGRND)  for  system 
operation. 


BACKGROUND 

The  Background  Task  takes  up  all  computer  time  not  spent  in  the  Foreground. 
The  primary  function  of  the  Background  is  to  perform  I/O.  The  several 
functions  are:  typewriter  input/output,  card  reader  input,  line  printer  ou^ut, 
and  the  CRT  display.  Other  functions  are  also  performed;  the  calculation  of 
dwell  time  (amount  of  time  spent  in  the  background)  and  the  vehicle  generated 
wave  hei^  profile  calculations. 

Typewriter  input  allows  the  operator  to  change  the  contents  of  any  memory 
location,  the  time  of  day,  and  the  date.  Typewriter  output  displays  the 
coiitents  of  any  memory  location  in  any  of  several  formats.  The  typewriter 
may  also  display  software  status  messages  to  the  operator.  Cards  may  be 
read  to  perform  modification  to  data  within  memory.  The  line  printer  may  be 
cued  to  display  the  state  of  the  ACV-LC  in  a clear  format  or  the  contents  of 
the  DATAPCXJL  may  be  dumped. 
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PROGRAMMING.  Figures  45  through  49  suow  tne.  ilow  of  module  execution  of  the 
Background  Task.  These  flow  diagrams  are  functional  in  nature  and  do  not  show 
the  minute  programming  details.  The  iiiajor  modules  iiuiking  un  the  Background 
Task  are: 

BACKGRND 

TYPE 

TYPE  MSG 

CARDS 

CARDIO 

PRINT 

These  modules  are  cycled  indefinitely  and  are  subject  to  interruption  by  traps, 
but  primarily  by  the  Foreground  Task.  Other  modules  are  called  as  needed 
to  perform  specific  functions  (i.  e. , CARDREAD(M:READ),  PRNTLINE(M:PRINT). 

Dwell  Time 

The  background  timer  is  a counter  that  keeps  track  of  the  length  of  time 
that  the  background  task  is  in  operation.  The  value  of  the  counter  is  saved 
and  reset  each  cycle  by  the  foreground.  Subtracting  the  background  time 
from  the  time  of  a complete  ACV  simulation  sampling  period  will  provide 
the  time  ^ent  in  the  foregrotmd.  The  timer  (BGTIME)  is  incremented  in  units 
of  12  microseconds.  It  is  saved  in  location  BGTIMSAV.  With  the  background 
timer  acti^,  the  peripheral  equipmeit  and  the  CRT  and  the  VGW  cannot 
operate.  Tme  background  timer  occupies  the  entire  time  spent  in  the  Backgroimd 
Task. 

Switches 

There  are  2 sets  of  operator  controlled  switches  which  control  the  simulation 
program,  the  console  sense  switches  (CSW)  and  the  system  switches  (SSW). 

The  switches  and  their  functions  in  the  Background  are: 


CSW  1 deactivate  VGW  & CRT 

3 deactivate  CRT  display 

SSW  8 Read  cards 


call  SYSPRINT 


dwell  time 
CaU  PRTDATA 


s 
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Figure  46.  Background  Task  Flow  - TYPE 
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Figure  47  , Background  Task  Flow-TYPEMSG 
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CARDS 


PRINT 

Figure  48.  Background  Task  Flow-CARDS 


194 


1 


NAVTRAEQUIPCEN  75-C-0057-1 


BACKGRND 


Figure  49.  Background  Task  Flow-PRINT 


VGWAVE 

SYSPRINT 

PRTDATA 

DLOOPER 

This  module  does  the  necessary  data  conversion  to  be  used  for  generation 
of  the  CRT  display  and  calls  DLCXDP.  DLCXDP  generates  the  CRT  display. 


VGWAVE 


This  module  calculates  the  vehicle  generated  wave  height  at  one  point  under  the 
hull  planform.  Successive  calls  calculate  the  wave  height  at  the  next  point. 

SYSPRINT 

This  module  prints  the  ACV  state  vectors,  forces  and  moments,  and  the 
status  of  effectors  at  their  current  values.  This  printout  is  formatted  for 
readability. 

PRTDATA  * 

This  module  dumps  the  DATAPOOL  contents  onto  the  line  printer  by 
hexidecimal  address  in  floating  point  short  (FS)  format. 

FLAGS 

CRTFLAG  is  a software  flag  to  activate  the  CRT  display.  This  flag  is 
set  in  the  foregroxmd. 

VGGOFIAG  is  a software  flag  to  execute  the  VGW  profile  eauations.  This 
flag  is  set  in  the  foreground. 


i 

i 


! 


i 

I 

! 

! 

i 

j 

I 

1 

1 

1 

\ 


196 





NAVTRAEQUIPCEN  75-C-0057-1 


FOREGROUND 

The  Foreground  Task  is  activated  from  the  Counter  3 Zero  interrupt  at 
regular  intervals.  This  task  updates  the  time  of  day,  saves  the  computed 
•dwell  time  from  the  Background  and  then  resets  the  dwell  time  counter. 

The  interrupt  is  reset  and  if  execution  in  the  foreground  exceeded  the 
proper  interval,  an  error  indicator  is  incremented.  If  the  error  occurred, 
immediate  return  is  made  to  that  point  of  interruption  for  processing 
continuation.  The  previous  cycle  will  be  oo  mpleted  and  the  new  cycle 
will  be  lost.  Normally,  this  ei  i‘or  condition  would  not  occur  and  the 
simulation  -would  be  csiUed.  However,  if  console  sense  switch  (CSW)  1 is  set, 
the  simulation  executive  would  not  be  called.  Return  to  the  Background  is 
made  at  the  conclusion  of  the  Foreground.  See  Figure  50  for  Foreground 
Task  Flow. 

The  foreground  task  is  primarily  intended  to  provide  the  timing  for  the 
simulation  of  the  ACV-LC.  This  timing  is  necessary  to  provide  the  proper 
time  intervals  for  the  mathematical  modelling  techniques  of  real-time 
simulation.  Although  the  software  description  to  this  point  has  only  been 
concerned  -with  the  timing  structure  and  overhead,  the  primary  objective 
of  simulation  is  accomplished  in  the  Foreground  Task.  The  mathematical 
modelling  of  the  ACV-LC  is  controlled  by  the  simulation  executive  ACVSIM. 
Consequently,  operation  in  real-time  necessitates  a call  to  ACVSIM  at  every 
time  period. 

The  Counter  3 Zero  Interrupt  occurs  at  regular  intervals  of  50  msec  (20Hz). 

This  interriqjt  activates  the  Foregroimd  Task  and  is  handled  by  the  routine 
EXECACVl.  The  timing  of  the  foreground  occimrence  cai  be  altered  by 
changing  the  variables  RUNFREQ  and  EXCOUNTI.  RUNFREQ  is  the  integer 
value  of  the  frequency  (20  Hz).  EXCOUNTI  is  the  integer  value  -8000 /RUNFREQ 
(-400).  In  order  that  real-time  be  realized,  RUNFREQ  should  divide  into 
8000  an  integer  number  of  times  exactly.  EXCOUNTI  is  used  to  reset  the 
Counter  3 Interrupt.  RUNFREQ  and  EXCOUNTI  are  the  only  variables  used  by 
the  program  to  determine  -the  foreground  timing.  However,  the  sampling  period 
for  real-time  simulation  should  correspond  to  the  interrupt  period  (1 /RUNFREQ). 
DELTAT  is  a floating  point  variable  and  is  the  sampling  period  used  in  the 
simulation  mathematical  model.  The  value  equals  . 05  sec. 


ACVSIM 

The  simulation  executive  (ACVSIM)  controls  the  execution  of  the  mathematical 
model.  (See  Figure  51).  Upon  entry,  the  CRT  display  flag  (CRTFLAG)  and 
the  VGW  execution  flag  (VGGOFLAG)  are  reset.  Simulation  moding  switches 
are  tested  to  determine  the  path  of  execution: 

SSW  1 initial  conditions  (IC)  or  reset 

2 hold  or  freeze 

3 operate  or  run 
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EXECACVl 


Figure  50.  Foreground  Task  Flow 
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ACVSIM 


Figure  51.  Simulation  Exectrtive  Flow-ACVSIM 
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Any  combination  of  mere  than  1 mode  s\vitch  on  or  all  off  results  in  the 
hold  mode,  which  is  a donothing  mode. 

The  reset  mode  causes  initialization  of  the  ACV-LC  to  a stable  condition 
at  zero  speed  and  specified  position.  (See  Figure  52  for  flow^  The  time 
overflow  indicator  cannot  be  incremented  and  the  cockpit  displays  are 
updated  to  correspond  to  the  initial  conditions. 

The  operate  mode  is  the  read -time  simulation  of  the  ACV-LC.  The  vehicle  is 
operational  and  maneuverable  in  this  mode.  Each  cycle  in  the  operate  mode 
first  tests  a flag  (CIND)  for  a collision  with  the  mother  ship.  If  a collision 
has  occurred,  the  hold  mode  is  entered.  Second,  execution  time  (time  since 
entering  operate  from  IC)  is  updated.  Third,  the  input  discrete  word  is  read 
to  determine  coclq>it  switch  status.  Fourth,  if  execution  time  is  less  than 
10  seconds  or  if  console  sense  switch  (CSW)  2 is  set,  the  pilot  commands 
are  bypassed.  If  execution  time  is  greater  than  10  seconds  or  CSW  2 is 
reset  (the  normal  case),  then  the  pilot  commands  from  the  cockpit  are 
read  into  the  simulation  for  use  in  the  model  equations.  Fifth,  if  CSW  3 
is  not  set,  CRTFLAG  is  incremented  to  enable  the  CRT  display.  Sixth, 
the  mathematical  model  equations  (OPB^NS,  see  Figure  53  )*are  called 
for  exectition,  and  then  the  displays  in  the  cockpit  are  updated. 

The  variable  MODEWORD  is  set  to  a value  of  characters  corresponding  to  the 
mode  of  operation; 

*INIT'  IC,  reset 

‘HOLD'  hold,  freeze 

'CRSH'  hold,  mother  ship  collision 

'OP*  operate  without  pilot  commands 

*OPGO*  Qperate  normally 

This  set  of  characters  when  printed  is  an  immediate  indicator  of  system  status 
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ICEQNS 


Figure  52.  Initialization  Equations  Flow-  ICEQNS 
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OPEQNS 

l:all  CVIHULL 
OCEAN 


VGWFLAG=^ 


VGGOFLAG=VGWFLAG  i 


Call  WINDPAR 
EFPOS 
TURBINES 
FILTER 
HULLHT 
SKAREA 
VOLCUSH 
WVPUMP 
CUSHP 

airflow 

SEA 

CUSH 


VGWFLAG^ 


Call  SKIRT 


Call  SKIRTVG 


Call  AIR 
AERO 
DUCTS 
DAMP 
RUDDERS 
NOZZLES 
PROPS 
INTGRT 
HISTORY 


Figure  53.  Operational  Equations  Flow  - OPEQNS 
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Traps  cause  an  interruption  in  program  flow.  Returning  to  flow  continuation  is 
optional  and  depends  upon  its  implementation  (and  operation,  in  some  cases). 
The  Foreground  has  already  been  discussed,  but  is  a trap  (Counter  3 Zero) 
interrupting  the  Background  and  does  return  for  continuation.  Other  traps  are 
activated  for  a simple  task  and  return  immediately  for  continuation: 

Counter  3 

Control  Panel 


The  Power  Up  interrupt  causes  an  unconditional  branch  to  the  entry  INITSTRT. 

The  Power  Dovm  interrupt  causes  an  infinite  donothing  lo.op.  Restart  is  not  really 
a'  trap',  but  simply  a means  of  initializing  and/or  restarting  the  simulation  program. 
When  the  operator  pushes  SYSTEM  or  CPU  RESET  the  computer  initializes  itself 
to  memory  location  X'26'.  If  the  operator  pushes  RUN  an  unconditional  branch 
to  entry  INITSTRT  is  caused. 

CALI  is  a software  caused  trap  which  is  used  to  cue  I/O  procedures.  CALI 
is  only  used  for  those  procedures  involving  the  reading  of  the  V.GW  table.  This 
trapTT  activated  to  perform  an  I/O  function  and  then  returns  to  the  program 
for  continuation. 

All  abnormal  errors  are  hardware  cued  and  are  caused  by  software  and/or 
hardware  faults.  These  faults  are  handled  by  traps.  If  an  abnormal  trap 
occurs,  the  CAL4  trap  is  activated  (see  Figure  54  ) to  record  the  error, 
turn  on  the  audible  alarm,  and  then  places  the  simulation  program  into  a WAIT 
state  for  operator  intervention.  The  fault  traps  which  are  handled  in  this  manner 
are: 

0 Memory.parity 

1 Non -allowed  operation 

2 Unimplemented  instruction 

3 Push  down  stack 

4 Fixed  point  arithmetic  overflow 

5 Floating  point  fault 

6 Decimal  arithmetic  fault 

7 Watch  Dog  timer  runout 

The. error  number  is  recorded  in  location  C4ERR#  and  the  address  of  the 
error  is  recorded  in  C4ERRLOC,  If  the  Memory  Parity  trap  occurred, 
the  memory  location  which  has  incorrect  parity  is  recorded  in  C4ERRMEM. 
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CAL4  SAVE=  All  Registers 
Deactivate  Interrupts 
C4ERR#=  Error  Number 
C4ERRL0C=  Address  of  Error 
C4ERRMEM=  Address  of  Memory  Parity 


Figure  54  . 


CAL4  Error  Processing  FI 
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COMMON  STORAGE 

The  DATAPOOL  acts  as  common  storage  for  the  complete  simulation 
program.  Most  modules  interface  with  the  DATAPOOL  in  order  to  facilitate 
transfer  of  data  from  one  module  to  another.  This  centralized  location  of 
data  is  also  for  the  operator's  convenience  and  is  intended  to  minimize 
external  definitions  of  variables  within  program  modules. 

The  DATAPOOL  is  formatted  for  readability  in  the  following  vtay; 

P’orces  and  Moments 
Summation  Forces,  Moments 
State  Vectors 
Transformation  Matrix 
State  for  Line  Drawing 
Integration  Increments  (Not  Used) 

Change  in  Inertial  Position 

Constants  about  CG 

Environmental  Wind  Factors 

Special  System  Parameters 

Equations  of  Motion 

Propellers 

Thrust  Nozzles 

Rudders 

Aerodynamics 

Propeller  Ducts 

Skirt 

Manifolds 

Compartment  Volumes 
Wave  Pumping 
Cushions 
Gap  Areas 

Gas  Turbines  and  Power  Shafts 

Seaway 

VGW 

Wave  Heights 

Hull  Point  Locations 

Height  of  Hull  above  Water 

Skirt  Clearances 

Input  Discretes 

Analog  Inputs 

Analog  Outputs 

I/O  Parameter  List 

Typewriter  Status  Messages 
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MODULARITY  AND  TIMING 

As  has  already  been  stated,  the  simulation  program  consists  of  a itiodular 
breakdown  by  function.  That  is,  each  module  does  a specified  function  or 
task  mthin  the  system.  Most  modules  are  executed  as  subroutines  with 
calls  to  them  and  returned  from  them.  The  structure  has  already  been 
shown,  so  that  it  is  obvious  that  EXECACV  and  READVGW  are  entered  by 
uncon^tional  branches.  The  major  breakdown  of  the  Background  Task 
has  transfers  from  major  modules  by  unconditional  branches. 

All  subroutines  link  by  using  Register  15  except  those  whose  names  begin 
with  a digit  (See  Table  10).  Subroutines  v/hich  have  names  beginning  with 
a digit  link  by  Register  6 and  are  library  functions.  The  I/O  library 
functions  link  by  Register  15. 

The  digit  6 preceding  the  name  signifies  that  the  module  was  taken  from  the 
Sigma  7 FORTRAN  library.  The  digit  7 preceding  the  name  signifies  that  a 
special  numerical  technique  was  developed  for  rapid  execution  as  a library 
function.  The  digit  1 preceding  the  name  signifies  that  the  module  was 
specially  developed  to  solve  a particular  numerical  task  within  the  mathem.^tical 
model  and  are  not  library  functions  bvit  use  Register  6 as  the  linkage. 

Secondary  entry  points  within  the  major  modules  of  Table  10  follow  the 
same  rules.  See  Table  11  for  the  list. 

The  timing  of  most  real-time  subroutines  has  been  indicated  in  the  tables.  These 
numbers  were  obtained  under  the  monitor  and  are  not  necessarily  accurate  as  the 
monitor  does  not  always  respond  at  the  proper  time.  Nevertheless,  these  timing 
numbers  are  good  indicators  of  relative  timing  as  they  were  obtained  by 
executing  each  subroutine  1000  times.  The  timing  for  the  library  functions  was 
obtained  by  either  taking  the  timing  indicated  from  the  FORTRAN  library  listings 
or  by  counting  instructions. 

Individual  timing  of  each  subroutine  is  not  important  by  itself,  but  the  proper 
sum  in  the  operate  mode  is.  From  the  table,  the  timing  of  OPEQNS  is 
29.372  m sec.  This  number  represents  the  system  off-line  without  FILTER 
and  without  the  real-time  simulation  overhead  of  mode  determination,  reading 
A/D's  and  discretes,  and  outputting  D/A's.  From  actual  dwell  time  of  the 
real-time  system,  29.  2 msec  represents  the  time  spent  executing  equations  in 
the  deep  ocean  ( 29.2/50=  58%)  and  40.5  msec  represents  the  time  spent  in  the 
ocean  shoaling  region  (40.  5/50=  81%). 
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Major  Simulation  Modules 
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TABLE  10.  Major  Simulation  Modules 
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OPERATIONS 


SYSTEM  DESCRIPTION 

The  ACV  Simulation  can  logical  Iv  be  de.^cribed  in  three  .sectiop.s:  the  real- 
time  equation.s,  the  background  l/O,  and  the  traps. 

The  largest  section  of  core  is  taken  b;  the.  real-time  .\CV  ecuatioip  . The 
timing  for  these  equations  is  provided  by  the  Counter  Cero  Lnterruii.  Ai 
present  the  Count  3 Zero  Interrupt  is  pr.-grarired  te  run  at  20  it::,  '^he  pro- 
cessor for  the  Counter  3 Zero  Interrupt  is  contained  in  the  urogram  hXEC.ACV 
which  contains  all  the  ACV  interrupt  processors.  The  Counter  3 Zero  Interrupt 
calls  the  simulation  program  ACVSIM.  ACVSIM  drives  ali  tiie  ACV  r ath  mod-.'i 
programs. 

The  timing  of  the  ACV  simulation  can  be  changed  from  20  Ha  bv  altering  the 
locations  RUNFREO,  EXCOUNTI,  and  DEl.TAT.  RUNFREQ  is  the  integer  value  of 
the  frequency  which  is  20  at  present.  EXCOliNTJ  is  -8000.  RiJNF.RECi  and  i.s  - iLD 
at  present.  EXCOUNTI  is  used  to  reset  the  Counter  3 Interrupt.  Dlii.TAT  is 
a floating  point  single  value.  It  is  the  integration  time  constant;  at 
present  set  to  .05  (1/20).  In  order  not  to  lose  time  in  the  real  time  simu- 
lation, RUNFREQ  should  divide  into  8000  an  integer  number  of  times. 

T!ie  ACV  system  has  been  programmed  to  . llov’  the  real  time  equations  to  inter- 
rupt themselves.  If  this  occures,  the  equations  will  not  be  re.started  during 
the  new  cycle.  The  last  cycle  will  be  completed  from  the  point  of  inter- 
ruption and  the  rest  of  the  cycle  will  be  spent  in  the  background  loop. 

This  event  will  be  recorded  by  the  Counter  3 Zero  Interrupt.  The  interrupt 
will  increment  the  location  EXSEI.FIN  each  time  it  interrupts  the  real  tire 
programs . 

« 

TKe  background  I/O  programs  are  continually  being  cycled  through  whenever 
the  real  time  equations  are  not  being  executed.  The  background  programs 
are  made  up  of  six  subsections.  These  subsections  are;  the  background  timer, 
the  CRT  driver,  the  typewriter  input  processor,  the  tyj ewriter  error  message 
processor,  the  card  reader  T-rocessor,  and  the  line  printer  processor.  The.st 
sections  contain  all  the  I/O  for  the  respective  peripheral  equipment. 

The  background  timer  is  requested  by  the  operator  by  setting  .sy.stem.  stnse 
switch  ll(SSWll).  The  background  timer  is  a counter  that  keeps  track  of 
the  length  of  time  that  background  is  in  operation.  The  value  of  t!ic  counter 
is  saved  and  reset  each  cycle  by  the  real  time  interrupt.  Subtracting  the 
background  time  from  the  time  of  a complete  ACV  cycle  will  provide  the  tire 
that  the  real  time  equations  are  using.  The  timer  is  incremented  in  unit.s 
of  12  microseconds.  It  is  saved  in  location  BGTIKSAV.  As  long  as  the  back- 
ground timer  is  being  calculated  (i.e.,  SSWll  is  set),  the  perinheral  equif- 
ment  and  the  CRT  cannot  be  used.  The  background  timing  program  occupies 
•the  entire  tine  the  simulation  is  in  the  background  state. 

The  CRT  driver  will  display  visual  cues  to  the  cockpit  operator  during  t'.u 
ACV  simulation.  It  will  not  be  executed  if  console  sense  switch  3(CSi'3) 
is  put  on. 
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Ine  typewriter  input  processor  reads  and  r recesses  opf ratnr  re; uc‘  l*  to 
the  typewriter.  Requests  can  be  made  tc*  the  typewriter  rrot  ra-  a.-;  J-  ng  as 
the  typewriter  input  light  is  lit.  This  light  is  set  on  dun  nr.  syntein  ini- 
tialization. It  can  be  put  on  or  ofi  alternately  by  triggering  the  console 
interrupt.  Ivhen  type  requests  are  not  being  made,  this  light  siiould  be  put 
off. 

Th.e  typev/riter  error  message  processor  will  print  error  messages  on  tlie 
typewriter.  These  messages  c.annot  be  typed  if  the  typewriter  input  light 
is  lit. 

The  card  reader  processor  reads  cards.  Card  reading  is  requested  bv  the 
operator  by  setting  SSW8. 

The  line  printer  processor  prints  syscera  parameters  onto  the  line  urinter. 
This  is  requested  by  setting  SSV’  9. 

The  third  section  of  the  ACV  simulation  is  the  error  trap  prcce.‘;sor.  vrnen 
an  error  trap  occurs,  the  trap  in  turn  calls  the  trap  CAi,  A.  CAL  will 
save  ail  the  register  values  plus  some  diagnostic  information.  Seven  types 
of  errors  can  occur.  (Memory  pari tv  is  not  really  a trap  but  is  treated 
1 Ike  one  in  the  error  routine) . 


Error  V Error 

0 Memory  parity 

1 _ Non-allowed  operation 

2 Unimplemiented  Instruction 

3 rush  down  stack  fault 

A Fixed  point  arithm.etic  overflow 

5 rioating  point  fault 

6 . Decimal  arithm.etic  fault 

7 Vlatch  dog  timer 

Three  words  are  saved,  along  with  the  data  registers,  when  one  c£  these 
errors  occurs.  CA  ERR?  contains  the  number  of  the  error  that  occurred  in 
Che  right  most  hexadecimal  character.  CA  ERRLOC  is  the  PSl’l  of  the  L’oD  o£ 
the  trap  that  occurred.  This  word  contains  the  memory  location  that  caused 
the  error  (except  memory  parity  error)  and  the  condition  code  bits  that 
further  define  the  error.  CA  ERRMEM  has  the  memory  module  fault  bits  in 
bits  2A-31.  If  an  error  occurs  the  console  alarm  is  triggered.  The  com- 
puter will  go  into  the  wAlT  state.  There  are  three  ways  to  recover  from 
the  W.AIT  condition.  The  program  will  continue  from  the  location  chat  caused 
■Che  error,  if  the  operator  sets  console  sense  switch  1 (CSR’l)  off  and  puts 
the  computer  in  the  IDLE  and  then  RUM  modes.  The  system  will  be  reinitialized 
if  the  operator  sets  CSWl  on  and  puts  the  comruter  in  the  IDLE  and  RUN  state. 
An  alternate  method  of  reinitializing  the  system  can  be  done  by  putting  the 
computer  in  IDLE,  pressing  CPU  RESET,  and  putting  the  computer  in  the  RUK 
mode.  The  last  method  will  reinitialize  the  ACV  system  at  anv  time  during 
the  .ACV  Simulation. 
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lAL  LOAD 


The  ACV  Similation  Systen  is  contained  on  R^f:netic  tape-  Boot  the  < v 
into  the  coE^puter  with  console  sense  switch  1 (CS«‘J)  on  and  all  thi 
sense  switches  (SSW)  off.  CSUl  inhibits  the  real  tire  iiiCerrupt  fror 
calling  the  prograr  ACVSIM.  CSV;  1 is  in  eitect  an  alternate  frec-xo  ">odr-. 


The  typewriter  will  print: 

TYPE  IN  TIME  (HH:Mri)  AND  DATE  (MOM  DD.’YV) 


and  the  typewriter  input  light  will  be  rut  on.  The  tiise  and  date  should 
be  typed  in.  The  simulation  will  run  if  these  values  are  not  entered, 
however . 


The  format  of  the  time  input  is: 

TIME  KH:{94 

KH  is  the  time  in  hours  (0-23)  and  Hy*  is 
spacing  is  free  format  except  there  must 
between  HH  and  121.  (;  is  optional.  Any 

space) . 

The  format  of  the  date  input  is: 

BATE  MON  BD.’YY 

HON  is  the  three  character  abbreviation  of  the  current  month  (e.g.,  JAN, 
FEB,  MAR,  etc.)  DD  is  the  day  of  the  month  (01  to  31);  ’YY  is  the  last 
two  digits  of  the  year  (e.g.,  *75).  The  spacing  of  the  data  input  is  fre 
format  but  the  format  of  12)N  DDj'YT  is  fixed.  There  must  be  one  and  only 
one  space  between  the  month  and  the  day. 

Unless  other  type-in  requests  are  cade,  the  typewriter  input  light  should 
be  put  off  by  activating  the  console  interrupt.  Other  type-in  requests 
will  be  discussed  below. 


the  tine  in  minutes  (0-59) - The 
be  one  and  only  one  character 
character  will  do,  including  a 


The  ACV  simulation  is  controlled  bv  the 


stem  sense  switches  1 
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SIMULATION  OPERATION 

The  ACV  simulation  is  controlled  by  the  svstem  seir-e  iW-Lciic-  j . . i.s.  1. 

SSWl  is  on  for  RESET  mode 
SSW2  is  on  for  FREEZE  mode 
SSU'3  is  on  for  OPERATE  mode 

One  and  only  one  switch  should  be  on  at  any  one  time.  Any  other  ccmbLna- 
tion  (two  or  three  of  the  SSWs  on  or  all  of  them  off)  of  switches  1,  2,  and 
3 will  put  the  simulation  in  the  FREEZE  rode. 

When  CSWl  is  on  (as  given  in  the  load  rroccdure  above)  none  oi  t!i..  -;y.:lcr' 
sense  switches  are  tested.  Their  positions  at  this  time  are  imi.aterial. 

To  start  the  ACV  sim.ulation,  put  the  system  into  RESET  by  nutting  SyV.'i  on 
and  SSW3  off.  Put  CSWl  off.  CSW2  and  CSV.'3  should  also  he  off  at  this 
time.  The  computer  is  now  continually  resetting  the  simulation  parameters. 

To  put  the  ACV  simulation  in  OPERATE  set  SSWl  off  and  SSW3  on.  The  FREEZE 
mode  will  be  entered  momentarily  until  the  p'  per  moding  is  achieved. 

To  freeze  the  simulation  SSi.’2  should  be  set  to  on.  The  state  of  the  other 
moding  switches  is  not  important.  FREEZE  mode  is  in  effect  a S;  steti'  nc-np. 

No  simulation  I/O  or  system  equations  will  be  executed.  I/O  on  peripheral 
equipment  such  as  the  typewriter  can  be  done,  however. 

CSWl,  CSW2,  and  CSVJ3  are  used  to  inhibit  different  parts  of  tlie  ACV  system. 
CSWl  inhibits  the  entire  simulation  (ti.e  call  to  ACVSJM  and  the  CP'i  driver). 
CSW2  inhibits  the  commands  from  the  cockpit  to  the  simulation.  CS'.'3  inhibits 
the  call  to  the  CRT  driver.  CSW2  and  CSW3  are  only  effective  in  OPEPJ'TE 
mode  because  the  cockpit  parameters  and  the  CRT  driver  (DLOOP)  are  i.ot 
referred  to  in  RESET  or  FREEZE  modes.  The  cockpit  input.s  and  CRT  outputs 
will  remain  at  their  last  values  until  RESET  is  entered.  Thev  wi]i  be  re- 
initialized at  this  time.  CSWl,  CSW2,  and  CSW3  do  not  effect  the  operation 

of  the  typewriter,  the  card  reader,  or  the  line  printer. 

The  status  of  the  ACV  simulation  can  be  monitored  periodically.  Momen- 
tarily setting  SSW9  prints  various  simulation  parameters  on  the  line  printer. 
IF  SSW9  is  left  on,  the  data  will  be  printed  continuously.  A description 
of  the  printed  data  can  be  found  in  ACV  System  Printout  Description.  Printing 
can  be  done  in  any  simulation  mode  but  the  printed  data  would  be  more  useful 
if  printing  is  done  in  the  FREEZE  mode  only.  Printing  the  state  variables 

cannot  be  done  if  the  background  timer  is  in  operation  (SSWll  is  on)  or  if 

card  reading  is  taking  place. 

The  ACV  system  has  the  capability  of  altering  or  printing  tlie  contents  of 
any  location  in  memory.  Locati''ns  can  be  altered  with  either  the  card 
reader  or  the  typewriter.  Locations  can  be  printed  on  the  typewriter. 

The  capability  to  print  random  memory  location  is  an  additional  feature 
separate  from  the  line  printer  status  printout.  The  typewriter  and  card 
reader  I/O  are  Inhibited  if  SSWll  is  on  (i.e.,  the  background  timer). 
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CARD  MOniFICATION 

To  alter  locations  from  the  card  reader,  ; ;h'8  should  1)0  set.  Cards  '-iJ  1 
be  read  and  printed  on  the  line  printei . No  cards  will  be  processtd  untia 
the  first  keyword  card  is  read.  This  card  has  the  word: 

ACVTN  LT 

on  it,  starting  in  column  one.  Cards  w:]J  continue  to  be  read,  regardless 
of  the  state  of  SSW8,  until  the  second  keyword  card  is  read.  This  card 
has  the  word: 

EKDACVIN 

on  it,  starting  in  column  one.  All  the  cards  from  ACVINIT  to  EKDACVIN 
will  be  processed  for  valid  memory  alterations. 

There  are  two  ways  to  change  memory  locations  using  the  card  reader.  The 
first  is  by  using  the  keyword  IfliMORY.  The  second  is  by  using  the  keyv/ord 
PARAM.  Every  word  in  memory  can  be  changed  with  tlie  MEMORY  keyv.’ord.  The 
PARAM  card  uses  pre-assigned  parameter  names  to  access  mem.ory.  Only  a 
small  number  of  often  used  locations  can  be  accessed.  PARAM  cards  do  not 
have  to  be  changed  after  reassemblies. 

The  formats  of  the  MEMORY  card  and  the  PARAM  card  are: 

MEMORY  ADDRESS,  DATA,  DATA,... 

PARAM  PARAMETER,  DATA,  DATA,... 

In  both  cases  the  kejn^ord  must  start  in  column  one.  After  the  keyword  the 
spacing  is  free  format.  ADDRESS  is  a nup;eric  value  and  is  assumed  to  be  an 
absolute  memory  address.  PARAMETER  is  an  alphanumeric  word  internally 
associated  with  a specific  memory  address.  The  address,  whether  absolute 
or  narametric,  will  be  the  starting  address  in  which  to  store  the  data 
that  follows  the  address.  The  DATA  values  are  stored  in  consecutive  memory 
locations,  starting  with  the  first  address.  The  format  of  the  DATA  values 
includes  a value  refered  to  as  a multiple.  It  is  used  to  repeatedl''  store 
a single  DATA  value  in  a number  of  consecutive  locations.  That  is,  if  n 
DATA  value  has  a multiple  of  10,  10  values  of  DATA  will  be  stored  into 
locations  ADDPJESS  to  ADDRESS  + 9 (or  PARAMETER  to  FAPd^^ETER  + 9) . 

ADDRESS  has  the  same  format  as  the  DATA  words.  However,  an  multiple,  that 
is  provided  with  the  ADDRESS  value  is  ignored.  PARAMETER'S’ are  preassigned 
alphanumeric  words.  A list  of  those  v/ords  that  can  be  used  as  PAPjYlil.TER's 
is  given  in  Parameter  Values  for  Initialization.  The  format  of  the  DATA  value 
is  given  below. 

If  the  ilEMORY  card  or  the  PARAM  card  ends  v/itli  a comma,  data  can  be  con- 
tinued onto  the  next  card.  Data  on  a continuation  card  can  start  on  anv 
column.  Continuation  cards  can  also  have  continuation  cards.  A card  that 
ends  with  a comma  does  not  necessarily  have  to  have  a continuation  card. 

A valid  keyv/ord  in  column  one  will  terminate  a cor.rinuation  string.  Onlv 
DATA  x^ords  can  be  on  continuation  cards.  The  ADDRESS  or  PARAMETER  word 
followed  by  a comma  must  be  on  the  first  card. 
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If  any  part  of  a card  cannot  be  tran.-jiated,  an  error  ine.-ij'.ajje  \.'i  Lj  be  ir.rt-.-d 

on  the  line  printer  underneath  the  card't>  ; rir.ted  imaee.  The  t rate  ■ • .1' 

have  to  localise  the  error  on  the  card.  I'n  error  dot  ^-  "ot  caiv  t;  e 

tire  card,  however.  If  a valid  addresr  aisa  a ■•.•.idif  icat  <011  va]  i..,  afi  a* 

determined  by  the  card  processor  before  tb.c  error  l < found,  that 
modification  will  be  performed.  Data  or.  the  card  Iron  tht  crroi  to  the 
end  of  the  card  will  be  ignored  and  any  continuation  cards  wilJ  be  iynir.d. 


DATA  words  have  one  of  the  following  forr'ats: 


MMM  D’iDDDDDD' 

MMM  D't.DDDDDD 
MMM  F'D.DDDDDEdrDD’ 
MMM  X'KHHHKH’ 


decimal  integer 
decimal  fraction 
floating  point 
hexadecimal 


In  all  cases  MMII  is  the  multiple  value  mentioned  previously.  Tliest*  cliuracters 
must  be  decimal  digits  and  can  be  any  number  of  characters  as  long  as  thr- 
word  does  not  overflow  a 32  bit  w’ord  magnitude.  The  multiple  cannot  have 
a sign  and  is  always  positive.  There  must  be  at  least  one  blank  after  the 
multiple. 

The  decimal  integer  and  fraction  formats  (D)  renresent  cither  a 32  bit  deci- 
mal integer  scaled  B31  or  a 32  bit  decimal  fraction  scaled  BO.  The  iol- 
lowing  rules  apply. 

1.  Data  must  start  with  D' 

2.  The  sigh  (+  or  -)  is  optional.  Plus  is  assum.ed. 

3.  The  first  character  can  be  a decimal  character  (0-9)  followed  b- 
a single  quote  or  a decimal  point  (.). 

4.  The  remaining  characters  must  be  declm.al  characters  (0-9)  follcvod 
by  a single  quote  ('). 

The  decimal  point  cannot  be  preceded  by  any  decimal  characters.  Too  man” 
characters  in  an  integer  v;ill  cause  an  arithmetic  overflow.  Extra  diglL.s 
in  a fraction  will  be  rounded. 

The  floating  point  format  (F)  represents  a Sigma  7 single  floating  point 
number.  These  rules  apply. 

1.  Data  must  start  with  F*. 

2.  The  sign  of  the  mantissa  (•!•  or  -)  is  optional.  I’lus  is  assumed. 

3.  The  first  character  must  be  a decimal  digit  (0-9) . 

4.  The  second  character  can  be  a decimal  point  or  a decimal  digit. 

If  the  decimal  point  is  left  out,  it  will  be  assumed  to  bo  after 
the  first  digit  (3  above). 

5.  The  mantissa  can  be  continued  to  any  number  of  digits.  Jt  will 
be  rounded. 
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6.  The  exponent  (E±Dl))  is  a I'ov.c  r o£  10.  it  i-  ontional  .”10  ji  Ifit 
out,  00  is  assumed.  Its  sien  i,  optional.  It  c.iP  ’iave  pjtlicr 
one  or  two  decimal  dipits. 

7.  The  floating  format  must  end  witii  a single  quote  ('). 

The  hexadecimal  format  (X)  repre.sents  a 32  bit  hexadecim.il  vaiiu  . ’’he  ruJe;- 
for  the  hex  format  are. 

1.  Data  must  start  with  X'. 

2.  Any  number  of  hex  characters  can  be  u.sed. 

3.  Data  must  end  v;ith  a single  quote  ('). 

The  data  is  stored  right  justified  witii  left  ch.aracters  zeroed.  :i  the 
data  word  is  more  than  eight  characters  Jong,  the  right  most  eight  uiuuac- 
ters  will  be  used.  Hexadecimal  data  words  have  no  sign. 

In  all  DATA  formats  there  can  be  no  embedded  blanks.  DATA  values  cannot 
be  continued  from  one  card  to  another.  Some  examples  of  accentable  and 
non-acceptable  DATA  words  are  given  below. 
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■ 

1 

i 

Acceptable  Data  f 

1 ’ 

Decimal  Integer 

i 

iVleta*Sym  hnl  l''f|uiY2.1ont 

1 

D 

o 

1 

D'2'  2 

m 

D*+48'  48 

H 

K 

D'-3976'  -3976 

B 

Decimal  Fraction 

K 

DATA  Meta-Symbol  Equivalent 

1 

FX'OBO' 

FX\5B0* 

FX'.623B0'  1 

i 

FX'-.9289B0'  ; 

1 

D'-.  62314383678901543'  FX'-.  623143837B0'  i 

1 

Floating  Point  ; 

f: 

Meta-Symbol  Equivalent  = 

M. 

FS'O'  |j 

1 

FS»5'  jj 

E 

F'+5'  j-g,-,  1 

K 

F»-6.2'  FS»-6.2'  1 

1 

F'-628E2'  FS'-628'orFS'-6.28E2'  | 

I 

F'-7.45E-10'  FS'-7.45E-10' 

1 

F'4. 666666666666'  FS'4.  66636667' 

1 

F'-^1800E3'  FS'1800' or  FS'l.  800E3' 

m 1 

F'2. 78E+28'  FS'2.78Er28'  I 

m 


1 f 

X"  X'O'  ( J ll 

1 1 

X'ABC  II 
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I 


i 

r 


f 


DATA 

X'ABCDE123456789' 

X'ABCDDCBA' 

Unacceptable  Data 
Decimal  Fraction 


Meta-S3Tnbol  Equivalent 

X'23456789' 

X'ABCDDCBA' 


DATA 
D'6.  ' 

D'.  999999999999' 
Floating  Point 


Reason  For  Error 

decimal  point  in  wrong  place 
data  overflow  after  rounding 


DATA 


Reason  For  Error 


F'65. 7' 

F'7.  82E' 

F'.35E2' 

F'5.32E-106' 

F'9.999999999999E-1' 


decimal  point  in  wrong  place 
no  exponent  characters 
decimal  point  in  wrong  place 
too  many  chare  ters  in  exponent 
data  overflow  alter  rounding 


None  of  the  examples  above  were  given  multiples  but  there  is  no  reason  why 
they  could  not  have  had  them.  The  following  lines  are  examples  of  MEMORY 
and  PARAM  cards  having  DATA  values  with  multiples  attached. 


MEMORY  X'2A01',F'2.0',  10  D'48' 
MEMORY  X'IBCD',28  X'1234*,D'14', 
F'6.E20',  4284  D'-64' 

PARAM  U0,F'450E2' 

PARAM  WIND,F'4.5E1' 

PARAM  MATRIX0.9  F'0» 

PARAM  MATRIX1,F'1',3  F'O',  F'l', 
3 F'O',  F'l' 
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TYPEWRITER  MODIFICATION  AND  DISPLAY 

The  typewriter  can  be  used  to  displav  aii.i  alter  i:-er.orv  locatiuns  and  Li 
set  or  change  the  date  and/or  time.  Tvnc-ins  can  be  done  anv  time  during 
the  ACV  simulation  as  long  as  the  typewriter  light  is  on.  Thir.  light  Js 
alternately  put  on  and  off  by  activatinf  the  console  interrupt. 

The  time  and  the  date  input  were  described  above.  Moivory  locations  can  bo 
altered  by  a MEMORY  or  a PARAli  type-in.  The  format  of  these  tyne-ins  are 
almost  the  same  as  the  card  l^EMORY  and  ‘'ARAM  inputs.  The  only  difference 
is  that  the  keywords  MEMORY  and  PARAM  do  not  have  to  start  cn  column  1 . 
Continuation  lines  can  be  used  with  the  i’APAl!  and  >’F.MOr.Y  tyro-ins. 

The  contents  of  memory  can  be  typed  onto  t!\e  typewriter  by  the  t\pe  keyword 
DUMP.  The  formats  of  the  DUMP  type-ins  are: 


DUMP  D DATA 
DUMP  X DATA 
DUMP  E DATA 


DUMP  F DATA 


decimal  type-out 

hexadecimal  type-out 

mantissa-exponent  type-outs  of 
single  floating  point  values 

integer  fraction  type-out  of 
single  floating  point  values 


The  spacing  of  the  type-ins  are  free  format.  As  many  spaces  as  desired 
(including  none)  can  be  typed  between  DUMP  and  the  forr.:at  characters  (D,X, 

E or  F)  or  between  the  format  character  and  the  DATA  value.  DATA  represents 
a multiple-value  pair  as  described  on  the  MEMORY  and  PARAM  DAT.A  rcords.  In 
this  case  the  multiple,  if  present,  is  the  number  of  word  locations  to  be 
typed  out.  The  value  of  the  DATA  word  is  the  starting  address  of  rhe  con- 
secutive locations  to  be  typed.  DUMP  type-ins  cannot  be  continued  cr.  a 
new  line.  The  DUMP  type-in  will  type  out  pairs  of  values  vjhich  renresen.t 
the  address  and  the  contents  of  the  address  of  the  reauested  locations. 





NAVTRAEOUirCKN  75-C-0057-1 


The  four  different  output  formats  are  as  follows: 


D The  decimal  equivalent  of  the  word  address  contents.  Leading 

zeros  are  suppressed  and  the  value  is  left  adjusted.  The  format 
is  SDDDDD.  S is  the  sign  + or -.  D's  are  decimal  digits  (0-9). 

X The  hexadecimal  equivalent  of  the  address  contents.  Leading 

zeros  are  typed.  The  format  is  HHHHHHHH.  H's  are  hex  digits. 

Eight  characters  are  always  typed. 

E The  mantissa  - exponent  equivalent  of  the  word  address  contents. 

The  format  is  SD.  DDDDDESDD.  S is  the  sign,  or  (»f  the  mantissa 
or  the  exponent.  D is  a decimal  character  (0-9).  E is  the  character 
'E'  and  represents  a power  of  10. 

F The  integer  fraction  equivalent  of  the  word  address  contents. 

This  is  for  a floating  point  value.  The  format  is  SDDDDDDD.  S 
is  the  sign,  + or  -.  D is  a decimal  character  (0-9)  or  a decimal 
point.  The  position  of  the  decimal  point  will  vary  depending  on  the 
magnitude  of  the  value.  Leading  zeros  are  suppressed  (except  the 
value  zero)  and  the  number  is  left  adjusted.  If  the  value  is  too  small, 
zero  will  be  typed.  If  it  is  too  large,  i 999999999  will  be  typed. 

For  all  type-ins  the  following  rules  apply.  This  includes  TIME,  DATE, 

MEMORY,  PARAM,  and  DUMP  type-ins. 

1.  Type  lines  should  end  with  a new  line  character  (NL). 

2.  A typed  line  will  be  terminated  after  85  characters  and  will  be 
rejected  by  the  1/  O processoro 

3.  To  cancel  a line,  type  in  an  end  of  message  character  (EOM) 

4.  Tlie  horizontal  tab  character  (HT  or  TAB)  acts  like  a NL  character. 
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Incorrect  characters  (except  for  NL,  EOM,  or  HT  characters)  can 
be  corrected  by  typing  a ^ (cent  mark;  for  each  character  that  was 
previously  typed  and  is  to  be  deleted  (e.g.,  DUNPC^MP  is  equi- 
valent to  DUMP). 
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ACV  SYSTEM  PRINTOUT  DESCRIPTION 
Heading  Line: 

ACV-LC,  Time  of  Day,  Date,  Current  System  Status 

The  current  system  status  consists  of  a ~ character  code  with  values: 


INIT 

HOLD 

0? 

OPGO 

CRSH 


Reset  Mode- 
Freeze  Mode 

Operate  Mode  without  Pilot  (or  < 10  sec) 
Operate  Mode  (Simulation  Operable) 

ACV-LC  Hit  Mother  Ship 


In  body  axis  coordinates,  the  forces  and  moments  are  shown  that  act  on  the 
vehicle  broken  into  the  various  components.  The  total  forces  and  moEents 
•are  shov/n  and  through  Newton’s  Second  Law,  the  body  axis  acceleratio:is  are 
given.  Initial  velocities  are  printed  and  by  integration  the  change  iii 
velocity  and  the  current  velocities  are  printed  in  body  axis  coordinates. 

Perforndng  a coordinate  transformation  yields  velocities  in  earth-fixed 
coordinates.  Integration  yields  change  in  position  and  summing  with  initial 
positions  yields  the  current  position  of  the  vehicle.  Angular  position  is 
dispiaved  in  both  degrees  and  radians. 

The  effectors  consist  of  the  various  physical  devices  that  directly  cause 
the  vehicle  to  maneuver,  pilot  coiranaiids  to  these  devices  where  applicable, 
and  other  conditions  which  are  measurable,  and  environmental  effects.  Refer 
to  Figure  55  for  a sample  printout. 


if 

I ’1 
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The  list  is  as  follows: 

STBD  and  PORT  propeller  speeds  {NPS,  NPP)  (RPM) 

STBD  and  PORT  propeller  pitch  (BETAPS,  BETAPP)  (DEG) 

Rudder  Angle  (PSIR)  (DEG) 

Thrust  Nozzle  Angle  (PSIN)  (DEG) 

STBD  and  PORT  propeller  thrust  (TPROPS,  TPROPP)(LBS) 

STBD  and  PORT  propeller  pitch  commands  (BETAPSC,  BETAPP C) (DEG) 


Rudder  Command  (PSIRC)  (DEG) 

Nozzle  Command  (PSINCA)  (DEG) 

Cushion  Compartment  Pressures  (PCUSH  1,PCUSH2,PCUSH  3,PCUSH4) 
(LBS/FT^) 

Manifold  Pressures  (PMANS,PMANP)  (LBS/FT^) 

Air  Flow  through  STBD  and  PORT  fans  (QFANS,  QFANP) 

(FT^/SEC) 

Air  Flow  through  STBD  and  PORT  Nozzles  (QNOZS.QNOZP) 

(FT^/SEC) 

STBD  and  PORT  fan  speeds  (NFANS.NFANP)  (RPM) 

Wind  Velocity  (WIND)  (FT /SEC) 

Wind  Direction  - Set  (PSI WIND) (DEG) 

Cushion  Compartment  Volume  (VCUSHl.V  CUSH2  , V CUSH  3,  V C USH4) 
(FT^) 

Air  Flow  due  to  wave  pumping  into  Cushion  Compartments 
{QPUMP1,QPUMP2,  QPUMP3,QPUMP4)  (FT^/SEC) 
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Skirt  Gap  Area  under  each  Compartment  (AREA  1,  AREA  2,  AREA  3^ 
AREA  4)  (rf ) 

STBD  and  PORT  power  turbine  speed  (NS3,NSP)  (RPM) 

STBD  and  PORT  gas  turbine  speeds  (NTS  1,  NTS  2,  NTS  3,  NTP  1, 
NTP  2,  NPT  3)  (RPM) 

STBD  and  PORT  power  turbine  commands  (NSCS,  NSCP)  (RPM) 

STBD  and  PORT  gas  turbine  commands  (NTCS  1,  NTCS  2,  NTCS  3, 
NTCP  1,  NTCP  2,  NTCP  3)  (RPM) 

Offshore  wave  height  (HEIGHT)  (FT) 

Offshore  wave  period  (PERIOD)  (SEC) 

Ocean  bottom  slope  (SLOPE ) 

Offshore -shoaling  boundary  (XSHOAL)  (FT) 

Offshore  wave  length  (LAMBDA)  (FT) 

Wave  encounter  frequency  (OMEGAE)  (RAD/SEC) 

System  parameters  are  displayed  to  show  specific  status. 

Execution  Time  (TIMEX)  (SEC) 

Divell  Time  (/SEC) 

Time  Overflow  Counter  (EXSELFIN)  (COUNTS) 


227 


I 


NAVTRAEQUIPCEN  7-)-C-00  57-l 


«o«oco 

»n  •Ck*o»r 

• o*«***  •^v*^***^ 

-I  ••  *«r.  o Of?*^  *«*0‘*‘*' 

*-f*  <3  * aO^v^OO 

«***«Ofu  «•«.>•■•  •««<jni  MMiu  • « • * • « 


;w 

•«i*  M ■» 

<ct^  « n 
van  « V O 
^ n < & tO  M 
« >•  OOOM 
•-«*  41  MM#!.' 


• ^ • oo 

% «a  • ^ • wt  tp 

• • • ^ • *r.  c» 

O M ^ • Af  fW 
>a  — 

t «0  • m ^ J' 

O ■«  • et  9 O O 

r»  ••  • O • o o 


4ia>  4im9n  cf^ « 

•*t  n ^••a  *b«41 

^«>  «'<«o 

• • (Mf^PkAt  0A4S 

fv  A Ok 

• Ci  lA  • A»p»a» 

^ rw  «A«o>0k  -xrto* 

I pk  o A*C«* 

' MAf  rkM«^ 

• t •««#  «•« 


^ PI « ^ «• 

9rift-^9K  «l 

«*  a 

AJ  CU  |H 

• <UAI^  01 

e\  A*  Sin  c 

••O0%O>  M 

• ^o»  ^ «• 


o»  • O'  -»  « 

«rt  * "» Pk 

p*i  O • — « <1 

'O  fw  P> 

« • , « « • 


• 9‘  * 0tf>0» 

• fs  • o«*n 

• «.  • C nC  4» 

• ««  « O ^ ^ 

• ^ 9k 

• ^ « OOO 

• A<  • U • • 

• • • o ♦ r 


« W W fitrfk****  Pk.^lAe 

X ■•  IkAsA*}***!:.:.! 

*-  ^ •-  •*  asji-^-a^-^^-aaAa. 

u «4  u « w u 14 14  u « ^ 7 :r 


Ok  • i'A*  -f  « 
4r  • O >«  «•  « 

«».  « rf*  lA 

* -n  * 
pk  • r*  ^ p»  -• 

•>  • ^ rfi  p» 


X*  A OUP  A< 

*W  A ^ J ISd 

O A (Tk  7 «P 

X • #k.  ^ 

t>i  • O C 


>A  AOOfW^  Aia 

‘A  Owf«  AOO*A 

kA  C^0k«O  »«*0  •« 

« A kO^A  # »• 

k • » mm  UI^^AffU  pltf 

>|/»  fSP»*  tfAfAi**# 

> ^ #••  OAA*«  «l 


> O Af 

• Pk  P 
k • U«<k  ^ 


A Pk  • 

A — A 9 V « 
A «»  A 9fn  -» 
A P.  A 9<V  A| 
A P A 9 • * 
A ^ A ,9AI  « 
AAA  909 
AO*  OOO 
A •■  A « »««  ar 


iPkPk-SM  OI-*P 

tiomoiok  aaop 

^cvpkOk  o 

^ «U  PI  A*  Ok  A 

» « 4>  Ak  « • lO  .« 

• A>  « * * M A < 

» A A O*” 

• p»«Xmaa  JiCkA 

» XKWaa^  M-*  • 

> A A A * AAA 


•VPk  Ji 

O-A 

• A A 9 


A A 

r o-W-Atfl-s  C O AC«»»A 

tPk  W>  A-WI  0*U«A«Pk  O AUTIOfA 

I eiAi  Aj3k#AA 

ttflA^fkOlAA*  A A « AOA*OA 


AOtAX^O  A o Oin#  0»  AA  !■»  « IMMV  ^ o A o O—  A ?• 

A«A9I««  A01P-»*»«'^0>-«>»-APkO  AUO<UAn 

AAOAJkOlAI  A tXUn»C^  *>  A C,  A OAe  ^ A 4!  O O * ^ 

• 4rA*«*  AMfW«r>APiAk«AI*A*A  ACrVoAA 

• «>kAA.#Bk««  AAAkVr‘AAXf«<>AAA*C  AAkOf^A** 

ArA'^A'A  AA«AAk^O«»U^^9  *900aa 

A aAOv-^I^  AAACI^7«*^VV'^^0  A09A<A'X 

ACAASkAlAI  A A(  0«*«*  0 9k  «Af*l  « « « « A C A ** 

AAk.AfV«A^  A-A«— ^^  «•«««*  « * A ^ 


A ^ < a • 
A W X «k  A 
> OO  OO  A 


^Al  Al  4^Ai.^a*U.AAkkMA*A 

O U A*  <0  a.  atx  aO^^^VIa 

^A*  sc'^^A.**  —V*  A 

< ww  kkt^AO  v»*p*a,c^=;a 

^ i.  J*n  ^ yt  4,«ia«^«*47a 

«•••  e^7^A>W«A^kM  aZw^XIA  a 

— &a  •^>t.X  A 

07  aw«-;ik4.i^7AVT  o^wiA 
VO  k«/kfa.;.74^«  « ir.  A 

O W«WA 

lAw  s>cLtfaA'^Poo»  ^ mi  9 .J 
mf  » O « i •--•«>-  A Wj«J  J 

• > V A JL  ^ s o ^pa^k— •AAAa 

>0  7 •“  •••-  A-  .«.  ^ 7 Al  A •-X  A A-WI^ 

V*  O ««  A*  %«  Xf  C' W Mi— ^ 7 A W =«»C  W o 


70^  ^ A • 
W ^ A*  A U A 

«*  * w**o^ 


A 7<PaO  a • « • 

A , A&«a»*a  m * • 

iM  u * V • • • « 

^ X7  ;cw<iiA  « • A • 

O A>  A>  ' M|X  WOl  • • • A 

(9  « ^ tfl^SAA  « • V A 

*A  •Am  7(090  A A A A 

iO  XX  *a9l^  A m A A 

«0  7 77  AS/I«7  • A A A 

V A A.  ^ -ao  a — ki*  • • • • 

AA  < V*  A a777^AAA  A o 

* •m  M— AA^  a y»  AO  A ^ 

X V A-  ^V*  AOaJa  aaOA 

M j AAAin  WAAca77*w*  wxx 

U •>  MMMO.  . a<V^^«7  A A au^a  ^ 

tf  A AaO*wA7A^«.XO 

o aa  a a Aca  <saOIo  aa  w«ax  w 

r*  X aaaamA  0»A-~  X a*  7 

(X  7 MM^tr  UOlWWw 

V AAAa  «s;«x.  ^o«  X X 

0~*AA  AAa^  ^AXk9A>^^A^AA^AA 

;m  V»C  AAAjXAfari  tf«OAA  a7X7- 

i«A.(dV«AA,i  A X (/is:7ktfAAf^>aa.  aa(a  mic 

A.  AAA.  ^AAAA»  XaA 

«^a-aa0^4  a a £«7  (O  X O a«  ^ A*Otf  **  a 
• -aaaaax  9 V lAk^A^OMaf^x  ^ ^ 

iw— ^«*wwo  wxc.aa«^a>aa  ^ X a 

— « -A  «A»A  « X X A.  A C A.  AMX  M ^ 

a ^ Mkinx  «:7X4*PAAVvX»0  * 

X A«AA  #*4XOX  (AaA^^AAA^XV-^OXSf 

a>i^A7a«a*^7  a«X  a*  W 7 7 O X (A  is  7 07 


* X P W X A 
X X A ^ A 

— %>  X Ai  ^ A. 
X X X X A>  4. 


NAVTRAEQUIPCiiN  75-C-0057-  1 


PARAMETER  VALUES  FOR  INTtiaLIZATICN 


uo 

Velocity  (linear)  in  direction  X (Ft/sec) 

b 

VO 

Velocity  (linear)  in  direction  Y (Ft/sec) 

b ■ 

wo 

Velocity  (linear)  in  direction  Z (Ft /sec) 

b 

Po 

Velocity  (angular)  about  axis  X (rad/sec) 

b 

QO 

Velocity  (angular)  about  axis  Y (rad /sec) 

b 

RO 

Velocity  (angular)  about  axis  Z (rad/sec) 

b 

XI  0 

Inertial  position  (North)  of  pilot  cabin  (Ft) 

yi  0 

Inertial  position  (East)  of  pilot  cabin  (Ft) 

ZIO 

Inertial  position  (Vertical)  of  pilot  cabin  (Ft) 

PHIO 

Euler  Angle  Roll  (Deg) 

THETA  0 

Euler  Angle  Pitch  (Deg) 

PSIO 

Euler  Angle  Heading  (Deg) 

WIND 

Wind  Velocity  (Ft/sec) 

PSIWIND 

Set  of  Wind  (Deg) 

DEETAT 

Sampling  Period  (Sec) 

PGusmo 

Pressure  within  cusMon  compartment  1 (Lbs/Ft“) 

229 


NAVTRAEQULPCKN  75-C-0057-1 


PCUSH20  Pressure  within  cushion  compartment  2 (Lbs /Ft  ) 

2 

PCUSH30  Pressure  within  cushion  compartment  3 (Lbs /Ft  ) 

2 

PCUSH40  Pressure  within  cushion  compartment  4 (Lbs /Ft  ) 

2 

PMANSO  Pressure  within  ST BD  manifold  (Lbs/Ft  ) 


PMANPO 


Pressure  within  PORT  manifold  (Lbs/Ft'^) 


NTSIO 


NSS0 


Gas  Turbine  Speed  - STBD  1 (RPM) 
Power  Shaft  Speed  - STBD  1 (RPM) 


VGWFLAG 


Set  = 1 to  run  VGW 


BETAPS0  Propeller  Pitch  - STBD  (Deg) 

BETAPP0  Propeller  Pitch  - PORT  (Deg) 


PSIN0 


Thrust  Nozzle  Angle  (Deg) 


PERIOD 


Offshore  Wave  Period  (Sec) 


SLOPE 


HEIGHT 


Ocean  Bottom  Slope 
Offshore  Wave  Height  (Ft) 


CSKIRT 


Skirt  Stiffness  Coefficient 
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SUMMARY  OF  ACV  OPEFATINO  PROCEDURES 
System  Sense  Switches  (SSW) 

1.  RESET  (Simulation  mode) 

2.  EREEZE  (Simulation  mode) 

3.  OPERATE  (Simulation  mode) 

8.  Read  Cards 

0.  Print  State  Parameters 

11.  Time  Background  (Units  of  12  U secs) 

15.  Print  Datapool  Contents  (FS) 

Console  Sense  Switches  (CSW) 

1.  Inhibit  ACV  Simulation  (Do  not  run  program  ACVSIM) 

2.  Inhibit  ACV  Cockpit  Command  Inputs 

3.  Inhibit  CRT  Driver  (DLOOP) 


Console  Interrupt 

Alternately  allow  or  disallow  type-ins 
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SUMMARY  OF  ACV  OPERATING  PROCEDURES  (CONTINUED) 


Traps 


Error  Locations 

C4ERR#  - right-most  hex  character  equal  to  trap  error  number 
C4ERRLOC  - PSWl  of  trap  PSD 

C4SRRMEM  - memory  module  fault  bits  (24-31) 

Trap  Error  Number 

0 Memory  Parity 

1 Non-allowed  Operation 

2 Unimplemented  histruction 

3 Push  Down  Stack  Fault 

4 Fixed  Point  Arithmetic  Overflow 

5 Floating  Point  Fault 

6 Decimal  Arithmetic  Fault 


7 Watch  Dog  Timer 

Trap  Restarts 

1.  CSWl  Off- IDLE -RUN 

2.  CSWl  On-IDLE-RUN 

3.  IDLE-CPU’RESET-RUN 


(pick  up  from  error) 
(reinitialize  - only  after  trap) 
(reinitialize  - anytime) 
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SUMMARY  OF  ACV  OPERATING  PROCEDURES  (CONTIMUED) 


Card  Inputs 


ACVINIT  Start  card  processing 

MEMORY  ADDRESS, DATA, DATA, . . . Absolute  memory  alteration 

PARAM  PARAMETER,  DATA,  DATA, . . . Parametric  memory  alteration 
ENDACVIN  End  card  processing 

(SSW8  can  be  put  off  after  ACVINIT  is  read.  Card  reading 
will  continue ‘.until  ENDACVIN  is  read) 


Type  Inputs 


TIME 

HH:MM  Set  time  (HH  = 

0-23,  MM  = 0-59) 

DATE 

MON  DD.  'YY  Set  date  (MON 

= JAN,  FEB,  etc. . . , one  space 

DD  = 01-31,  'YY  = '74,  '75, '76) 

MEMORY 

ADDRESS,  DATA,  DATA, . . . 

Absolute  memory  alteration 

PARAM 

PARAMETER,  DATA,  DATA , . . . 

Parametric  memory  alteration 

DUMP 

D 

DATA 

Dump  integer  in  decimal 

DUMP 

X 

DATA 

Dump  integer  in  hex 

DUMP 

F 

DATA 

Dump  floating  point  word  as 

integer-fraction 

DUMP 

E 

DATA 

Dump  floating  point  word  as 

mantissa-exponent 

Type  Control  Keys 


NL 

.HT 

EOM 

85  characters 
n(? 


End  and  process  message 
End  and  process  message 
End  and  ignore  message 
End  and  ignore  message 
Ignore  last  n characters 


3| 
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SUMMARY  OF  ACV  OPERATING  PROCEDURES  (CONTINUED) 


Input  DATA  Format 


MMM  D'+DDDDDD' 

MMM  D>± . DDDDD* 

MMM  X'HHHHHH’ 

MMM  F'±D.  DDDDDE+DD' 

MMM  - represents  a 

or  dumped. 


Decimal  integer  B31 
Decimal  fraction  BO 
Hexadecimal  integer 
Floating  point  single  value 

multiple  for  consecutive  locations  stored 
It  is  assumed  to  be  1 if  not  given. 


ADDRESS  Format 

The  same  as  DATA  but  the  value  is  considered  a memory  address  and 
the  multiple  (MMM)  is  ignored,  if  given. 


I i 


« 


ii 

i 


1 1 
i I 
i 


e 

I 
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SECTION  vr 
DERIVATIONS 

Thuse  equations  {sections  of  the  model)  that  changed  to  a degree  that  might 
not  be  recognizable  from  those  shown  in  the  report  of  the  math  model,  the 
derivations  are  shown  in  this  section.  The  only  model  that  does  not  either 
coincide  v/ith  the  math  model  report  or  is  not  shown  here  is  the  compartment 
pressure  model.  The  air  pressures  are  solved  in  a similar  manner  (that  is 
shown),  however,  tlie  derivation  was  performed  by  Damon  Cummings, 

CHANGE  IN  INERTIAL  POSITION 

This  dei'ivation  was  performed  on  the  second  order  predictor -corrector 
integration  scheme  to  integrate  velocity  to  a change  in  inertial  X and  Y 
position  for  every  sampling  period.  This  was  done  to  improve  accuracy  on 
a delta  X calculation,  subtracting  one  large  value  from  another  large  value. 
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APC  =4/5  (C^  - P^) 

=P^  +APC 
= P^  + 4/5  (C^-P^) 

= Ct  + 1/5  (P^-Cp 

° ^t-  At  «t  At*  ■"  k-2  At  2 ^t- 

-X  ..-^(X  --  '1 


t- At 


V % f \\, 


At 


*^t  - ^t-  At  = + 1/2  At  1/5 

+ 2/5  At  ^ - 1/10  At  X,  - 1/lOAtX. 

t " /\X 

A}q  = 4/5AtX,_^^ + 2/5AtX^  - 1/5AX^_^^ 

= 1/5  [2At(2i^.^^+i^,-AX..^l 
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SEAWAY 


The  wind  generated  waves  have  complex  equations  in  the  shoaling  region 
for  the  first  order  amplitude  and  the  second  order  Stokes'  amplitude. 

Both  equations  taken  from  the  math  model  contained  hyperbolic  fimctions. 

The  hyperbolic  identifies  and  the  wave  number  equation  were  substituted 
appropriately  to  obtain  simpler  equations  thft  would  decrease  computation  time. 

Reduction  of  the  complex  shoaling  equatioii  is  the  deep  ocean  boundary  is 
show'nto  vastly  simplify  the  offshore  wave  ec.-.ation. 
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2 Cosh  “ 

2 K D + sinh  2 K D 

XX  X >: 


1 - tarih'^  K..  D. 


K D (1-tanh  K . D„)  tanh  K D 


D -K^ 

X X ox  o 


2 K D„,  1/2 


tanh  Kj. 

(1-tanh^  K.  D ) 

X x' 

[3-tanh^  K.  D 


1-tavh  ^ K D 
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SECTION  VII 

HULL  BOTTOM  PLATIN’C  - 35  POINTS 

This  section  shows  the  equations  in  present  form  of  the  original  set  from  the 
"Mathematical  Model  of  an  Air  Cushioned  Vehicle",  These  equations  show  the 
mathematics  of  the  forces  and  moments  due  to  the  wave  heights,  cushion 
compartment  voliunes,  and  skirt  escape  areas  for  the  original  35  points 
under  the  hull  planform.  The  number  of  points  was  decreased  in  an  attempt 
to  get  the  VGW  to  run  real-time.  These  equations  are  presented  here  for 
completeness  only. 

Three  subroutines  and  the  DATAPOOL  changed  as  a result  of  changing 
from  35  points  to  25  points.  Other  equations  and  subroutines  depend  upon  these 
htiU  point  poeitions  and  the  numbering  scheme,  however,  that  is  simply  handled 
by  changing  the  number  of  points  and  their  locations  within  the  DATAPOOL. 

See  Figure  56  and  Table  12  for  hull  point  definitions. 

These  equations  necessitated  a change  as  special  numerical  techniques  were 
developed  for  solution  purposes. 
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Subroutiiic:  Si-IA 


P 

(n+n 

+n  - n *i 

3 

a’Si  1 

1 2 

3 6 

+p 

( n+n 

+n  -n  - 

n - n 

CUSII2  6 7 

8 21 

27  26 

+P 

( n + 

n “ 

n - n 

CUSH3  8 

16  ' 15 

19  20 

■'■P  (n  +n  +n 
CUSH4  1 11  12 


" 

n - n -r.  ) 

8 16  IS 


-4^ 

4 


P (n  +n  +0  +n  -n  -ri  -n  -n  ) 

CUSH  1 3 4 5 6 1 10  9 8 


( n + n +n  +n  -r,  -n  -n  -n  ) 
aJSH2  6 24  25  26  S 23  22  21 


+P  C n + n + n +n  - r,  -n  -n  -n 
CUSH5  21  22  23  8 15  17  18  ]( 


(n  + n + n+n-n-n  -n  -n' 
4 8 9 10  1 12  15  14  ’ 15' 


CUSHl 


(20)  (8)  ( n n + n -n  -n  - n ) 

5 6 7 8 12  3 


+(40) CIO)  (n  +n+n+n-n-n-n  -n  ) 
4 654389  10  l' 


*P  (40) (50)  f n + n + . . 

CUSH2  [-4—  26  25  •24  % '21" 

( n + n + n -n  - n - n )] 

4 21  27  26  8 7 6 


■^P  (20)C28)  fn  +n  +n  „ 

CUSHsL'-I 19  20  ■ 21  15  16  8 


+ (401(50)  (t,  + n +n  + -;-r  - n - n 

^ 21  22  23  S 19  IS  17 


+P 


CUSH4 

+ f iin'i  fi  n"!  f 


(20(28)  ( n + n + n - _ 

.3  8 15  Ih  I 11  "1- 


+ n 


-p 
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H 

CUSH  1 


H 

CUSH  2 


H 

CUSH  3 


H 

CUSH  4 


V = 

CUSH  n 

n = 1-— 4 

h , 
wk 

^"CUSH 
^CUSH  n 

’^CUSH  n 
k = 1---35 


Subroutine:  VOLCl'SIl 


.020833  (h  +h  +h  +h  ) + . 0625  (h  'h  Hi  +h  ) 
wl  w3  w6  w8  w4  w5  \v9  wlO 

+ . 08333  (h  +h  ) + . 25  (h  + h ) 
w2  w7  w30  w31 

= .020833  (h  +h  +h  +h  ) i . 0625  (h  +h  +h  ^ +h 

w6  w26  w22  w23  w24  w25 

+ .08333  (h  +h  ) + . 25  (h  +h 

w7  w27  vv28  w29 

. 020833  (h  +h  +h  +h  ) +.  0625  (h  +h  +h  +h  ) 
w8  w2l  wl9  wl5  wi7  wl8  w22  wr.3 


+ . 08333  (h  +h  ) + . 25  (h  <h  ) 
wl6  w20  w34  w35 


. 020833  (h  +h  +h  +h  ) + . 0625  (h  +h  +h  +h  ) 
wl  w8  wl2  wl5  w9  wlO  wl3  wl4 


+ . 08333  (h  +h  ) + .25  (h  +h  ) 
wll  wlO  'w32  w33 


A H 
CUSH  CUSH  n 


Number  of  Cushion  Compartments 

Height  of  hull  points  over  water  (Ft) 

2 

Cushion  compartment  area  = 800  Ft 

Average  Height  above  Water  of  Cushion  Compartment  (Ft) 

3 

Air  Volume  of  Cushion  Comp-xrtment  (Ft  ) 

Number  of  Hull  Bottom  Points 
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SUBROUTINE;  SKAKEA 


CUSHION  COMPARTMENT  ESCAPE  AREA 


CLR  = h - h >0 
k wk  s 


AREA  - 62.76  [CLR  + 4 CLR  + CLR  ] 

1 1 2 3 

H 94.  16  [CLR  + 3 CLR  + 3 CLR  .'CLR  ] 

3 4 5 6 

AREA^=  62.76  [CLR^^  + 4 CLR2.7  + CLR^j] 

I 94.  16  [CLR  + 3 CLR^  4 3 CLR,^  + CLR^/] 

6 24  25  26 

AREA  = 62.76  [CLR  + 4 CLR  + CLR  ] 

3 21  20  19 

+ 94.  16  [CLR  + 3 CLR  + 3 CLR  +CLR,^] 
19  18  17  15 

AREA  " 62.76  [CLR  4 4 CLR  + CLR  1 

4 12  11  1 

T 94.  16  [CLR  + CLR  -f-SCLR  4 CLR  J 
12  13  14  15 


k = 1 — 27  Number  of  Cushion  Periphery  Points 

h HEIGHT  OF  HULL  POINTS  OVER  WATER  (FT) 

wk 

h HEIGHT  OF  SKIRT  = 4.  5 FT 

s 

CLR  AIR  GAP  BETWEEN  SKIRT  AND  WATER  OF  HULL  POINTS  (FT) 
k 

AREA^  CUSHION  COMPARTMENT  1 ESCAPE  AREA  (FT^)  .Constant 

AREA^  CUSHION  COMPARTMENT  2 ESCAPE  AREA  (FT^)  .Constant 

AREA  CUSHION  COMPARTMENT  3 ESCAPE  AREA  (FT^)  . Constant 
3 

AREA^  CUSHION  COMPARTMENT  4 ESCAPE  AREA  (FT^)  .Constant 


! 

i 


1 
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TABLE  12  . PLAI'JFORM  POINT  LOCATIONS  (32) 


X 

Y 

Z 

k 

HULL 

HULL 

IIULI 

1 

10 

-18 

12 

2 

10 

-8 

12 

3 

10 

2 

12 

4 

-3  1/3 

2 

12 

5 

-16  2/3 

2 

12 

6 

-30 

2 

12 

7 

-30 

-8 

12 

8 

-30 

-18 

12 

9 

-16  2/3 

-18 

12 

10 

-3  1/3 

-18 

12 

11 

10 

-28 

12 

12 

10 

-38 

12 

13 

-3  1/3 

-38 

12 

14 

-16  2/3 

-38 

12 

15 

-30 

-38 

12 

16 

-30 

-28 

12 

17 

-43  1/3 

-38 

12 

18 

-56  2/3 

-38 

12 

19 

-70 

-38 

12 

20 

-70 

-28 

12 

21 

-70 

-18 

12 

22 

-56  2/3 

-18 

12 

23 

-43  1/3 

-18 

12 
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TABLE  12.  PLANl'ORM  POINT 

LOCATIONS 

(35) (CONI' U) 

k 

X Y 

HULL 

HULL 

Z 

HULL 

24 

-43  1/3 

2 

12 

25 

-56  2/3 

2 

12 

26 

-70 

2 

12 

27 

-70 

-o 

12 

28 

-56  2/3 

-8 

12 

29 

-43  1/3 

-8 

12 

30 

-16  2/3 

-8 

12 

31 

-3  1/3 

-8 

12 

32 

-3  1/3 

-28 

12 

33 

-16  2/3 

-28 

12 

34 

-43  1/3 

-28 

12 

35 

-56  2/3 

-28 

12 

fFT) 
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ACV 


PROGRAM  DESCRIPTION 


PROGRAM  NAME: 
PROGRAM  TYPE: 


INITIATED  BY: 


PURPOSE: 


CALLING  SEQUENCE: 
DESCRIPTION: 


(7/30/75) 


ATD 


Analog  Input  Driver 
Program  ACVSIM 

To  read  and  format  analog  input  date 
BAL,  15  ATD 


ATD  1 
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PROGRAM  NAME: 


PROGRAM  TYPE: 


INITIATED  BY: 


PURPOSE: 


CALLING  SEQUENCE: 


MVTRAEQUIl’CE.N  75-C-0057-  1 
ACV 

PROGRAM  DESCRIPTION 
(7/30/75) 

BACKGRND 
Background  Timer 
Program  PRINT 

To  provide  timing  for  the  background  and 
indirectly  the  real  time  simulation. 

B BACKGRND 


DESCRIPTION: 

The  program  BACKGRND  will  determine  how  long  the  background  idle 
loop  is  in  operation.  This  value  is  used  to  determine  the  execution  time  of  the  real 
time  programs.  BACKGRND  increments  a counter  which  is  scaled  at  IZ  psec  per 
unit.  The  timing  is  dependent  on  the  instructions  used  in  the  increment  loop.  The 
timing  is  done  only  at  the  option  of  the  operator.  The  time  value  is  stored  in  location 
BGTIMSAV  by  the  real  time  interrupt,  EXECACV.  If  the  timing  is  stopped  by  the 
operator  BACKGRND  is  terminated  by  branching  to  the  program  CRT. 


NAVTRAEOUll'CEN  75-C-0057-1 
ACV 

PROGRAM  DESCRIPTION 
(7/30/75) 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 


BCD AT A 

I/O  Subroutine 

EBCDIC  Input  Routines 

Translate  EBCDIC  character  strings  into 
one  word,  numeric  values. 


CALLING  SEQUENCE:  LW,6  BA(STRTADDR) 

LW,7  BA(ENDADDR) 

BAL,  15  BCDATA 
STW,8  WORD 
STW,5  MULT 

DESCRIPTION; 

The  subroutine  BCDATA  will  translate  a specified  byte  field  from 
EBCDIC  character  to  a numeric  word  value.  The  values  STRTADDR  and  ENDADDR 
are  the  starting  and  ending  byte  addresses  of  the  field  to  be  scanned.  The  output 
word  value  (WORD)  will  be  in  register  8 when  the  caller  resumes.  A value  called 
the  "multiple"  (MULT)  will  be  in  register  5.  This  word  will  be  explained  below. 

The  contents  of  register  9 will  be  destroyed. 

The  input  format  can  be  one  of  the  following: 

MMM  X'HHHHHH'  Hexadecimal  integer 

MMM  D'jpDDDDD'  Decimal  integer 

MMM  D'+ . DDDDDD'  Decimal  fraction 

MMM  D^D.DDDDDE+DD'  Floating  point,  single  value 

The  value  MMM  is  optional.  If  included,  it  is  translated  as  a decimal 
integer  and  is  returned  in  register  5,  This  value  is  used  to  specify  how  many  times 
the  value  that  follows  it  will  be  used  by  the  calling  program.  It  is,  therefore  called 
a "multiple".  If  the  multiple  is  not  present,  the  value  1 will  be  placed  in  register  5. 


The  format  of  the  remainder  of  the  field  is  described  in  the  ACV  Console 
Operating  Procedures  and  in  the  program  listing  for  BCDATA.  They  will  not  be 
repeated  in  detail  here.  Suffice  it  to  say  that  the  hexi decimal  format  is  translated  to 
a hexidecimal  integer  value  right  adjusted.  The  decimal  integer  format  is  translated 
to  a decimal  integer  value  scaled  B31.  The  decimal  fraction  format  is  translated  to 


BCDATA  1 
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a decimal  fraction  scaled  BO.  The  floating  point  format  is  represented  as  a lens 
. power,  mantissa-exponent  value.  It  is  translated  to  a standard  hex  power,  single 
word,  floating  point  value. 

If  any  irregularity  in  the  format  is  found  that  cannot  be  translated,  a - 1 
is  placed  in  the  return  value  of  register  5.  No  value  is  resolved. 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 


NAVTRiXEQMlPCKN  ’j-{ 

ACV 

PROGRAM  DESCRIPTION 
(7/30/75) 

BCDSUBS 

Input  Utility  Subroutines 

The  Programs  TYPE  and  CARDS 

To  Decode  input  EBCDIC  Data 

LI. 6 BA  (STRTADDR) 

LI,  8 BA{ENDADDR) 

BAL.  15  BCDCOMMA 
STW,7  COMMADDR 


LI,  6 BA  (STRTADDR) 
LI,  7 BA(ENDADDR) 
BAL,  15  BC  DP  ARAM 
STW,  5 ADDRESS 


LI,  6 BA(STRTADDR) 
LI, 8 BA(ENDADDR) 
LI,  5 CODE 
BAL,  15  BCDSTORE 
STW,  5 ERROR 


NAVTRAEQUJ  TCKN  7 5-(  0057-  1 


DESCRIPTION: 

The  program  BCDSUBS  is  made  up  of  three  subroutines.  These  subroutines 
are  used  to  decode  EBCDIC  character  strings  from  the  typewriter  and  the  care!  reader. 

The  subroutine  BCDCOMMA  has  the  sole  function  of  locating  the  character 
in  a character  buffer.  The  caller  must  specify  the  starting  and  ending  byte  location 
of  the  buffer  (STRTADDR  and  ENDADDR,  respectively).  The  address  of  the  byte 
previous  to  the  first  comma  encountered,  when  scanning  from  left  to  right,  is  stored 
in  register  7 (COMMADDR).  This  specifies  a field  that  contains  valid  input  data. 

There  are  two  alternate  conditions  specified  by  the  subroutine.  If 
no  comma  is  found  but  non-blank  characters  are  found,  a-2  is  placed  in  register  7. 

This  indicates  the  final  field  in  the  buffer.  If  no  comma  and  no  non-blank  characters 
are  found,  a- 1 is  placed  in  register  7.  This  indicates  either  that  a continuation  line 
follows  or  that  an  error  has  occurred. 


The  subroutine  BCDPARAM  will  search  a specified  byte  field  for  a valid 
parameter  that  can  be  used  on  a PARAM  input  line.  The  byte  field  is  specified  by 
the  byte  starting  address  (STRTADDR)  and  the  byte  ending  address  (ENDADDR)  in 
register  6 and  7,  respectively.  If  a valid  parameter  is  found,  its  associated 
memory  v/ord  address  is  placed  in  register  5 (ADDRESS).  If  no  valid  parameter 
is  found,  register  5 will  contain  a-1. 

The  subroutine  BCDSTORE  translates  the  address  and  data  fields  from 
PARAM  and  MEMORY  input  byte  strings  and  stores  the  data  into  the  given  address. 

The  caller  can  specify  one  of  three  types  of  lines  to  translate.  This  value  (CODE)  is 
placed  in  register  5 by  the  caller.  A 1 specifies  a PARAM  card.  A 2 specified  a 
MEMORY  card.  A 3 specifies  a continuation  card  from  a previous  PARAM  or  MEMORY 
card.  The  byte  starting  address  (STRTADDR)  and  the  byte  ending  address  (ENDADDR) 
of  the  input  buffer  are  stored  in  register  6 and  8,  respectively-  When  the  caller  of 
BCDSTORE  is  resumed,  register  5 (ERROR)  will  be  either  0 or  -1.  AO  specifies  a 
good  input  line.  A-1  specifies  that  an  error  was  found  somewhere  in  the  input  line. 

The  location  of  the  error  in  the  line  is  not  specified. 

The  PARAM  and  MEMORY  lines  are  processed  by  decoding  the  first  field 
of  the  line.  This  should  contain  the  starting  address  into  which  the  proceeding  data 
is  to  be  stored.  The  address  field  must  be  on  the  same  line  that  the  keyword 
PARAM  or  MEMORY  is  on  and  must  be  followed  by  a comma.  The  data  fields  are 
then  decoded  and  stored  one  after  the  other.  The  ending  address  +1  of  the  last  data 
stored  becomes  the  starting  address  of  the  next  data.  If  a buffer  is  ended  wdth  a 
comma  a-continuation  line  is  assumed  to  follow.  The  flag  (BCDSTWFL.)  is  set  to 
notify  the  caller  of  this  fact.  If  the  caller  finds  a continuation  line  (register  5=3),  the 
final  address  of  the  last  call  will  be  used  to  store  the  new  data. 

If  a data  field  cannot  be  decoded,  the  remainder  of  the  input  buffer  is 
ignored  and  register  5 will  be  set  to  -1.  Any  stores  that  were  made  up  to  the  point 
of  the  error  will  stay  in  effect.  The  continuation  flag  ,BCDSTWFL,  will  remain  reset. 
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If  the  input  line  was  intended  to  have  a continuation  line,  that  line  will  be  c onsidcrcd 
an  error,  since  it  will  not  contain  a valid  keyword.  Of  course,  if  ihe  aadre.-s  fieki 
contained  an  error,  register  5 will  be  set  to  - 1 , but  no  memory  stores  wj]]  bo 
performed. 


The  program  will  first  look  for  the  keyword,  ACVINIT,  starting  on 
column  one  of  each  card  read.  Cards  are  bypassed  until  this  card  is  read.  When 
an  ACVINIT  card  is  read,  only  one  of  the  following  cards  arc  valid:  a PARAM  card, 
a MEMORY  card,  or  an  ENDACVIN  card.  Alsc,  at  this  point  system  sense  switch 
8 can  be  put  off  and  cards  will  still  be  read  until  an  ENDACVIN  card  is  read. 

No  cards  \vill  be  processed  after  this  card  is  read.  If  sense  switch  8 is  off  when 
a valid  ENDACVIN  card  is  read,  card  reading  will  stop. 

The  P-A.RAM  and  the  MEMORY  cards  allow  the  computer  memory  to  be 
altered.  The  format  of  these  cards  is  discussed  in  the  program  listing  and  in  the 
Console  Operating  Procedures.  The  details  of  their  formats  -will  not  be  repeated 
here.  These  cards  can  have  continuation  cards.  Continuation  is  indicated  by  a 
comma  as  the  last  character  on  the  card.  Only  if  none  of  the  three  valid  keywords  ar' 
read  and  a continuation  is  indicated  from  the  previous  card,  will  the  card  be  treated 
as  a continuation  card-  Therefore,  a card  ending  wdth  a comma  does  not  have  to 
ha've  a continuation  card  after  it. 

Any  kind  of  format  error  will  cause  the  error  word  and  the  remainder 
of  the  card  to  be  bypassed  and  an  error  message  to  be  printed  on  the  line  printer 
and  the  typewriter.  Card  processing  will  continue,  however.  That  portion  of  the 
card  before  the  error  was  found  will  be  processed  and  any  valid  memory  alteration 
that  could  be  done  will  be  done.  Continuation  cards,  after  an  error  card,  will  also 
be  considered  errors. 
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PROGRAM  NAME; 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 
DESCRIPTION: 


NAVTRAF.QL  LPCK:  7 >•  » -OO  - /- 1 
ACV 

PROGRAM  DESCRIPTION 
(7/30/7^) 

DISCRETE 

Discrete  Driver 

Program  ACVSIM 

To  process  the  discrete  I/O  bits 

BAL,  15  DISCRETE 


The  program  DISCRETE  reads  the  discrete  bank  used  by  the  ACV  system. 
The  program  will  not  return  to  the  caller  until  all  the  discrete  I/O  is  completed. 


% 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 
DESCRIPTION: 


NAVTRilEQflPCKN  /'S-C-OOS/-  1 

ACV 

PROGRAM  DESCRIPTION 
(7/30/75) 

DLOOPER 
CRT  Driver 
Program  CRT 

To  set  up  ACV  position  variables  for  the 
CRT  driver  program  DLOOP 

BAL,15  DLOOPER 


The  program  DLOOPER  converts  the  ACV  simulation  X,Y,Z  position 
values  and  the  sines  and  cosines  of  the  pitch,  roll,  and  heading  from  floating  point 
values  to  integers  and  fractions.  These  conversions  are  made  for  a call  to  the  CRT 
driver  DLOOP.  The  X,Y,Z  position  values  are  converted  into  integer  inches.  The 
sines  and  cosines  are  converted  into  fractional  values  scaled  Bl.  Then  the  program 
DLOOP  is  called  to  drive  the  CRT. 


DLOOPER  1 
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PROGRAM  DESCRIPTION 

( 

(7/30/7H) 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 
DESCRIPTION: 


DTA 

Analog  Output  Driver 
Program  ACVSIM 

To  format  and  transmit  v'ariables  to  analog 
devices. 

BAL,  15  DTA 


The  program  DTA  drives  analog  output  devices  such  as  the  simulation 
cockpit  display  meters  and  the  chart  recorder.  DTA  takes  the  unformated  floating 
point  value  directly  from  simulation  parameters. 

There  are  two  types  of  analog  output  values.  Those  values  that  represent 
only  a voltage  and  those  that  represent  sines  and  cosines  to  a syncro  meter.  The 
syncro  needs  both  the  sine  and  cosine  of  its  associated  data  value.  Therefore,  these 
values  are  processed  twice. 

The  non  syncro  values  are  processed  as  follows.  They  are  first  adjusted 
to  subtract  any  bias  from  them.  Then  their  magnitude  is  limited  so  they  don't  go 
beyond  the  stops  of  their  associated  devices.  Finally  they  are  scaled  to  represent 
a fraction  of  10  volts  (10  volts  is  usually  full  scale  on  the  output  device).  The 
formatted  data  value  is  stored  in  an  analog  output  word  for  later  transmittal. 

Syncro  values  are  first  adjusted  for  any  bias.  They  are  then  limited  in 
magnitude  to  keep  their  value  in  meter  limits.  They  are  converted  to  degrees  and 
compared  to  their  last  value.  This  comparison  is  to  provide  a meter  change  limit 
so  that  the  meter  is  not  changed  to  a new  value  too  fast.  Finally  the  sine  or  cosine 
of  the  syncro  value  is  found  and  converted  from  floating  point  to  fixed  decimal  BU  . 
The  value  is  stored  in  an  analog  output  value  for  transmission. 

When  all  the  analog  output  values  are  converted,  they  are  all  sent  out 
to  the  analog  output  converter. 
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PROGRAM  DESCRIPTION 

* (7/30/75) 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 
DESCRIPTION: 


EXECACV 

Interrupt  Processor , Initialization 

Hardware  Interrupts,  Load  Procedure 

To  process  the  interrupts  and  to  initialize 
the  ACV  system 

(Interrupt  Linkage) 


The  program  EXECACV  contains  all  the  interrupt  processors  for  the 
ACV  system  and  also  contains  the  initialization  routine.  The  interrupts  that  are 
processed  are  the  Counter  3 Zero  Interrupt,  the  Control  Panel  Interrupts,  the 
Memory  Parity  Interrupt,  all  the  Error  Traps,  and  the  CAL4  Trap. 

The  Counter  3 Zero  Interrupt  provides  the  timing  for  the  ACV  simulation 
real  time  math  model.  It  maintains  the  internal  clock  used  to  print  out  the  time, 
it  processes  the  background  timer  counter  when  the  timing  of  the  background  is 
requested.  The  Counter  3 Zero  can  interrupt  itself.  If  it  does,  a record  of  this 
happening  is  kept  by  incrementing  a system  flag.  The  real  time  programs  are 
executed  with  Register  Pointer  set  to  0,  (Bits  55-59  of  the  PSD). 

The  ACV  simulation  is  initialized  at  load  time  and  any  other  time  upon  the 
request  of  the  operator.  During  initialization  the  interrupt  locations  are  reset,  various 
system  parameters  are  reset,  and  all  I/O  is  reset.  Location  X'  26'  is  set  to  a 
branch  to  initialization  routine.  This  adlows  the  operator  to  restart  the  system  by 
pressing  IDLE-CPU  RESET-RUN.  The  initialization  routine  is  terminated  by 
branching  to  background. 

The  Control  Panel  Interrupt  does  a logical  exclusive  or  to  a typewriter 
input  flag.  This  will  alternately  set  the  flag  to  zero  and  non-zero.  The  flag  is  tested 
by  the  typewriter  input  I/O  processor  to  allow  and  disallow  the  operator  to  perform 
input  type  commands. 

The  Power  Up  Interrupt  will  branch  to  ACV  initialization.  The  Power 
Down  Interrupt  will  execute  an  infinite  loop. 
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The  Memory  Parity  Interrcq^t  and  all  the:  Error  Traps  ar-  pr<.i.essod  in 
the  same  way.  They  execute  a CAL4  Trap.  CAl. 4 will  save  the  number  of  the  error, 
the  location  and  the  status  of  the  condition  code  at  the  time  oi  the  error,  the  status 
of  the  memory  fault  bits  (for  memory  parity),  and  the  16  registers  being  used.  The 
error  can  be  recovered  by  either  going  back  to  the  error  location  plus  one  or,  if 
the  console  switch  is  on,  going  to  the  initialization  routine. 
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PROGRAM  DESCRIPTION 
(7/30/70) 


PROGRAM  NAME: 


FIXPRNT 


PROGRAM  TYPE: 


Data  Formatter 


INITIATED  BY: 


Typewriter  and  line  printer  programs 


PURPOSE:  To  convert  floating  point  values  to  EBCDIC 

fixed  point  values 

CALLING  SEQUENCE:  LI,  2 BA{BUFFER) 

LW,  6 DATA 
BAL,  15  FIXPRNT 


DESCRIPTION: 

The  program  FIXPRNT  converts  a floating  point,  single  value  into 
EBCDIC  print  bytes.  The  format  of  the  output  is 

SDDDDDDDDD 

where  S is  the  sign  of  the  value,  +or  - . D is  either  a decimal  character  (0-9) 
or  a decimal  point.  The  position  of  the  decimal  point  depends  on  the  magnitude  of 
the  input  value.  If  the  value  is  smaller  than  tIO  -7  or  larger  than  _+  10  the 
data  is  converted  to  zero  or  + 999999999,  respectively.  Leading  zeros  are 
suppressed.  Ten  characters  are  always  stored  in  the  output  buffer  (BUFFER). 

The  calling  sequence  expects  register  2 to  contain  the  byte  address 
•of  the  buffer  in  which  to  store  the  output  bytes.  Register  6 should  contain  the 
data  value  to  be  converted.  When  the  calling  routine  is  resumed,  register  6 will 
be  unchanged  and  register  2 will  be  incremented  by  10. 

The  input  data  is  converted  from  a single  word,  hexidecimal  mantissa 
and  exponent  to  a double  word,  decimal  mantissa  and  exponent.  The  program 
FLTDEC  is  called  to  do  this.  The  mantissa  is  then  converted  into  EBCDIC 
characters  using  the  exponent  to  position  the  decimal  point. 
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PROGRAM  DESCRIPTION 
(7/30/75) 


PROGRAM  NAME: 


FLTDEC 


PROGRAM  TYPE: 


Data  Conversion  Subroutine 


INITIATED  BY: 


Caller 


PURPOSE: 


CALLING  SEQUENCE: 


DESCRIPTION: 


To  convert  a single  floating  point  value 
from  a hexadecimal  mantissa -exponent  to 
a decimal  mantissa-exponent 

LW,6  INDATA 
BAL,15  FLTDEC 
STW,6  OUTMANT 
STW.7  OUTEXP 


The  program  FLTDEC  converts  standard  single,floating  point  numbers 
from  base  16  to  base  10.  The  input  value  (INDATA)  should  be  in  register  6.  The 
output  mantissa  (OUTMANT)  and  the  output  exponent  (OUTEXP)  are  found  in  registers 
6 and  7, respectively.  The  converted  value  is  used  to  print  floating  point  values  onto 
the  line  printer  and  typewriter. 


The  method  used  is: 

16  **E16  = 10**E10 

Log  {16**E16)  = log  (10**E10) 

El6*log  16=E10 

16**E16  = 10**  (Ei6  log  16) 


E16**  log  16  is  broken  up  into  its  integer  and  fraction.  We'll  call  them 
lElO  and  FEIO,  respectively.  Therefore: 

• MANT  16*{16**E16)  = MANTl6*(10*=?FEl0)=t{10*^IEl0) 

Where  MANT  16  is  the  mantissa  of  the  input  floating  point  value.  lEIO 
becomes  the  output  exponent  to  the  base  10, 

MANT16*(10**FE10)  is  the  new  mantissa  for  the  decimal  floating  point 
number.  The  decimal  exponent  and  mantissa  may  get  adjusted  to  insure  that  the 
mantissa  is  between  +1  and  -1. 
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The  v'alue  oi  I0*^'FE10  is  calculated  with  a Iiastings  polvnoniial.  The- 
coefficients  are: 

a =1.0 
0 

a^  = 1.15129277 

a = .6627?o384 

2 

a^  = .254393574 

a = .072951736 

4 

a^  = .017421119 

5 

a = .002554917 

6 

a = . 000932642 

7 

The  polynomial  is: 

10*«X  = (I  a xV 

1 

(See  Appro.ximations  for  Digital  Computers  by  Cecil  Hastings,  Jr.  , 
Princeton  University  Press,  1955) 
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PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 


ACV 

PROGRAM  DESCKiPTlON 
(7/30/755 


FLTPRNT 


Data  Formatter 

T\T>ewritcr  and  line  printer  programs 


To  convert  floating  point  values  to 
EBCDIC  floating  point  values 

LI,  2 BA(BUFFER) 

LW,6  DATA 
BAL,  15  FLTPRNT 


DESCRIPTION: 

The  program  FLTPRNT  converts  floating  point,  single  word  data  values 
into  EBCDIC  print  bytes.  The  format  of  the  output  is; 

SD.DDDDDDD  ESDD 

where  S is  the  sign,  + or  -,  of  the  mantissa  and  the  e:<ponent.  The  D's  are  decimal 
characters  (0-9).  The  characters  '.  ’ and  'E'  are  constants.  The  exponent  represents 
a power  of  ten  exponent. 

The  calling  sequence  e:tpects  register  2 to  contain  the  byte  address  of  the 
buffer  in  which  the  output  characters  are  to  be  stored.  Register  6 contains  the  input 
floating  point  value.  The  calling  program  is  resumed  with  register  6 unchanged  and 
register  2 incremented  by  16. 

The  input  data  is  converted  by  the  program  FLTDEC  from  a hexi  decimal 
mantissa- exponent  to  a double  word,  decimal  mantissa  - exponent.  The  double  word 
is  then  converted  into  EBCDIC  characters. 
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NAVTRAEqri PCi;N  7 5-c-C)057- 1 
ACV 

PROGRAM  DESCRIPTION 

(7/ 3u/75 j 

INTEGER 

Data  Formatter 

Print  and  Typeout  Programs 

To  convert  computer  words  to  decimal  and 
hexi  decimal  integer  EBCDIC  formats. 

LI, 2 BA(BUFFER) 

LW,6  DATA 
BAL.15  DECIT 

LI,  2 BA(BUFFER) 

LW,6  DATA 
BAL,15  HEXIT 

DESCRIPTION: 

The  program  INTEGER  will  format  a data  word  into  EBCDIC  format.  It 
will  convert  the  value  into  EBCDIC  decimal  or  he.^  decimal  characters  depending  on 
whether  DECIT  is  called  or  HEXIT  is  called. 

The  word  to  be  converted  should  be  placed  in  register  6.  The  byte  address  I 

of  the  buffer  in  which  the  formated  bytes  are  to  be  stored,  should  be  placed  in  register  2,  * 

In  botn  cases  the  buffer  address  is  updated  by  the  number  of  bytes  used  in  the  formatting.  j 

will  always  use  eight  bytes  with  leading  zeros,  DECIT  suppresses  leading  zeros  j 

and  left  adjusts  the  characters.  However,  it  will  transfer  trailing  blanks  and  uses  11  1 

bytes.  } 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 


INTEGER  1 
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ACV 

PROGRAM  DESCRIPTION 
(7/30/75) 


PROGRAM  NAME: 
PROGRAM  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE: 


lOPERIF 

I/O  Subroutines  Drivers 


PRINT,  TYPE,  CARDS,  TYPEMSG 


To  provide  the  I/O  interface  between  the 
hardv'are  unit  and  its  software  processor. 

BAL,  15  CARDREAD 
DATA  BA(CARDBUF) 

BAL.  15  PRNTLINE 
DATA  BA(PRNTBUF) 

BAL,  15  PRNTPAGE 

BAL,  15  TYPEOUT 
DATA  BA{BUFFER) 

BAL,  15  TYPEIN 
DATA  DA(BUFFER ) 


DESCRIPTION: 


lOPERIF  contains  all  the  I/O  drivers  for  the  peripheral  equipment.  These 
subroutines  are  designed  to  be  used  in  background  only.  The  typewriter  routines 
can  be  used  by  the  typewriter  processors  only.  Each  subroutine  is  described  below. 


The  routine  CARDREAD  will  read  one  card  and  store  the  input  characters 
in  the  buffer  provided  by  the  caller.  The  program  will  always  store  80  characters. 

If  a card  reading  error  occurs,  the  program  will  wait  for  the  operator  to  correct  the 
error  on  the  reader  and  reread  a new  or  the  corrected  card. 


The  routine  PRNTLINE  will  print  one  line  of  data  onto  the  line  printer. 
The  input  buffer  should  have  as  its  first  character  the  size  of  the  remaining  buffer. 
This  first  character  is  not  printed.  The  input  buffer  can  be  any  size  but  only  132 
characters  or  less  will  be  printed.  When  the  calling  program  is  reentered,  it  is 
free  to  change  its  buffer  even  though  the  last  line  is  not  completed  printing. 


lOPERIF  1 
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The  routine  PHNTPAGE  will  skip  the  printer  page  to  the  head  of  paye. 

Nc  input  buffers  are  needed.  The  next  line  printed  will  be  put  on  the  first  line  of 
the  new  page. 

The  routine  TYPEOUT  \vill  type  one  line  of  data  onto  the  typewriter. 

The  first  character  of  the  input  buffer  should  be  equal  to  the  size  of  the  remaining 
buffer.  The  first  character  is  not  printed.  The  input  buffer  size  should  not  be 
longer  than  85  characters.  The  carriage  return  character  is  put  at  the  end  of  the  data 
line  by  the  TYPEOUT  subroutine.  TYPEOUT  returns  to  the  calling  program  before  the 
line  is  typed,  but  the  calling  buffer  is  immediately  available  for  reuse.  See  the  note 
below. 

The  subroutine  TYPEIN  reads  one  line  of  typed  data.  An  input  line  on  the 
typewriter  is  defined  as  one  that  ends  with  one  of  the  three  characters  NL(new  linej, 
EOM  (end  of  message),  or  HT  (horizontal  tab-T.A.B).  A line  that  exceeds  85  characters 
will  also  be  terminated  as  a complete  typed  input  line.  The  line  that  w'as  ended  by  an 
EOM  character  or  because  it  was  85  characters  long  w'^ill  be  rejected  by  TYPEIN  and 
TYPEIN  will  read  a new  line.  Also  if  the  line  was  terminated  with  an  I/O  error 
indication,  it  will  be  rejected  and  a new  line  will  be  read.  The  first  character  of  the 
resultant  input  line  will  be  equal  to  the  size  of  the  remaining  buffer  and  should  not  be 
considered  part  of  the  input  line. 

The  cent  sign  (^)  is  used  as  a backspace  by  the  operator.  The  input  line 
is  scanned  for  cent  signs  and  TYPEIN  will  correct  the  line  so  that  the  caller  will 
never  see  them.  If  the  ^ marks  backspace  the  entire  line,  a new  line  will  be  read. 

See  the  note  below. 

Note  - Since  the  typewriter  is  such  a slow  device  especially  on  input, 

•evice  busy  tests  use  the  routine  TYPEWAIT.  This  routine  will  allow'  other  devices 
besides  the  typewriter  to  be  referenced  even  though  the  typew'riter  is  in  the  middle 
of  an  I/O  request.  See  the  prograni  TYPE. 


lOPERIF  2 


26S 


NAVTRAEQUIPCHN  75-C-0057-1 
ACV 


PROGRAM  DESCI'iPTION 
{7/30/75} 


PROGRAM  NAME: 
PROGRAM  TYPE: 

T^TT*T»T  A 

X A U X 1 

PURPOSE: 

CALLING  SEQUENCE: 
DESCRIPTION: 


PRINT 


Background  I/O 


Program  CARDS 

To  print  the  state  parameters 

B PRINT 


SYSPRIN 


The  program  PRINT  is  a driver  for  the  program  SY<;pnT^” 
prints  the  simulation  state  parameters  Tf  T ® T . ^ 

not  be  done.  The  program  is  termin  k u fading  is  being  done, printing  uill 
program  ,s  termmatod  by  branching  to  the  program  BACKGRKD. 


PRINT  1 


NAVTRAEQUIPCEK  75-C-00f>/-i 
ACV 


PROGRAM  DESCRIPTiON 


PROGRAM  NAME: 


PROGRAM  TYPE: 


INITIATED  BY: 


PURPOSE: 


CALLING  SEQUENCE: 


DESCRIPTION: 


(7/30/75) 


PRNTTIME 


I/O  Subroutine 


SYSPRINT 


To  pro\'ide  the  time  and  date  in  EBCDIC  format 

BAL.  15  PRNTTIME 
DATA  BUFFER 


The  program  PRNTTIME  will  store  the  EBCDIC  time  and  date  into  the 
specified  buffer  (BUFFER  in  the  calling  sequence).  The  format  of  the  time  and 
date  is: 

HH:MM  MON  DD,  YY 

HH  is  the  military  hours 

MM  is  the  minutes 

MON  is  the  3 letter  month 

DD  is  the  day  of  the  month 

YY  is  the  last  two  numbers  of  the  year 

The  data  is  stored  on  a word  boundary  and  :ieeds  four  words. 

The  time  is  calculated  from  the  \’alue  RUNTIME  which  is  updated  by 
EXECACV.  RUNTIME  is  divided  by  RUNFREQ  which  equals  the  driving  frequency 
of  the  real  time  simulation.  This  calculates  the  elapsed  time  in  seconds.  The  minutes 
anj  hours  are  easily  derived  from  the  elapsed  seconds. 

The  hour  is  calculated  to  be  between  0 and  27  If  it  gets  to  be  24  or  larger, 
E4  hoars  is  subtracted  from  the  time  and  the  day  of  the  month  is  incremented.  The 
month  and  the  year  are  never  changed.  The  day  will  be  incremented  past  the 
proper  number  of  days  in  the  month  up  to  a maximum  of  99- 
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ACV 

PROGRAM  DESCRIPTION 
(7/30/75j 


P'"  . NAME; 

..\i  TYPE: 
INITIATED  BY: 
PURPOSE: 

CALLING  SEQUENCE; 


TYPE 

Background  I/O 
Program  CRT 

To  read  and  process  the  typewriter  input 
commands 

B TYPE 


DESCRIPTION: 

The  program  TYPE  will  read  and  process  typewriter  input  lines.  The 
operator  initiates  at',  input  line  by  putting  on  the  typewriter  input  light.  This  is  done 
with  the  Control  Panel  Interrupt.  The  interrupt  alternately  sets  and  resets  the  flag 
TYPEREQ.  When  the  flag  is  non-zero  (set),  a read  command  is  made  to  the  type- 
writer control  unit.  It  is  this  read  command  that  causes  the  typewriter  light  to  be 
put  on.  The  operator  then  types  a line  onto  the  typewriter,  "'he  program  will 
process  it  and  perform  a new  read  command.  When  the  flag  is  zero  (reset),  the 
pending  read  conamand  is  cancelled  and  the  light  is  turned  off. 

Because  the  typewriter  is  a slow  input  and  output  device,  a routine  has 
been  implemented  to  allow  other  background  tasks  to  be  performed  while  the  type- 
writer processors  are  waiting  for  the  typewriter  to  become  "not  busy".  This 
routine  is  TYPEWAIT.  The  routine  is  called  during  I/O  operations  by  either  of  the 
routines  TYPEOUT  or  TYPEIN.  Both  of  these  routines  are  in  the  program  lOPERlF. 
TYPEWAIT  saves  all  of  the  16  registers  and  branches  to  the  program  CARDS.  This 
allows  other  background  I/O  operations  to  be  performed.  Each  time  TYPE  is  entered 
a test  is  made  to  determine  if  TYPEWAIT  was  called.  If  it  was,  the  routine  that  called 
TYPEWAIT  is  resumed  with  all  of  its  registers  intact. 

The  program  TYPE  is  entered  from  the  program  CRT.  It  .-ormally  exits 
by  branching  to  the  program  TYPEMSG.  If  TYPEWAIT  is  in  process,  TYPE  will 
branch  to  CARDS. 

If  TYPEWAIT  is  not  in  effect,  TYPE  will  follow  its  normal  path  as  follows. 
I^  the  ACV  system  was  initialized,  TYtL'  will  first  type  a message  to  request  the 
operator  to  enter  the  date  and  time.  If  a read  request  is  indicated,  a line  is  read  and  it 
is  te.  ed  for  a valid  keyword.  The  keyword,  which  is  the  first  word  on  a typed  line, 
specifies  the  function  to  be  performea.  If  a valid  keyword  is  read, the  line  is  processed 
by  the  koyword  function  processor.  An  invalid  keyword  '-hich  is  not  a continuation  line 

TYPE  1 
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(see  MEMORY  and  PARAM)  will  cause  an  error  line  to  be  printed.  If  the  function 
requested  requires  type-outs,  the  output  operation  is  done  in  the  program  TYPE, 
not  in  TYPEMSG.  When  the  line  is  completely  processed  or  rejected  as  an  error 
a new  line  is  read  until  the  operator  resets  the  type  request  flag,  TYPEREQ. 

The  TYPE  program  can  recognize  five  keywords.  They  are  TIME,  DATE, 
MEMORY , PARAM,  and  DUMP.  The  format  and  function  of  each  keyword  request 
is  in  the  program  listing  and  the  Console  Operating  Procedures.  They  will  not  be 
repeated  here.  MEMORY  and  PARAM  are  duplicates  of  the  card  reader  inputs  of 
the  same  names.  They  can  have  continuation  lines.  None  of  the  other  requests  can. 
The  time  request  will  convert  the  input  time  into  units  of  the  real  time  update 
frequency.  The  time  in  seconds  is  multiplied  by  RUNFREQ  and  stored  in  RUNTIME. 
The  date  that  is  typed  in  with  the  DATE  keyword  is  stored  with  no  changes  made  to 
it.  The  DUMP  request  will  space  a line  and  type  one  location  per  line  until  all  the 
requested  locations  are  typed. 

If  any  of  the  keyword  processors  cannot  read  the  input  line,  an  error  line 
will  be  typed  and  the  line  will  be  ignored  from  the  point  of  the  error.  Except  for 
MEMORY  and  PARAM  this  means  the  whole  line  is  ignored.  The  error  is  not 
localized  by  the  software.  The  operator  has  to  determine  the  location  of  the  error 
in  the  typed  line. 


NAVTRAEQUIl'CEiN  75  {.-UOd/-  J 


ACV 

PROGRAM  DESCRIPTION 
(7/30/75i 


PROGRAM  NAME; 
PROGRAM  TYPE: 
INITIATED  BY; 
PURPOSE; 

CALLING  SEQUENCE; 
DESCRIPTION: 


TYPEMSG 
Background  I/O 
Program  TYPE 

Typeout  error  messages  onto  the  typewriter 
B TYPEMSG 


The  program  TYPEMSG  tests  positions  on  a message  queue.  If  any 
message  is  flagged  to  be  typed,  it  is  typed.  Because  of  the  length  of  time  it 
takes  to  type  messages  the  queue  is  examined  repeatedly  until  there  are  no  messages 
left  to  type.  The  program  ends  by  branching  to  the  card  reader  I/O  processor 
CARDS. 


TYPEMSG  I 


' < y 'liin  >"<  Iiillillllliilllll  <U  nil  • llliiil  Iii|,<|i||ii|ll>>il||l||liiril  ||||I' : 


SECTION  IX 

SUGGESTIONS  FOR  CUXi  JNL'IN'.:  DEVELOPMENT 


As  noted  in  the  main  body  of  this  report,  the  Vehicle  Generated  Wave  (VGW) 
model  does  not  run  in  real-time.  A considera'ole  amount  of  work  needs  to 
go  into  simulation  development  to  enhance  real-time  execution.  As  simple  as 
the  calculations  might  be,  the  sheer  volume-  of  computation  is  responsible 
for  the  non  real-time  operation  of  the  VGW  package.  The  Suggestions  for 
Improvement  presented  in  this  section  also  include  other  areas  which  came 
to  light  while  writing  this  report.  All  suggestions,  however,  are  aimed  at 
making  the  VGW  package  run  in  real-time. 


MODEL  UPDATES 

Updating  the  "Mathematical  Model  of  an  Air  Cushion  Vehicle"  is  a priority 
task.  This  task  will  resolve  differences  in  the  modelling  by  Cummings,  et 
al,  and  data  which  has  been  implemented  at  NT  EC.  The  simulation  program 
which  was  delivered  and  is  described  in  this  report  reflects  the  model  as 
given  and  does  not  necessarily  represent,  in  all  areas  of  modelling,  the 
true  vehicle.  The  mathematical  model  was  developed  from  Bell's  preliminary 
design.  Now  that  JEFF-B  has  been  built,  full-scale  actual  data  should  exist 
to  properly  update  the  model. 


Some  of  the  areas  to  be  investigated  are  as  follows; 

air  flow  rates 
air  pressures 

skirt  gap  areas  - skirt  clearances 
natural  frequency  of  vehicle 
drag  coefficients 
engine  characteristics 
wind  forces 

vehicle  weights  and  inertias 
fan  characteristics 

The  model  represents  the  best  estimate  from  Bell's  engineering  experience 
and  now  real  data  should  be  available  to  update  the  pilot  trainer  to  re- 
present the  actual  vehicle. 

CLEANUP  AND  TIMING  IMPROVEMENTS 

Simulation  program  cleanup  includes  tasks  as  diverse  as  coding  and  model- 
ling errors  which  have  been  identified  at  both  CSDL  and  NAVTRaequtpcen.  Although, 
these  have  been  properly  implemented  at  NTEC,  a proper  program  update 
should  be  accomplished  so  that  program  listings  and  documentation  can 
remain  current  at  CSDL  and  NAVTRAEOUIPCEN . 

There  are  several  areas  where  execution  Timing  Improvement  can  be  made. 
Perhaps  other  areas  exist,  but  it  would  necessitate  a close  look  by  an 
experienced  assembly  language  programmer.  The  Timing  Improvements 
described  here  are  exclusive  of  the  VGW  package,  which  is  described  in  the 
next  section. 
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A time  savings  can  be  realized  by  changing  'in-  model  of  tlio  gas  turbines 
Combining  the  three  port/star'noard  turbiii<-.s  uro  one-  engine  v, ill  save  considerable 
computation.  This,  of  course,  is  not  reali.siic  oui  there  is  o?)]\  one  pilot  con- 
trol each  for  port  and  starboard  simulaK'd 

A small  execution  time  savings  can  be  realizes!  using  a const;;int  wind  speed 
and  set  per  run.  This  would  necessitate  changjng  these  piirameters  du'dng 
initialization  only.  The  execution  time  saved  would  be  that  of  computing  the 
sine  and  cosine  functions. 

The  seaway  simulation  timi  ng  can  be  impro\  ed  hy  performing  algebraic  mani- 
pulation of  the  equations  and  by  using  the  piecewise  approximation  to  the 
cosine  function. 

There  is  a considerable  potential  of  execution  time  savings  by  improving  the 
cushion  pressure  solution  of  the  6 simultaneous  equations.  The  present  method 
uses  a Newton-Raphson  technique  to  compute  a 6 x 6 matrix  and  solution  by 
matrix  inversion.  To  create  the  matrix,  many  computation.s  are  involved  and 
32  function  calls  to  7ASQRT  and  7DASQRT  per  iteration.  .Although  this  i.s  an 
excellent  method  in  accuracy  and  convergence,  another  method  may  take  many 
more  iterations,  but  require  less  computational  time. 

A completely  different  method  for  solving  for  the  cushion  pressures  would 
completely  relieve  the  execution  duty  c^cle  by  implementing  the  6 simultaneoi’.s 
equations  on  the  analog  computer.  This  would  be  a fairly  simple  implementa- 
tion in  the  hybrid  facility  at  NTEC  and  should  guarantee  convergence. 

The  last  method  that  could  be  implemented  would  be  the  complete  programming 
in  fixed  point  arithmetic  of  the  all  digital  simulation.  This  would  be  a la.st 
resort  item  after  all  other  possibilities  have  been  explored.  There  is  great 
flexibility  now  by  using  floating  point  arithmetic,  and  fixed  point  arithmetic 
would  limit  that  flexibility.  Also,  this  would  necessitate  a nearly  complete 
redevelopment  of  the  real-time  simulation  program. 


VGW  TIMING 

The  Vehicle  Generated  Wave  (VGW)  model  presents  the  most  .significant  problem 
to  real-time  operation.  The  huge  number  of  calculations  to  be  performed  pre- 
sents the  greatest  difficulty.  At  the  present  time,  a complete  execution  to 
calculate  25  VGW  heights  takes  approximately  1.5  seconds.  For  each  wave 
height,  82  kernel  values  must  be  calculated  by  table  lookup  and  interpolated. 

Several  possibilities  exist  to  decrease  the  execution  time.  The  first  and 
foremost  must  be  to  eliminate  the  linear  interpolation  to  obtain  the  kernal 
value.  This  would  necessitate  taking  the  kernel  table  that  now  exists  and 
updating  it  with  constant  intervals  and  by  minimizing  kernel  error. 

Since  the  kernel  values  are  already  fixed  point  in  the  table,  the  floating 
conversion  could  be  eliminated  for  the  82  kernel  values  per  point  and  then 
the  convolution  integral  could  be  calculated  in  fixed  point  arithmetic.  By  doing 
this,  there  would  be  only  a total  of  25  floating  point  conversions  instead  of 
25  X 82  (2050), 
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After  all  modelling  and  timing  updates  have  been  performed,  a final  deter- 
mination needs  to  be  made  as  to  whether  the  VGW  model  can  be  run  in  real- 
time.  The  object  is  to  remove  the  first  order  lag  on  VGW  heights  (subroutine; 
FILTER).  By  this  time,  there  may  be  sufficient  execution  time  remaining  in 
the  duty  cycle  to  include  the  VGW  ''alculations.  Perhaps  not  all  25  points, 
but  5 points  a cycle  would  probably  be  fast  enough  for  true  vehicle  updates. 
This  would  yield  a total  update  rate  of  the  VGW  heights  of  250  msec  f4HZ). 
Studies  would  then  have  to  be  made  to  determine  the  adequacy  of  the  update 
rate. 

The  final  method  of  solution  would  be  to  off-load  all  VGW  height  calculations 
to  a fast  mini -computer  for  execution.  Then  the  only  load  on  the  simulation 
program  duty  cycle  would  be  the  data  transfer  between  computers,  which 
wovild  not  be  trivial. 


PILOT  VIRV.il  .Av^.KCTS  AND  vi'.tION  SI>!ULATI0KS 

For  a pilot  trainer  to  be  truly  effective,  proper  visual  affects  should  be 
simulated  for  the  pilot's  awareness.  Affixing  the  pilot  station  to  a motion 
platform  is  also  an  added  affect  on  the  pilot's  awareness.  Concurrent  deve- 
lopment of  the  hardware  ./software  interface  may  be  of  benefit.  At  the  very 
least,  program  listings  and  documentation  would"  be  jointly  held  and  current 


PEPORT  UPDATES 

In  no  case  should  changes  be  made  to  the  real-time  simulation  of  JEFF-B 
without  adequate  documentation.  This  documentation  should  not  consist  of 
merely  handwritten  notes  or  markups  within  this  report.  This  report  should 
be  properly  updated  at  proper  intervals  during  simulation  development. 
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